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Objectives: 

A  program  of  research  was  conducted  during  which  modeling  of  the  plasma  kinetics  of 
02('A)  generation  in  flowing  plasmas  for  use  in  electrically  excited  chemical  oxygen  iodine 
lasers  (eCOILs)  was  investigated.  The  goals  of  this  program  were  to  improve  the  knowledge 
base  in  this  important  technology  area,  provide  computational  support  to  ongoing  experiments 
and  to  propose  new  plasma  excitation  strategies  to  optimize  and  scale  02(!A)  production  and 
yields. 

Status  of  Effort: 

The  objectives  of  the  research  program  were  met.  Significant  progress  was  made  in 
improving  our  understanding  of  the  plasma  kinetics  of  O^'A)  production  for  pumping  of 
eCOILs.  Advances  were  made  in  3  areas:  1)  pulse  power  schemes,  2)  use  of  additives  and  3) 
high  pressure  operation.  Methods  and  scaling  laws  were  developed  to  optimize  the  yield  of 
O^'A)  in  flowing  plasmas  in  each  of  these  areas. 

Accomplishments  and  New  Findings: 

Computational  Issues: 

In  this  computational  investigation,  two  modeling  platforms  were  utilized.  nonPDPSIM 
is  a  two  dimensional,  plasma  hydrodynamics  model  in  which  the  fundamental  plasma  kinetics 
equations  are  solved  in  concert  with  the  Navier-Stokes  equations  for  flow,  plasma  chemistry  and 
electromagnetics  (or  electrostatics).  GlobalKin  is  a  global  or  plug  flow  model  which  uses  the 
same  reaction  mechanisms  as  nonPDPSIM.  Given  its  less-dimensional  nature,  Global  Kin  runs 
more  rapidly.  Significant  computational  improvements  were  made  to  nonPDPSIM  to  enable 
simulation  of  spiker-sustainer  circuitry,  high  pressure  flow  at  high  power,  and  having  multiple 
nozzles  for  injection  of  gases.  Similar  improvements  were  made  to  Global  Kin. 

Spiker  Sustainer  Circuitry: 

Operation  of  the  chemical  oxygen-iodine  laser  (COIL)  is  based  on  an  electronic  transition 
between  the  spin-orbit  levels  of  the  ground  state  configuration  of  the  iodine  atom  I(2Pi/2)  — > 
I(2P3/2)  (1.315  pm).  The  upper  level  is  populated  by  near  resonant  energy  transfer  (quantum 
defect  ~  219  cm'1)  from  02(‘A)  to  ground  state  I(2P3/2)-  The  COIL  generates  ground  state  I  by 
dissociative  excitation  transfer  from  O^'A)  to  I2.  Production  of  O^'A)  using  electric  discharges, 
eCOIL,  provides  a  means  for  an  all  gas  phase  system.  Much  of  the  development  of  eCOIL  has 
focused  on  efficiently  generating  02('A)  by  engineering  the  operating  E/N  (electric  field/gas 
number  density)  of  the  discharge  to  be  closer  to  the  optimum  value  for  exciting  02('A).  To 
maximize  the  fraction  of  discharge  power  that  is  directly  dissipated  in  the  electron-impact 
excitation  of  02('A),  the  electron  temperature,  Te,  should  be  near  1.2  eV,  which  corresponds  to 
an  E/N  of  *10  Townsend  [1  Townsend  (Td)  =  10'17  V-cm2].  Self-sustained  discharges  in  He/02 
mixtures  operate  at  least  a  few  tens  of  Td. 

One  of  the  methods  to  reduce  the  time  averaged  electron  temperature  is  to  use  a  pulsed 
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discharge  in  which  a  short,  high  power  pulse  (the  spiker)  is  followed  by  a  longer  period  of  lower 
power  (the  sustainer)  before  applying  another  high  power  pulse.  The  duration  of  the  high  power 
pulse  should  be  long  enough  for  Te  to  spike  and  for  the  gas  to  avalanche,  producing  an  electron 
density  in  excess  of  the  steady-state  value.  The  high  power  pulse  should  be  short  enough  so  that 
the  discharge  does  not  come  into  a  quasi-steady  state.  If  there  is  sufficient  excess  ionization 
during  the  period  of  lower  power  deposition,  there  may  be  an  extended  period  where  Te  falls 
below  the  self-sustaining  value,  Teo,  which  for  He/02  mixtures  enables  a  more  efficient 
production  of  Cty'A).  The  duty  cycle  of  the  spiker  and  length  of  the  sustainer  should  be  chosen 
so  that  the  average  value  of  electron  temperature  Te  <  Teo.  The  optimum  length  of  the  sustainer 

is  then  largely  determined  by  the  time  required  for  the  electron  density  to  decrease  to  its  steady 
state  value,  thereby  increasing  Te  towards  Teo.  This  method  for  engineering  Te  is  often  referred  to 
as  spiker-sustainer  (S-S)  excitation. 

In  previous  work  we  found  that  the  yield  of  (^('A)  generally  scales  linearly  with  energy 
deposition  for  moderate  loadings  up  to  5-8  eV/Ch-molecule.  Energy  deposition  beyond  these 
values  produces  excess  dissociation  which  ultimately  reduces  yield.  Initial  scaling  studies  of  S-S 
excitation  suggested  that  under  ideal  conditions,  as  are  represented  by  global  models,  (^('A) 
yields  approaching  30%  might  be  possible  using  S-S  methods.  In  this  work,  we  computationally 
investigated  S-S  methods  for  optimizing  Ch^A)  producing  using  nonPDPSIM.  Doing  so  more 
realistically  represents  the  electrical  circuitry,  electrode  losses,  uniformity  and  flow 
considerations.  In  particular,  we  investigated  capacitively  coupled  radio  frequency  (rf)  excited 
systems  in  He/02  mixtures. 

We  found  that  S-S  methods  do  hold  the  potential  for  increasing  Ch^A)  yields  above  that 
for  continuous  wave  (CW)  excitation  for  the  same  average  power,  though  the  amount  of  the 
increases  is  less  than  that  suggested  by  the  global  models.  S-S  techniques  generally  do  lower  the 
time  averaged  Te  compared  to  CW  excitation  and  so  improve  the  efficiency  of  excitation  and 
yield  for  C^'A).  This  advantage  to  S-S  excitation  is  diminished  at  higher  powers  where  Te  is 
naturally  lower  for  CW  excitation.  Lower  duty  cycles  (shorter  spiker  pulses)  are  generally  more 
advantageous  for  S-S  excitation  as  the  limit  of  a  delta-function  spiker  is  approached.  Increasing 
the  frequency  of  excitation  is  also  advantageous  due  to  the  lowering  of  Te  and  increasing  the 
electron  density,  ne.  Again,  the  advantage  of  S-S  excitation  is  diminished  as  frequency  increases 
due  to  the  natural  decrease  in  Te  with  CW  excitation.  The  length  of  the  sustainer  pulse  should  be 
short  enough  so  that  Te  is  largely  below  the  self  sustaining  value  but  long  enough  so  that  Te  and 
ne  begin  recovering  towards  their  steady-state  values.  This  recovery  indicates  a  good  utilization 
of  the  excess  ionization  produced  by  the  spiker. 

For  details,  see  the  Appendix  for  the  reprint:  Natalia  Yu.  Babaeva,  R.  A.  Arakoni  and  M. 
J.  Kushner,  "Production  of  (^('A)  in  Flowing  Plasmas  Using  Spiker-Sustainer  Excitation",  J. 
Appl.  Phys.  99,  113306  (2006). 

Use  of  Additives: 

In  eCOIL  systems,  (M’A)  is  typically  generated  upstream  of  the  laser  cavity,  h  is 
injected  immediately  prior  to  the  cavity  upon  which  the  flow  is  supersonically  expanded  to  lower 
the  gas  temperature  as  required  to  maximize  the  gain.  A  challenge  in  eCOILs  is  to  produce 
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sufficiently  high  yields  of  (^('A),  and  hence  laser  gain,  to  enable  the  laser  transition  to  be 
saturated  and  produce  high  power.  In  this  regard,  recent  research  has  focused  on  tailoring  the 
discharge  parameters  and  using  additives  such  as  CO,  H2,  D2  and  NO  to  improve  the  excitation 
efficiency  of  02^).  NO  and  NO2  are  also  used  as  additives  to  control  the  post-discharge 
chemistry.  In  fact,  all  demonstrations  of  laser  gain  and  oscillation  to  date  have  used  flowing 
plasmas  in  He/02  mixtures  with  NO  as  an  additive. 

eCOILs  differ  from  conventional  COILs  in  that  atomic  oxygen  is  also  produced  in  the 
electrical  discharge  by  electron  impact  dissociation  of  O2.  The  O  atoms  flow  downstream  where 
3-body  reactions  produce  O3  and  the  remaining  O  atoms  may  react  with  the  injected  h. 
Beneficial  reactions  of  O  with  I2  produce  ground  state  I(2P3/2)  atoms  (referred  to  as  I),  thereby 
eliminating  the  expense  of  02^)  molecules  for  the  initiating  dissociating  reactions.  Detrimental 
reactions  involving  O  atoms  include  quenching  of  I(2Pi/2)  (the  upper  laser  level,  referred  to  as  I  ) 
which  reduces  gain.  As  such,  management  of  the  O  atom  density  is  important  to  optimizing 
these  opposing  effects. 

Injection  of  NO  or  NO2  through  or  downstream  of  the  discharge  in  eCOIL  systems  has 
two  goals;  improving  the  efficiency  of  direct  production  of  02('A)  in  the  discharge  by  electron 
impact  and  management  of  the  O  atom  density  downstream  of  the  plasma.  Including  NO  in  the 
gas  stream  flowing  through  the  discharge  has,  in  part,  the  goal  of  improving  production  of 
02('A).  NO,  having  a  lower  ionization  potential  (9.26  eV)  than  either  O2  or  He,  is  likely  to 
provide  more  rapid  ionization  with  the  possibility  of  lowering  the  operating  E/N  (electric 
field/gas  number  density)  and  electron  temperature,  Te.  Lowering  Te  from  the  values  typical  of 
self-sustained  He/02  mixtures  is  advantageous  in  more  efficiently  producing  02('A)  by  direct 
electron  impact.  Discharges  in  NO  have  also  been  known  to  produce  02('A)  in  relatively  large 
amounts  even  though  O2  may  not  be  present  as  a  feedstock  gas.  However,  when  flowing  NO 
through  the  discharge,  some  of  the  power  that  would  otherwise  be  available  to  excite  O2  is 
dissipated  by  excitation  and  ionization  of  NO,  an  unwanted  consequence. 

In  addition  to  possibly  improving  the  production  of  02('A),  NO  and  NO2  are  potentially 
effective  in  managing  the  inventory  of  O  atoms  through  direct  and  cyclic  reactions  which  have 
the  effect  of  converting  O  atoms  back  into  O2.  In  the  context  eCOIL  systems  where  I2  is 
injected,  O  atoms  are  beneficial  by  dissociating  h  (and  product  species  10)  to  form  I  atoms 
which  are  then  pumped  to  I*  by  collisions  with  02('A).  The  O  atoms  are  detrimental  by 
quenching  I*.  Totally  eliminating  O  atoms  is  therefore  not  necessarily  beneficial.  NO  and  NO2 
also  have  secondary  effects  in  that  they  react  with  I  atoms  forming  intermediary  species  such  as 
INO,  and  INO2  which  further  react  with  I  to  reform  I2. 

A  computational  investigation  was  performed  on  the  consequences  of  NO  and  NO2 
additives  on  flowing  He/02  plasmas  and  their  afterglows  with  I2  injection.  These  investigations 
were  conducted  using  nonPDPSIM and  Global  Kin.  We  fund  that  although  for  most  conditions 
the  addition  of  NO  to  the  inlet  gas  stream  reduced  the  density  of  02('A),  it  ultimately  increased 
the  densities  of  I*  downstream  of  the  discharge  through  management  of  the  O  atom  density.  The 
inlet  NO  mole  fraction  also  typically  increased  the  extent  of  the  region  over  which  positive  laser 
gain  could  be  achieved.  This  can  be  particularly  useful  in  high  speed  flows  where  mixing 
lengths  are  longer.  Injection  of  NO2  in  the  post-discharge  flow  can  help  in  rapidly  scavenging  O 
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atoms  in  two  body  reactions  (as  compared  to  scavenging  by  NO  which  proceeds  by  a  3-body 
mechanism).  In  general,  addition  of  NO2  in  the  post-discharge  region  improves  laser  gain.  The 
downside  to  NO2  injection  is  a  rise  in  gas  temperature  due  to  the  exothermicity  of  the  reactions 
between  NO2  and  O. 

For  details,  see  the  Appendix  for  the  reprint:  R.  A.  Arakoni,  N.  Y.  Babaeva  and  M.  J. 
Kushner  "02('A)  Production  and  Gain  in  Plasma  Pumped  Oxygen-Iodine  Lasers:  Consequences 
of  NO  and  N02  Additives",  J.  Phys.  D  40,  4793  (2007). 

High  Pressure  Operation: 

Many  system  issues  motivate  operating  eCOILs  at  higher  pressures  to  obtain  larger 
densities  of  02('A)  for  a  given  yield  and  to  provide  higher  back  pressure  for  supersonic 
expansion  in  the  laser  cavity.  If  higher  order  effects  are  not  important,  and  parameters  such  as 
energy  deposition  per  molecule  are  maintained  constant,  it  is  expected  that  absolute  02(!A) 
production  (that  is,  the  total  number  of  (^('A)  molecules  produced)  should  scale  linearly  with 
pressure,  thereby  providing  additional  motivation  for  higher  pressure  operation. 

In  previous  work,  it  was  found  that  the  yield  of  (^('A)  generally  scaled  linearly  with 
energy  deposition  at  low  pressures  (<10  Torr)  until  the  ground  state  O2  is  depleted.  We 
computationally  investigated  radio  frequency  (rf)  discharge  excited  flowing  He/02  plasmas  using 
nonPDPSIM  and  Global  Kin  with  the  goal  of  developing  scaling  laws  for  02(  A)  production  in 
eCOIL  systems  when  operating  at  higher  pressures. 

We  found  that  the  densities  and  yields  of  02('A)  can  have  significant  sub-linear  scaling 
with  pressure.  Although  yields  may  decrease  with  increasing  pressure,  the  absolute  densities  of 
02('A)  typically  do  increase.  Although  these  results  depend  on  the  layout  of  the  electrodes  and 
the  aspect  ratio  of  the  flow  tube  more  general  conclusions  can  be  made.  Obtaining  high  yields  of 
02('A)  will  require  careful  management  of  the  O3  density.  Left  unchecked,  quenching  of  02('A) 
by  the  O3  produced  at  high  pressure  is  a  rate  limiting  step.  At  higher  energy  densities  and 
pressures,  quenching  by  O  atoms  and  energy  pooling  with  02('A)  may  also  become  important. 
This  is  particularly  the  case  for  possible  3-body  quenching  involving  O  atoms,  though  the 
importance  of  this  reaction  can  be  minimized  by  managing  the  density  of  O  atoms  with  additives 
such  as  NO.  Pressure  scaling  also  requires  management  of  the  gas  temperature,  as  exothermic 
recombination  reactions  rapidly  increase.  Although  the  destruction  of  O3  at  higher  temperatures 
is  beneficial,  in  general  an  intermediate  gas  temperature  is  likely  preferred. 

Discharge  stability  must  also  be  managed  and  this  will  be  highly  temperature  dependent. 
Non-uniform  excursions  of  gas  temperature  in  the  plasma  zone  will  likely  produce  non-uniform 
power  deposition  and  constrictions  near  the  electrodes.  As  such,  low  aspect  ratio,  transverse 
electrodes  allowing  for  more  aggressive  gas  cooling  and  more  distributed  electric  fields  will 
likely  be  required  for  high  pressure  operation 

For  details,  see  the  Appendix  for  the  reprint:  N.  Y.  Babaeva,  R.  A.  Arakoni  and  M.  J. 
Kushner  "(^('A)  Production  in  High  Pressure  flowing  He/02  Plasmas:  Scaling  and  Quenching", 
J.  Appl.  Phys.  101,  123306  (2007). 
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In  chemical  oxygen  iodine  lasers  (COILs),  oscillation  at  1.315 /im  in  atomic  iodine  (~^l/2 
-+2PVJ  is  produced  by  collisional  excitation  transfer  of  Oy^A)  to  I2  and  I.  Plasma  production  of 
02('A)  in  electrical  COILs  (eCOILs)  eliminates  liquid  phase  generators.  For  the  flowing  plasmas 
used  for  eCOILs  (He/02,  a  few  to  tens  of  torr),  self-sustaining  electron  temperatures,  Te%  are 
2-3  eV  whereas  excitation  of  02( 1  A)  optimizes  with  7,=  1-1.5  eV.  One  method  to  increase  02(*A) 
production  is  by  lowering  the  average  value  of  Te  using  spiker-sustainer  (SS)  excitation  where  a 
high  power  pulse  (spiker)  is  followed  by  a  lower  power  period  (suslainer).  Excess  ionization 
produced  by  the  spiker  enables  the  sustainer  to  operate  with  a  lower  Te.  Previous  investigations 
suggested  that  SS  techniques  can  significantly  raise  yields  of  02(‘A).  In  this  paper,  we  report  on  the 
results  from  a  two-dimensional  computational  investigation  of  radio  frequency  (if)  excited  flowing 
He/02  plasmas  with  emphasis  on  SS  excitation.  We  found  that  the  efficiency  of  SS  methods 
generally  increase  with  increasing  frequency  by  producing  a  higher  electron  density,  lower  Te,  and, 
as  a  consequence,  a  more  efficient  production  of  02('A).  ©  2006  American  Institute  of  Physics. 

[DOI:  10. 1063/1.2199387] 


I.  INTRODUCTION 

Chemical  oxygen  iodine  lasers  arc  being  investigated  be¬ 
cause  of  their  optical  fiber  deliverable  wavelength 
(1.315  /jl m),  highly  scalable  continuous  wave  (cw)  power, 
and  favorable  material  interaction  properties.16  Operation  of 
the  chemical  oxygen  iodine  laser  (COIL)  is  based  on  an  elec¬ 
tronic  transition  between  the  spin-orbit  levels  of  the  ground 
state  configuration  of  the  iodine  atom  I(2A>,/2)  — ►  I(2^>3/2) 
where  the  upper  level  is  populated  by  near  resonant  energy 
transfer  (quantum  defect  —219  cm-1)  from  02(jA)  to  ground 
state  l(2P3/2).  The  COIL  generates  ground  state  I  by  disso¬ 
ciative  excitation  transfer  from  O^A)  to  I2. 

Typically,  02(*A)  is  produced  in  an  external  chemical 
reactor  by  a  gas-liquid  reaction  between  gaseous  chlorine 
and  a  basic  hydrogen  peroxide  solution7  producing  yields 
approaching  100%  of  the  oxygen  emerging  in  the  02(!A) 
state.  The  long  lifetime  of  02(IA)  (60  min)  and  robustness 
against  quenching  enables  transport  over  long  distances  to 
the  laser  cavity.  There  are  many  system  issues  having  to  do 
with  weight,  safety,  and  the  ability  to  rapidly  modulate  the 
production  of  the  02(,A)  which  have  motivated  investiga¬ 
tions  into  methods  to  produce  the  precursor  02(!A)  using 
flowing  electric  discharges  and  so  produce  a  purely  electrical 
excited  laser  (eCOIL).80  Recent  and  ongoing  investigations 
have  shown  that  substantial  yields  of  CM1  A)  can  be  gener¬ 
ated  by  an  appropriately  tailored  electric  discharge,  typically 
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in  mixtures  with  rare  gas  diluents  such  as  He.l<>_22  Recently, 
positive  gain  was  reported  in  atomic  iodine  resulting  from 
electric  discharge  produced  0>( 1  A).23,24  followed  by  laser 

25 

demonstrations."' 

Much  of  the  development  of  eCOIL  has  focused  on  ef¬ 
ficiently  generating  O^C'A)  by  engineering  the  operating 
E/N  (electric  field/gas  number  density)  of  the  discharge  to  be 
closer  to  the  optimum  value  for  exciting  02(!A).  To  maxi¬ 
mize  the  fraction  of  discharge  power  that  is  directly  dissi¬ 
pated  in  the  electron-impact  excitation  of  CM1  A),  the  elec¬ 
tron  temperature,  Te .  should  be  near  1.2  eV.  which 
corresponds  to  an  E/N  of  ~  10  Townsend  [I  Townsend 
(Td)=  lO-17  V  cm2].  Self-sustained  discharges  in  Hc/02 
mixtures  operate  at  least  a  few  tens  of  Td.27 

One  of  the  methods  to  reduce  the  time  averaged  electron 
temperature  is  to  use  a  pulsed  discharge  analogous  to  that 
proposed  in  Refs.  12  and  17.  In  these  devices  a  short,  high 
power  pulse  (the  spiker)  is  followed  by  a  longer  period  of 
lower  power  (the  sustainer)  before  applying  another  high 
power  pulse.  The  duration  of  the  high  power  pulse  should  be 
long  enough  for  Tt  to  spike  and  for  the  gas  to  avalanche, 
producing  an  electron  density  in  excess  of  the  steady-state 
value.  The  high  power  pulse  should  be  short  enough  so  that 
the  discharge  does  not  come  into  a  quasi-stcady-statc.  If 
there  is  sufficient  excess  ionization  during  the  period  of 
lower  power  deposition,  there  may  be  an  extended  period 
where  Tr  falls  below  the  self-sustaining  value.  1\,0.  which  for 
He/02  mixtures  enables  a  more  efficient  production  of 
02(’A).  The  duty  cycle  of  the  spiker  and  the  length  of  the 
sustainer  should  be  chosen  so  that  the  average  value  of  elec¬ 
tron  temperature  fe<  Tto.  The  optimum  length  of  the  sus- 
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tainer  is  then  largely  determined  by  the  lime  required  for  the 
electron  density  to  decrease  to  its  steady-slate  value,  thereby 
increasing  Te  towards  Tto.  This  method  for  engineering  Te  is 
often  referred  to  as  spiker-sustainer  (SS)  excitation. 

In  previous  work,  the  scaling  of  production  of  02('A) 
was  computationally  investigated  using  global-kinetics19 
one-dimensional.20  and  two-dimensional  models.21  It  was 
found  that  the  yield  of  02(’A)  generally  scaled  linearly  with 
energy  deposition  for  moderate  loadings  (a  few  eV/02  mol¬ 
ecule)  up  to  5-8  eV/molecule.  Energy  deposition  beyond 
those  values  produced  excess  dissociation  which  ultimately 
reduced  yield.  Initial  scaling  studies  of  SS  excitation  meth¬ 
ods.  discussed  in  Ref.  20.  suggested  that  under  ideal  condi¬ 
tions,  as  are  represented  by  global  models,  02(,A)  yields 
approaching  30%  might  be  possible  using  SS  methods.  In 
this  work,  we  computationally  investigate  SS  methods  for 
optimizing  02('A)  using  a  two-dimensional  (2D)  model  that 
more  realistically  represents  the  electrical  circuitry,  electrode 
losses,  uniformity  and  flow  considerations.  In  particular,  we 
investigated  capacitively  coupled  radio  frequency  (rf)  ex¬ 
cited  systems  in  He/02  mixtures.  We  found  that  SS  methods 
do  hold  the  potential  for  increasing  02(,A)  yields  above  that 
for  cw  excitation  for  the  same  average  power,  though  the 
amount  of  the  increases  is  less  than  that  suggested  by  the 
global  models. 

The  model  and  reaction  mechanism  are  briefly  described 
in  Sec.  II.  The  results  from  our  investigation  are  discussed  in 
Sec.  Ill  where  yields  and  efficiency  of  O^A)  production  are 
compared  for  cw  and  SS  excitations  for  rf  carrier  frequencies 
of  13.  27,  and  40  MHz.  Optimization  of  02(*A)  production 
using  the  SS  technique  is  discussed  in  Sec.  IV.  Concluding 
remarks  are  in  Sec.  V. 

II.  DESCRIPTION  OF  THE  MODEL 


temperature  by  solving  the  electron  energy  equation  for  av¬ 
erage  energy  e. 


f)(nee) 

fit 


—  q\\,  •  E  -  N )Kj  -  V  • 

j 


(4) 


where  jkTr=t:.  The  terms  in  Eq.  (4)  are  for  the  contribution 
from  Joule  heating,  summation  of  elastic  and  inelastic  impact 
processes  with  heavy  neutrals  and  ions  with  energy  loss  Kr 
and  electron  heat  flux  consisting  of  terms  for  electron  energy 
flux  (T,)  and  a  conduction  (\r  is  the  electron  thermal  con¬ 
ductivity).  The  electron  transport  coefficients  and  rate  coef¬ 
ficients  for  bulk  electrons  as  a  function  of  Tr  arc  obtained  by 
solving  the  zero-dimensional  Boltzmann’s  equation  for  the 
electron  energy  distribution  to  capture  the  non-Maxwcllian 
nature  of  the  electron  swarm.  These  values  arc  stored  in  a 
tabular  form  and  interpolated  during  execution  of  the  code. 
The  tables  are  periodically  updated  to  reflect  changes  in  spe¬ 
cies  densities.  These  updates  are  then  followed,  in  a  time 
splicing  manner,  with  an  implicit  update  of  neutral  particle 
densities. 

The  fluid  averaged  advcctivc  velocity  v  is  obtained  by 
solving  a  modified  form  of  the  compressible  Navicr-Stokes 
equations  in  which  momentum  transfer  from  ion  and  electron 
collisions  and  acceleration  by  the  electric  field  arc  included 
in  the  momentum  equation,  and  Joule  heating  is  included  in 
the  energy  equations, 


r9p 

in 


=  -  V  •  (pv)  +  (inlets,  pumps). 


(5) 


(Xpv) 

fit 


=  -  Sip  -  V  (f>\\)  -  V  •  t+  X  E. 


The  model  used  in  this  study,  nonPDPSIM,  is  a  multi¬ 
fluid  2D  hydrodynamics  simulation  in  which  transport  equa¬ 
tions  for  all  charged  and  neutral  species  and  Poisson's  equa¬ 
tion  are  integrated  as  a  function  of  time.  nonPDPSIM  is 
described  in  detail  in  Ref.  28.  Poisson’s  equation  [Eq.  (1)], 
transport  equations  for  conservation  of  the  charged  species 
[Eq.  (2)].  and  the  surface  charge  balance  equation  [Eq.  (3)] 
arc  simultaneously  integrated  using  a  Newton  iteration  tech¬ 
nique, 

-  V  •  (e()e,  V  <I>)  =  2  Njcfj  +  pv,  ( I ) 

j 


rlN i 
dt 


=  - V 


r/  +  sr 


(2) 


T1  =  S  </,(-  V  ■  I,  +  5,)  -  V  •  M-  V4>)].  (3) 

Here  e«.  F.r,  pv.  Nr  F7,  o\  Sr  and  </,  arc  the  permittivity  of 
free  space,  dielectric  constant,  electric  potential,  surface 
charge  density,  conductivity  of  solid  materials,  sources,  and 
charge,  respectively.  The  subscripts  j  denote  gas-phase  spe¬ 
cies.  Updates  of  the  charged  particle  densities  and  electric 
potential  are  followed  by  an  implicit  update  of  the  electron 


(6) 

=  -  V(-  K  V  T+  pvcpT)  +  2  j,  •  E  -  2  K,A H, 

*  j 

+  /;V  v.  (7) 

Here  p  is  the  total  mass  density,  p  is  the  thermodynamic 
pressure,  r  is  the  viscosity  tensor,  cp  is  the  heat  capacity,  k  is 
the  species  averaged  thermal  conductivity,  fi  is  the  mobility, 
and  M  the  molecular  weight.  The  subscripts  j  are  for  sum¬ 
mations  over  species.  A H,  is  the  change  in  enthalpy  due  to 
reaction  /  having  total  rate  Rr  The  reactions  include  Frank- 
Condon  heating  from  electron-impact  dissociation  of  mol¬ 
ecules  as  well  as  conventional  chemical  reactions.  The  sums 
(other  than  for  reactions)  arc  over  all  charged  and  neutral 
species.  The  contributions  to  momentum  from  charged  par- 
tides  include  those  of  electrons.  The  contributions  to  the 
energy  equation  from  Joule  heating  include  contributions 
from  ions.  The  heal  transfer  from  electrons  is  included  as  a 
collisional  change  in  enthalpy.  The  relationship  among  pres¬ 
sure.  density,  and  temperature  is  given  by  the  ideal-gas  law 
The  contributions  to  momentum  due  to  charged  panicle 
transport  assume  that  the  collision  frequencies  with  electrons 
and  ions  arc  large  compared  to  the  time  rate  of  change  in  the 
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electric  field,  and  so  there  is  little  momentum  that  is  instan- 
taneously  stored  in  the  charged  particles." 

The  reaction  mechanism  lor  He/02  plasmas  used  here  is 
essentially  the  same  as  that  described  in  Refs.  19-21,  and 
involves  reactions  in  the  gas-phase  discharge  and  afterglow 
as  well  as  recombination  and  quenching  reactions  on  the 
discharge  tube  walls.  The  species  in  the  model  are  ground 
state  neutrals  02,  O,  O3,  and  He;  vibrationally  excited  02(i;) 
(representing  the  total  vibrational  population  consisting  of 
the  first  four  vibrational  levels  of  02);  electronic  states 
03('A).  CM'S),  0('D).  O('i’),  and  He(JS);  and  ions  02+, 
O*.  O,’,  O',  O,'.  and  He+. 

The  states  of  primary  interest,  the  precursors  to  the 
cCOIL.  are  02('A)  (0.97  cV)  and  02('X)  (1.6  eV).  which 
are  dominantly  produced  in  the  discharge  region  by  direct 
electron  impact  with  the  ground  state. 


e  +  02  -*02('A)  +  £\ 

(8) 

<?  +  o2  —  Oj('2)  +  *. 

(9) 

Production  of  generally  also  results  in  generation  of 

02( 1  A)  through  rapid  collisional  quenching  reactions  with 
atomic  oxygen  O.  The  atomic  oxygen  is  dominantly  pro¬ 
duced  by  electron-impact  dissociation  of  02, 


e  +  02->0  +  0  +  e. 

(10) 

0  +  02('X)-*0  +  02('A). 

(11) 

Note  that  the  02(,S)  and  02(,A)  states  have  radiative  life¬ 
times  of  12  s  and  approximately  I  h.  These  states  are  only 
nominally  quenched  by  collisions  with  species  other  than 
electrons  for  the  conditions  of  interest  although  at  higher 
pressures  quenching  of  02('A)  by  O3  becomes  important. 
These  electron  collisions  arc.  in  order  of  importance,  super- 
elastic  relaxation,  dissociative  excitation,  and  ionization. 
Other  channels  for  02(IS)  production  (although  less  impor¬ 
tant  for  the  conditions  of  interest)  are  through  production  and 
quenching  of  O('O). 


e  +  02-*0(lD)  +  0  +  e, 

(12) 

e  +  O  —  >  0( } D)  +  ey 

(13) 

0(]  D)  +  02—>  O  +  02(,S). 

(14) 

The  generation  of  02('i)  in  these  reactions  is  then  followed 
by  quenching  by  O  to  CM1  A).  As  the  density  of  02('A) 
increases  to  the  many  percent  level,  losses  to  upper  elec¬ 
tronic  states,  superelastic  deexcitation  to  the  ground  state, 
and  dissociation  begin  to  become  important. 


e  +  02(  *A)  — ►  02('X)  +  e. 

(15) 

e  +  Oit'A)  — »  02  +  e. 

(16) 

<?  +  02('A)  —  O  +  O +  e. 

(17) 

<*  +  02('A)  —  0('D)  +  0  +  e. 

(18) 

The  02( 1  A)  then  persists  far  into  the  afterglow  due  to  its  long 
radiative  lifetime,  where  the  most  significant  quenching 


mechanism  is  through  collisions  with  O  atoms, 

0  +  02('A)-.02  +  0.  (19) 

At  this  point  energy  pooling  and  quenching  reactions  with 
02(!A)  and  03  also  begin  to  reduce  02('A)  yields,  though 
not  significantly  for  our  conditions.  The  rate  of  quenching 
02(‘A)  by  collisions  with  the  walls  is  uncertain  due  to  the 
variability  of  the  quenching  probability  with  temperature  and 
conditions  of  the  wall.  Based  on  estimates  for  similar  condi¬ 
tions,  we  have  assigned  a  wall  quenching  coefficient  of  10"\ 
which  is  insignificant  for  typical  eCOIL  conditions. 

The  effective  yield  of  02(‘A)  is  defined  as  the  ratio  of 
the  combined  02(!A)  and  CM1^)  densities  to  the  sum  of  the 
densities  of  all  oxygen-containing  species  on  a  molecular  02 
basis, 

_ [02('A)  +  Q2('X)] _ 

(LOjJ  +  102C  A)J  +  |02('X)]  +  0.5[OJ  +  1 .510, J} ' 

(20) 

This  choice  of  yield  was  made  with  the  prior  knowledge  that 
the  majority  of  02('Z)  is  quenched  directly  to  02('A).  As 
such,  Eq.  (20)  is  the  most  relevant  for  the  best  case  energy 
scaling. 

III.  CONTINUOUS  WAVE  AND  SPIKER-SUSTAINER 
EXCITATIONS 

A  schematic  of  the  idealized  cCOIL  device  we  modeled 
is  shown  in  Fig.  1.  A  Hc/02  =  70/30  mixture  is  flowed 
through  a  quartz  tube  60  cm  in  length  and  6  cm  in  diameter 
at  3  Torr.  A  rf  electric  discharge  is  operated  between  two 
ring  electrodes  2  cm  wide  with  centers  separated  by  13  cm. 
The  electrodes  are  powered  up  to  a  few  hundred  watts  at 
13.56,  27,  or  40  MHz.  The  flow  rate  of  6  lpin  corresponds  to 
an  average  axial  inlet  speed  of  985  cm/s.  Our  investigations 
were  limited  to  the  region  of  the  reactor  prior  to  supersonic 
expansion  and  injection  of  I2.  The  numerical  grid  uses  an 
unstructured,  cylindrically  symmetric  mesh  with  triangular 
elements.  The  wall  temperature  was  held  fixed  at  300  K  as¬ 
suming  active  water-jacket  cooling  (not  included  in  the 
mesh). 

A.  Continuous  wave  excitation 

As  a  point  of  departure  and  to  provide  a  basis  of  com¬ 
parison.  02(,A)  production  in  a  flowing  afterglow  using  cw 
rf  capacitive  excitation  will  be  discussed.  The  base  case  op¬ 
erating  conditions  are  He/02  =  70/30.  3  Torr.  flow  rate  of 
6  1pm,  inlet  gas  temperature  of  300  K.  and  power  deposition 
of  40  W  with  a  rf  of  13.56  MHz. 

Power  deposition,  electron  temperature  (7*,).  and  elec¬ 
tron  density  (n,)  for  the  base  case  are  shown  in  Fig.  2(a). 
Values  are  shown  averaged  over  the  rf  cycle.  Power  deposi¬ 
tion  is  peaked  just  off  axis  with  a  maximum  of  0.2  W/cm  \ 
The  peak  value  of  T,  of  2.9  eV  is  near  the  upstream  electrode 
where  the  electron  density  is  low.  Tt  drops  to  2.2  eV  where 
the  electron  density  is  maximum  at  about  9.8  X  I0g  cm"3  and 
where  the  majority  of  the  CM1  A)  is  produced.  The  region  of 
elevated  Tr  extends  significantly  beyond  the  rf  electrodes 
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FIG.  1 .  Schematic  of  the  cylindrical  discharge  tube,  (a)  Close-up  showing 
ring  electrodes  and  (b)  complete  computational  domain  showing  the  numeri¬ 
cal  mesh.  The  flow  enters  uniformly  from  the  right  (or  top)  where  pressure, 
speed,  and  temperature  are  specified.  The  walls  are  held  at  constant  tem¬ 
perature  and  mass  flux  is  conserved  at  the  outlet.  Boundary  conditions  for 
solution  of  Poisson’s  equation  require  the  extension  of  computational  do¬ 
main  beyond  the  tube. 


|  Inlet 


(upstream  and  downstream)  due  to  the  large  electron  thermal 
conductivity  and  gas  flow  which  entrains  ions  through  their 
large  momentum  transfer.  Although  the  electrons  transfer 
little  momentum  to  or  from  the  gas,  the  electrons  are  pulled 
by  the  ions  in  the  downstream  direction  through  the  ambipo- 
lar  electric  field.  The  electron  density  is  as  large  as  I09  cm-3, 
a  few  centimeters  ahead  of  the  upstream  electrode  due  to 
thermal  conduction  and  diffusion  of  electrons  against  the  ad- 
vcctivc  flow  and  some  local  ionization.  The  steep  rise  in 
electron  density  at  this  point  marks  the  front  end  of  the  re¬ 
active  plasma  zone.  The  same  action  occurs  on  the  down¬ 
stream  side  of  the  discharge.  With  diffusion  now  in  the  di¬ 
rection  of  the  advcctive  flow,  the  length  of  the  afterglow  is 
extended.  The  downstream  extension  of  the  plasma  zone  is 
aided  by  rarefaction  of  the  gas  downstream  of  the  electrodes 
that  provides  for  more  rapid  charge  particle  diffusion,  and 
dissociation  of  02  that  reduces  the  rate  of  loss  of  electrons  by 
dissociative  attachment. 

Gas  temperature,  Tr  and  densities  of  02,  O.  02(IS),  and 
CM1  A)  are  shown  in  Fig.  2(b)  for  the  base  case.  T v  increases 
by  33  K  above  ambient  by  Joule  and  Frank-Condon  heating. 
The  peak  in  the  gas  temperature  is  shifted  downstream  by  the 
gas  flow.  Heating  occurs  primarily  near  the  electrodes  where 
the  electric  field  is  largest  (see  the  local  peak  in  TK  near  the 
upstream  electrode)  and  secondarily  off  axis  in  the  bulk 


Power  (0.2  W/cm3  (2  dec))  1 3  MHz  40  W  CW 

B  II 

T.  (2.9  eV) 

I  /r-  Cll1  1"i  II 

Iel  (9.8  x  109cm  3  (3  dec)) 

. "''"kggBgEB^  II 


T  ^333  K1  i  q  uu->  ac\  va/  r'vA/ 


min  mm.m: t  km  max 

FIG.  2.  (Color)  Plasma  properties  and  species  densities  for  Ihe  base  case 
conditions  (He/O2=7()/30,  3  Torr,  6  Ipm.  40  W)  when  using  cw  excitation 
at  13.56  MHz.  (a)  Power  deposition,  electron  temperature,  and  electron  den¬ 
sities:  and  (b)  gas  temperature  and  densities  of  0»,  O.  0:(  ’S).  and  0>{  ’A) 
All  values  are  averaged  over  one  rt  period.  The  flow  is  from  the  right.  The 
scales  are  linear  with  zero  minimum  values  with  the  exception  of  TfJS  unless 
the  number  of  decades  (i.e.,  “decs”)  is  indicated  for  log  plots  The  maximum 
value  is  indicated  in  each  figure.  Depletion  in  ground  state  0:  is  due  to  gas 
heating  by  the  localized  power  deposition.  0:('A)  maximizes  downstream 
in  part  due  to  quenching  of  02('~). 

plasma  where  the  ion  density  is  largest.  Thermal  conduction 
and  heat  transfer  at  the  walls  shift  the  peak  in  7\,  to  the  center 
of  the  tube.  The  gas  returns  to  near  ambient  values  at  the 
exit. 

Dissociation  and  excitation  processes  contribute  to 
depletion  of  ground  state  02  in  addition  to  that  due  to  rar¬ 
efaction  by  gas  heating.  Electron-impact  dissociation  pro¬ 
duces  about  6%  of  the  depletion  with  an  additional  5%  com¬ 
ing  from  electron-impact  excitation  to  02('A)  and  O^1!). 
The  density  of  O  atoms  peaks  26  cm  downstream,  decreas¬ 
ing  thereafter  due  to  recombination  on  the  walls,  and  three- 
body  association  reactions  in  the  gas  phase  to  form  02  and 
Ov  The  density  of  02(,!£)  peaks  in  the  plasma  zone  where 
production  is  largest.  As  the  density  of  O  atoms  increases. 
02('X)  is  rapidly  quenched  to  02('A)  by  collisions  with  O 
atoms  in  the  flow  direction  after  the  discharge.  After  the 
02(,S)  is  converted  to  CM1  A),  the  density  of  CM1  A)  only 
moderately  decreases  along  the  tube  due  to  quenching.  Ex¬ 
cited  states  of  02  and  atomic  O  also  extend  5-  1(1  cm  up¬ 
stream  of  the  electrodes,  a  consequence  of  the  plasma  ex¬ 
tending  upstream  and  back  diffusion  against  the  flow. 

The  combined  yields  of  02('A)  and  02('l)  on  axis  for 
cw  excitation  are  shown  in  Fig.  3  for  different  powers  but  for 
otherwise  the  base  case  conditions.  The  yield  is  maximum  at 
the  edge  of  the  plasma  zone  near  the  downstream  electrode. 
The  yield  is  nearly  constant  thereafter,  only  slowly  decreas¬ 
ing  due  to  heavy  particle  quenching.  Yields  increase  linearly 
with  lower  power  deposition  and  begin  to  saturate  at  high 
powers.  This  is  due  to  increased  power  losses  to  the  clcc- 
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FIG.  3.  Combined  yield  of  02('A)  and  (^(‘X)  along  the  axis  of  the  tube 
(He/02=70/30,  3  Torr.  6  Ipm)  at  13.56  MHz  for  different  powers.  The 
flow  is  from  the  right.  Maximum  yields  occur  within  the  range  of  the  elec¬ 
trodes.  Yields  scale  linearly  with  lower  power  deposition  and  begin  to  satu¬ 
rate  at  high  powers. 

trodes,  principally  ion  acceleration,  typical  for  capacitive  dis¬ 
charges.  For  example,  for  40  W  (specific  energy  deposition 
of  0.31  eV/02),  the  yield  is  about  4.5%  at  the  exit  of  the 
tube  or  an  efficiency  of  0.15  02('A)  molecules/cV  of  energy 
deposition.  For  340  W  (specific  energy  deposition  of 
2.95  cV/02),  the  yield  is  19.5%  at  the  end  of  the  plasma 
zone  and  18%  at  the  exit  of  the  tube,  an  efficiency  of  0.04 
CLl'A)  moIcculcs/cV. 

B.  Spiker-sustainer  excitation 

cw  self-sustained  discharges  require  that  the  sources  ol 
ionization  and  charged  particle  loss  be  equal.  For  moderate 
pressure  mixtures  of  He/02.  this  translates  into  a  balance 
between  electron-impact  ionization  and  electron  loss  by  dis¬ 
sociative  attachment  to  02  which  requires  2-3  eV.  Un¬ 
fortunately,  the  value  of  T,  at  which  power  dissipation  into 
excitation  of  02(‘A)  and  O^1^)  is  maximum  is  about 
1-1.5  eV.  This  value  is  largely  independent  of  gas  mixtures 
as  it  is  determined  by  the  values  of  the  excitation  cross  sec¬ 
tions.  As  such,  significant  efforts  have  been  expended  in  us¬ 
ing  gas  additives  to  lower  the  E/N  of  the  discharge  so  that  Tf 
is  1-1.5  eV.  In  most  cases,  these  values  are  inaccessible  to 
self-sustained  electric  discharges.  In  SS  systems  for  eCOIL 
applications,  the  discharge  parameters  are  selected  so  that  a 
significant  fraction  of  the  power  deposition  occurs  at  Te  be¬ 
low  self-sustaining  so  that  the  rates  of  direct  electron-impact 
excitation  of  CM1  A)  and  02('S)  are  maximized. 

Lowering  the  time  averaged  Te  to  match  the  cross  sec¬ 
tions  of  desired  excited  states  is  a  well  known  strategy.  For 
example,  externally  sustained  electric  discharges,  such  as 
electron  beam  sustained  discharges  (EBSDs).  have  long^bccn 
used  to  excite  C02  lasers  in  He/N2/C02  mixtures.  0  The 
self-sustaining  E/N  and  Tt.  of  these  mixtures  typically  ex¬ 
ceed  that  which  optimizes  excitation  into  the  upper  laser 
level,  either  by  direct  electron-impact  vibrational  excitation 
of  CO?  or  through  excitation  transfer  from  N2(u=l).  The 
external  ionization  provided  by  the  electron  beam  enables 
the  E/N  provided  by  the  discharge,  for  a  given  power  depo¬ 
sition.  to  be  below  the  self-sustaining  value  and  be  a  better 
match  to  vibrational  excitation  of  C02  and  N2.  The  disad- 


FIG.  4.  Plasma  properties  when  the  He/02  gas  mixture  at  3  Toit  is  excited 
by  a  single  spiker  pulse,  (a)  Power  deposition,  (b)  electron  density,  and  (c) 
electron  temperature.  The  excess  electron  density  produced  by  the  spiker 
enables  the  electron  temperature  to  decrease  to  a  more  optimum  value  lor 
excitation  of  02( *21 ).  After  the  pulse  the  steady-state  level  is  regained  after 
15-20  /xs. 

vantage  of  this  approach  is  the  added  complication  and  ex¬ 
pense  of  the  electron  beam  hardware.  EBSD  production  of 
02(!A)  was  investigated  by  Ionin  et  al . 13 

The  SS  technique  attempts  to  achieve  the  same  ends  as 
the  EBSD  while  using,  in  principle,  a  single  discharge  appa¬ 
ratus.  Using  the  SS  technique,  the  excess  ionization  provided 
by  the  spiker  enables  the  E/N  and  Te  during  the  suslaincr 
period,  for  a  given  time  averaged  power  deposition,  to  be 
below  self-sustaining  similar  to  the  EBSD.  The  excess  ion¬ 
ization  in  the  SS  technique  is  provided  in  situ  and,  unlike  the 
EBSD.  is  transient.  That  is,  the  discharge  will  eventually 
recover  to  a  steady,  self-sustaining  state  having  a  larger  T,.. 

The  SS  method  we  are  investigating  consists  of  pulsed 
modulated  rf  excitation;  a  scries  of  high  power  rf  pulses  fol¬ 
lowed  by  a  period  of  lower  power  rf  excitation.  As  a  point  of 
departure,  plasma  properties  during  a  single  spiker  pulse  fol- 
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FIG.  5.  Typical  values  of  (a)  voltage,  and  (b)  power  for  SS  excitation  with 
I  MHz  repetition  rate.  The  carrier  frequency  is  27  MHz  The  Vss/Vcw  ratio 
is  5/2.  Power  is  120  W  when  averaged  over  a  SS  period  (I  /us). 
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FIG.  6.  (Color)  Plasma  properties  and  species  densities  for  the  base  case 
conditions  (He/02  =  70/30,  3  Torr,  6  Ipm.  40  W)  when  using  SS  excitation 
at  13.56  MHz  with  a  duly  cycle  of  20%.  (a)  Power  deposition,  electron 
temperature,  and  electron  density;  and  (b)  gas  temperature  and  densities  of 
02.  O.  02( '2).  and  02('A)  All  values  are  averaged  over  one  SS  period. 
Labeling  is  the  same  as  in  Fig.  2.  The  flow  is  from  the  right.  With  SS 
excitation  the  electron  temperature  is  lower  and  maximum  electron  density 
higher  compared  with  cw  excitation. 


lowed  by  a  long  sustainer  period  will  be  described.  Using  the 
base  case  conditions,  the  discharge,  which  is  first  operated  in 
a  steady  state  at  40  W,  is  then  excited  by  a  single  short  high 
power  pulse  of  270  W  of  duration  of  200  ns  and  is  then 
followed  by  a  (sustainer)  period  of  40  W.  The  resulting 
power,  and  ne  and  Te  on  axis  midway  between  the  electrodes 
(approximately  at  the  peak  plasma  density)  are  shown  in  Fig. 
4  as  a  function  of  time.  The  spiker  pulse  avalanches  ne  to 
values  as  high  as  2.1  X  10,u  cm"3  above  the  steady-state 
value  of  ne(,=  1.1  X  10,()  cm"3.  As  recombination  and  attach¬ 
ment  consume  the  excess  ionization.  nr  decreases  towards 
neo  over  a  period  of  15-20  /zs.  The  quasi-steady-state  elec¬ 
tron  temperature  is  7^=2.!  cV.  Tc  increases  during  the 
spiker  to  avalanche  the  gas,  thereby  creating  excess  ioniza¬ 
tion  and  then  decreases  below  Ten  after  the  spiker  pulse  ter¬ 
minates.  As  long  as  nr  is  above  nro  then  is  below  Teo, 
falling  to  as  low  as  1.6  eV.  As  the  electron  density  decays 
towards  neo,  then  Tr  increases  towards  Teo.  These  lower  val¬ 
ues  of  7,,  coupled  with  the  excess  ionization,  have  the  po¬ 
tential  to  increase  the  production  of  Oi^A). 

Typical  voltage  and  power  wave  forms  for  a  discharge 
operating  in  a  SS  mode  for  a  time  averaged  power  of  120  W 
are  shown  in  Fig.  5.  The  conditions  are  otherwise  the  same 
as  the  base  case.  Excess  ionization  is  provided  by  means  of  a 
high  voltage  spiker  pulse  ( Vss)  of  200  ns  duration  repeated 
at  1  MHz  at  the  carrier  frequency  of  27  MHz  yielding  a  duty 
cycle  of  20%.  In  between  the  spiker  pulses,  the  plasma  is 
sustained  by  a  cw  rf  voltage  (Vcw).  The  choices  of  VsS/Vcw 
=5/2  and  the  triangular  spiker  pulse  shape  are  somewhat 
arbitrary  and  have  not  been  optimized  for  maximum  yield  of 


02(iA).  Power  is  averaged  over  the  full  SS  period  of  I  /zs 
and  the  voltages  V'ss  and  Vcw  are  adjusted  to  maintain  the 
total  time  averaged  power  deposition  at  the  specified  value 
Typically  up  to  15 — 20  (15-20 /xs)  of  SS  periods  were 
simulated  in  order  for  the  system  to  come  into  a  pulse- 
periodic  steady  state.  For  the  results  in  Fig.  5  (time  averaged 
power  of  120  W)  the  power  during  the  spiker  reaches  320  W 
with  a  peak  ^ss  =900  V  resulting  in  the  power  during  the 
sustainer  being  approximately  90  W  with  Vcw=360  V. 

Time  averaged  plasma  properties  and  species  densities 
for  the  base  case  conditions  when  using  SS  (duty  cycle  20%. 
Vss/ Pcw=5/2,  1  MHz  repetition  rate,  and  1 3.56  MHz  carrier 
frequency)  arc  shown  in  Fig.  6  for  40  W  average  power. 
Values  have  been  averaged  over  the  entire  1  /zs  SS  period. 
(Compare  these  results  to  Fig.  2  for  cw  excitation.)  With  SS 
excitation,  the  time  averaged  Tt.  is  lower  in  the  volume  where 
the  electron  density  is  maximum  while  the  lime  averaged  nr 
is  larger.  With  SS  excitation,  the  time  averaged  Te=  1.7  eV 
where  n(.  is  maximum  whereas  with  cw  excitation.  7', 
=  2.2  eV.  The  Te  remains  near  its  cw  value  adjacent  to  the 
electrodes  where  sheath  heating  dominates.  The  time  aver¬ 
aged  maximum  value  of/?,  increases  from  9.8  x  109  cm  with 
cw  excitation  to  2  X  10,0cm"3  with  SS  excitation.  The  de¬ 
crease  in  Te  and  increase  in  nt  are  both  favorable  for  more 
efficient  CM1  A)  production.  The  relative  spatial  distributions 
of  species  densities  [O.  02.  02('A).  and  02('l)]  do  not  sig¬ 
nificantly  differ  from  the  cw  case.  The  exception  is  the  spa¬ 
tial  distribution  of  TK  that  has  more  localized  heating  at  the 
electrodes. 
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FIG.  7.  Oxygen  atomic  and  molecular  densities  (a)  and  gas  temperature  (b) 
along  the  axis  for  the  base  case  conditions  (He/O2  =  70/3().  3  Torr,  6  Ipm, 
40  W)  when  using  cw  (dotted  line)  and  SS  (solid  line)  excitations.  The  flow 
is  from  the  right.  The  degree  of  dissociation  decreases  and  02('A)  produc¬ 
tion  increases  with  SS  excitation. 

Time  averaged  species  densities  and  gas  temperature  on 
the  axis  of  the  tube  are  compared  for  cw  and  SS  excitations 
in  Fig.  7.  With  SS  excitation,  the  lower  value  of  Te  results  in 
less  electron-impact  dissociation  of  02,  producing  a  smaller 
density  of  O  atoms,  while  increasing  the  rate  of  excitation 
02(‘A).  The  contribution  of  quenching  of  02(]~)  to  produc¬ 
tion  of  02(,A)  by  collisions  with  O  atoms  is  about  20%.  This 
quenching  is  somewhat  more  rapid  in  the  cw  case  due  to  the 
larger  density  of  O  atoms.  The  production  of  03  is  also  larger 
in  the  cw  case  due  to  the  larger  production  of  O  atoms,  as 
shown  in  Fig.  7(b).  The  peak  gas  temperature  is  400  K  with 
SS  excitation,  whereas  with  cw  excitation  the  peak  TH  is  only 
about  333  K.  This  is  counterintuitive  as  the  smaller  amount 
of  Frank-Condon  heating  due  to  there  being  less  dissociation 
with  SS  excitation  should  produce  a  lower  TK.  The  larger 
time  averaged  value  of  nr  and  lower  value  1\  obtained  with 
SS  excitation  produce  more  vibrationally  excited  02  that  un¬ 
dergoes  V-T  relaxation  and  more  charge  exchange  healing. 
As  a  result,  Ty  increases  with  SS  excitation. 

The  time  averaged  n(J  and  Te  at  the  location  of  the  maxi¬ 
mum  electron  density  are  shown  in  Fig.  8  as  a  function  of 
power  for  cw  and  SS  excitations.  For  powers  between  10  and 
150  W.  Tf  is  approximately  0.5  cV  lower  with  SS  excitation 
compared  to  cw  excitation,  while  the  electron  density  is  ap¬ 
proximately  a  factor  of  2  larger.  Note  that  with  increases  in 
power.  Te  decreases  for  both  cw  and  SS  excitations.  This 
results  from  the  larger  degree  of  dissociation  of  02  and  larger 
density  of  excited  states.  The  larger  dissociation  of  02  results 
in  there  being  less  dissociative  attachment  and  fewer  electron 
losses,  and  so  a  smaller  Te  is  required  to  sustain  the  dis¬ 
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FIG.  8  Plasma  parameters  as  a  function  of  power  for  cw  and  SS  excitations, 
(a)  Electron  temperature  and  (b)  peak  electron  density.  The  electron  tem¬ 
perature  is  measured  at  the  location  of  peak  electron  density.  The  values  are 
averaged  over  a  SS  period  of  I  /its. 

charge.  The  larger  excited  state  densities  enable  a  larger  rate 
of  multistcp  ionization  which  also  allows  for  a  smaller  value 
of  Tc. 

Atomic  oxygen  and  02(*A)  densities  on  axis  for  40  and 
120  W  are  shown  in  Fig.  9.  The  exit  density  of  02(,A)  in¬ 
creases  with  SS  compared  to  cw  excitation  whereas  the  O 
atom  density  decreases.  The  location  of  the  maximum  O 
atom  density  progressively  moves  downstream  with  increas¬ 
ing  power  due  in  large  part  to  the  increasing  gas  temperature, 
rarefaction,  and  extension  of  the  plasma  zone. 

The  combined  yields  of  O^'A)  and  02('S)  along  the 
axis  and  at  the  exit  of  the  tube  for  cw  and  SS  excitations  for 
different  powers  are  shown  in  Fig.  10.  The  efficiency  of  pro¬ 
duction  of  02(,A).  molcculcs/cV  of  energy  deposition,  is 
also  shown.  For  these  discharge  parameters,  yields  arc  typi¬ 
cally  higher  for  SS  excitation  for  low  powers  and  converge  at 
higher  powers  ( 150  W).  For  example,  at  low  power  (20  W) 
the  energy  deposition  required  to  produce  one  02('A)  mol¬ 
ecule  is  7.3  eV  for  cw  and  5.6  eV  for  SS  excitation,  netting  a 
30%  improvement  in  02('A)  yield.  This  improvement  in 
yield  diminishes  as  the  total  power  is  increased  until,  for 
these  conditions,  the  advantage  of  SS  excitation  is  lost  for  a 
power  of  120  W.  The  decrease  in  improvement  using  the  SS 
is  likely  the  natural  decrease  in  Te  that  occurs  when  increas¬ 
ing  power  with  cw  excitation,  as  shown  in  Fig.  8.  As  1\ 
decreases  towards  I  eV,  the  improvement  in  efficiency  of 
exciting  02('A)  obtained  by  lowering  Tr  another  0.5  eV  with 
SS  excitation  is  diminished. 
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FIG.  9.  Atomic  oxyger  and  02('A)  along  the  axis  for  40  and  120  W  (a)  cw 
and  (b)  SS.  The  flow  s  from  the  right.  Dissociation  of  Ch  is  significantly 
suppressed  with  SS  excitation. 


The  choice  of  carrier  frequency  is  important  due  to  the 
dependence  of  ionization  rate  and  Te  on  frequency.  For  ex¬ 
ample,  the  combined  yields  of  02('A)  and  along  the 

axis,  exit  yield,  and  efficiency  for  a  carrier  frequency  of 
27  MHz  using  cw  and  SS  excitations  are  shown  in  Fig.  1 1 
for  different  powcis.  As  with  13.56  MHz,  yields  are  typically 
higher  for  SS  compared  to  cw  excitation  at  low  powers.  The 
efficiency  obtained  with  two  methods  also  converges  at  high 
powers.  At  20  W,  the  energy  deposition  required  to  produce 
an  Oil1  A)  molecule  with  SS  excitation  decreases  to  4.6  cV 
at  27  MHz.  This  enables  a  42%  improvement  in  efficiency 
compared  to  cw  excitation,  while  at  13.56  MHz,  the  im¬ 
provement  is  30%  for  the  same  conditions. 

O^'A)  production  efficiency  increases  with  increasing 
carrier  frequency  with  cw  excitation  as  well.  For  example, 
02( 1  A)  production  efficiency  as  a  function  of  power  with  cw 
excitation  for  13.56  and  27  MHz  is  shown  in  Fig.  12(a).  The 
efficiency  is  larger  at  27  MHz  for  powers  up  to  150  W.  At 
both  frequencies,  efficiency  first  increases  with  power  and 
then  saturates.  The  increase  with  power  results  from  the  de¬ 
crease  in  Tt.  (more  dissociation  of  02  and  less  dissociative 
attachment)  that  enables  better  matching  with  cross  section 
for  excitation  of  CM1  A).  This  advantage  is  offset  by  the  dis¬ 
sociation  of  02  anJ  rarefaction  that  occurs  at  higher  power. 
Although  increasing  frequency  also  increases  the  efficiency 
of  excitation  with  SS  excitation,  the  trend  with  power  is  op¬ 
posite  to  that  with  cw  excitation,  as  shown  in  Fig.  12(b). 
Since  Tt,  is  already  low  [and  near  optimum  for  CM1  A)  pro- 


FIG.  10.  Plasma  parameters  as  a  function  of  power  for  cw  and  SS  excita¬ 
tions  at  13.56  MHz.  (He/0:  =  70/30.  3  Torr,  6  Ipm).  (a)  Combined  yield  of 
02(’A)  and  Ojf'X)  along  the  axis  for  flow  from  the  right,  (b)  combined 
yields  at  the  outlet  of  the  tube,  and  (c)  efficiency  of  0:('A)  production 
Yields  are  typically  higher  for  SS  regime  for  low  power  but  approach  cw 
excitation  at  high  powers. 


duction]  with  SS  excitation,  the  decrease  in  02  density  that 
results  from  electron-impact  dissociation  with  increasing 
power  serves  to  lower  the  fractional  power  expended  in 
02('A)  excitation. 

The  scaling  of  yield  with  carrier  frequency  depends  on 
the  method  of  excitation,  as  shown  in  Fig.  13(a).  When  in¬ 
creasing  frequency,  Tf  generally  decreases  due  to  the  more 
efficient  electron  ionization  and  reduction  in  power  loss  to 
ion  acceleration  that  occurs  at  higher  frequencies,  as  shown 
in  Fig.  13(b).  With  cw  excitation,  the  decrease  of  Tr  with 
increasing  frequency  up  to  40  MHz  brings  Tf  closer  to  the 
optimum  range  of  1-1 .5  eV.  As  a  result,  the  yield  of  02( 1  A) 
monotonically  increases  attaining  a  maximum  at  about 
33  MHz.  With  SS  excitation.  Te  is  lower  than  with  cw  power 
and  near  the  optimum  value  at  27  MHz.  Increasing  fre- 
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FIG.  1 1.  Plasma  parameters  as  a  function  of  power  for  cw  and  SS  excita¬ 
tions  at  27  MHz.  (He/02=70/30,  3  Torr,  6  Ipm).  (a)  Combined  yield  of 
02('A)  and  02(l~)  along  the  axis  for  flow  from  the  right,  (b)  combined 
yields  at  the  outlet  of  the  tube,  and  (c)  efficiency  of  02('A)  production. 

quency  from  13.56  to  27  MHz  decreases  Te  into  the  opti¬ 
mum  range  producing  an  increase  in  yield.  Further  increases 
in  frequency  decrease  Te  to  being  below  optimum  for  02('A) 
production  and  so  yield  decreases. 

IV.  OPTIMIZING  SS  EXCITATION 

Optimizing  O-^A)  yield  using  SS  excitation  depends 
upon  the  plasma  dynamics  that  result  from  the  details  of  the 
pulse  shape,  as  shown  in  Fig.  14(a).  These  characteristics 
include  the  carrier  frequency,  duty  cycle  (fraction  of  the  SS 
cycle  for  the  spiker),  SS  frequency  (time  between  spiker 
pulses),  spiker  pulse  shape,  and  value  of  Vss/Vcw.  For  ex¬ 
ample.  combined  yields  as  function  of  Vss/Vcw  are  shown  in 
Fig.  14(b).  Maximum  yields  are  obtained  for  Vss/Vcw 
=  2.5-3.  Vss/Vcw=l  corresponds  to  cw  excitation.  As 


FIG.  12.  Comparison  of  efficiency  of  02('A)  production  for  13.56  and 
27  MHz  with  (a)  cw  and  (b)  SS  excitations.  Using  higher  earner  frequencies 
is  favorable  for  both  cw  and  SS  excitations. 


FIG.  13.  Plasma  parameters  as  a  function  of  rf  for  cw  and  SS  excitations,  (a) 
Combined  yield  at  the  end  of  the  lube  and  (b)  electron  temperature.  Electron 
temperatures  generally  decrease  with  increasing  carrier  frequency  and  may 
become  too  low  for  optimum  0>('A)  production 
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FIG  14.  Optimization  of  CM1  A)  production  with  SS  excitation,  (a)  Voltage 
wave  form  for  a  duty  cycle  of  20%.  Combined  yield  as  a  function  of  (b)  the 
ratio  of  spiker  to  sustaircr  voltage,  (c)  duty  cycle  for  SS  excitation  at  10  and 
27  MHz  carrier  frequencies,  and  (d)  sustainer  pulse  length.  Increasing  the 
duty  cycle  generally  decreases  yield  as  the  discharge  appears  to  operate  in  a 
cw  mode. 

Vss/V cw  increases  from  one,  the  amount  of  excess  ionization 
increases,  thereby  enabling  a  lower  value  of  Tc  during  the 
sustainer.  When  keeping  the  total  power  constant,  as  Vss/Vc„ 
increases  a  larger  raction  of  the  power  is  expended  during 
the  avalanche  when  Te  is  large  and  the  efficiency  of  exciting 
CM1  A)  is  low.  Therefore,  very  large  values  of  Vss/Vcw  are 


not  optimum.  Vss/Vcw  should  be  large  enough  to  produce 
sufficient  excess  ionization  to  last  through  the  sustainer  pe¬ 
riod.  but  small  enough  that  the  majority  of  the  power  is  dis¬ 
sipated  during  the  sustainer  when  Tr  is  low. 

Combined  yields  as  a  function  of  spiker  pulse  length 
(duty  cycle)  for  10  and  27  MHz  carrier  frequencies  are 
shown  in  Fig.  14(c)  for  a  total  SS  period  of  I  /zs.  Maximum 
yields  are  obtained  for  duty  cycles  of  l0%-20%  depending 
on  carrier  frequency.  The  ideal  situation  would  be  a  delta- 
function  spiker  that  produces  a  large  electron  temperature, 
efficiently  creating  large  amounts  of  excess  ionization  and 
depositing  little  energy.  The  reality  is  that  a  finite  time,  even 
for  large  electron  temperatures,  is  required  to  produce  the 
excess  ionization.  For  these  conditions,  that  finite  lime  is 
100-200  ns.  There  are  practical  limitations  for  the  length  of 
the  spiker  since  at  least  a  rf  cycle  or  two  arc  required  during 
the  spiker.  At  27  MHz  and  20%  duty  cycle,  there  arc  five  rf 
periods  during  the  spiker,  whereas  at  10  MHz  and  10% 
spiker  duty  cycle,  there  is  only  one  rf  period  during  the 
spiker.  At  the  other  extreme,  as  the  duty  cycle  approaches 
unity  SS  excitation  appears  to  be  cw  excitation  as  quasi¬ 
steady-state  conditions  are  obtained  during  the  spiker.  Yield 
therefore  decreases  and  approaches  the  cw  value. 

Increasing  the  sustainer  pulse  length  from  small  values 
(<1  /zs)  increases  O^C'A)  production.  For  these  conditions, 
02(,A)  yield  optimizes  for  a  sustainer  pulse  length  of  1 .8  /zs, 
as  shown  in  Fig.  14(d).  The  longer  sustainer  pulse  enables  a 
longer  period  of  operating  at  low  Tt,.  While  keeping  the  total 
power  constant,  extending  the  sustainer  reduces  the  amount 
of  power  dissipated  during  the  spiker  and  so  reduces  the 
amount  of  excess  ionization.  With  sufficient  ionization  dur¬ 
ing  the  spiker,  the  optimum  sustainer  length  should  approach 
the  10-15  /zs  typically  required  for  the  plasma  to  recover  to 
a  steady  state. 

Based  on  these  and  other  scaling  studies  we  find  that 
02(,A)  production  optimizes  using  SS  techniques  for  the  fol¬ 
lowing  conditions. 

•  Operate  with  a  high  carrier  frequency  to  lower  Tf  and 
increase  proportion  of  power  dissipated  by  electrons. 

•  The  spiker  pulse  length  should  be  as  short  as  possible 
while  still  providing  sufficient  excess  ionization  to  last 
through  the  entire  sustainer  period. 

•  The  spiker  voltage  should  not  elevate  the  electron  tem¬ 
perature  above  that  required  for  rapid  ionization. 

•  The  sustainer  pulse  length  should  be  long  enough  that 
the  majority  of  power  is  dissipated  during  the  sustainer 
but  not  so  long  that  the  plasma  fully  recovers  to  its  cw 
conditions. 


V.  CONCLUDING  REMARKS 

The  consequences  of  spiker-sustain  excitation  on  yield 
of  02('A)  in  flowing  He/02  plasmas  were  investigated  with 
2D  plasma  hydrodynamics  model.  SS  techniques  generally 
do  lower  the  time  averaged  Tt  compared  to  cw  excitation  and 
so  improve  the  efficiency  of  excitation  and  yield  for  02(‘A). 
This  advantage  to  SS  excitation  is  diminished  at  higher  pow¬ 
ers  where  Tt,  is  naturally  lower  for  cw  excitation.  Lower  duty 
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cycles  (shorter  spiker  pulses)  are  generally  more  advanta¬ 
geous  for  SS  excitation  as  the  limit  of  a  delta-function  spiker 
is  approached.  Increasing  the  frequency  of  excitation  is  also 
advantageous  due  to  the  lowering  of  Te  and  increase  in  ne. 
Again,  the  advantage  of  SS  excitation  is  diminished  as  fre¬ 
quency  increases  due  to  the  natural  decrease  in  Te  with  cw 
excitation.  The  length  of  the  sustainer  pulse  should  be  short 
enough  so  that  Te  is  largely  below  the  self-sustaining  value 
but  long  enough  so  that  Te  and  ne  begin  recovering  towards 
their  steady-state  values.  This  recovery  indicates  a  good  uti¬ 
lization  of  the  excess  ionization  produced  by  the  spiker. 
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Abstract 

The  1.315 /tm  [I(2Pi/2)  -*  I(2P3/2)]  transition  of  atomic  iodine  in  the 
chemical  oxygen-iodine  laser  (COIL)  is  pumped  by  sequential  reactions  of 
I2  and  I  with  02(1  A).  In  electrically  pumped  systems  (eCOILs),  electron 
impact  excitation  of  02  produces  the  02(*  A)  and  also  produces  O  atoms 
through  dissociative  excitation.  The  O  atoms,  through  reactions  with  I2, 

I(2P,/2)  and  I(2P3/2),  lead  to  dissociation  of  I2,  quenching  of  the  upper  laser 
level  and  removal  of  the  lower  laser  level.  While  dissociating  I2  is 
potentially  beneficial,  quenching  of  the  upper  laser  level  is  detrimental  and 
so  management  of  the  O  atom  density  is  necessary  to  maximize  laser  gain. 

In  this  regard,  NO  and  N02  additives  have  been  used  to  manage  the  O  atom 
density  by  cyclically  reacting  with  O  and  I.  In  this  paper,  results  from  a 
computational  investigation  of  eCOIL  systems  using  plug  flow  and 
two-dimensional  models  are  discussed  where  NO  and  N02  additives  are 
used.  The  system  is  a  flowing  plasma  sustained  in  He/02/NO  mixtures  with 
downstream  injection  of  N02  followed  by  injection  of  I2.  We  found  that 
addition  of  NO  and  N02  is  effective  in  managing  the  density  of  O  atoms  and 
maximizing  gain  by  minimizing  quenching  of  the  upper  laser  level.  We 
found  that  by  optimizing  the  additives,  laser  gain  can  be  maximized  even 
though  02('  A)  densities  may  be  lower  due  to  the  management  of  quenching 
and  dissociation  reactions. 


1.  Introduction 

Chemical  oxygen-iodine  lasers  (COIL)  operating  on  the 
1.315  fxm  [I(2Pi/2)  -*  I(2P3/2)]  transition  of  atomic  iodine  are 
being  investigated  due  to  their  high  efficiency  and  potential 
for  multi-kilowatt  CW  operation  [1-3].  A  series  of  collisional 
transfer  reac  ions  between  O^1  A),  I2  and  ground  state  l(2P3/2) 
(to  be  referred  to  as  I)  result  in  excitation  of  the  upper 
laser  level  I(2Pi/2)  (to  be  referred  to  as  P).  Typically 
02(!A)  is  generated  upstream  of  the  laser  cavity.  I2  is 
injected  immediately  prior  to  the  cavity  upon  which  the  flow 
is  supersonically  expanded  to  lower  the  gas  temperature  as 

3  Author  to  whom  any  correspondence  should  be  addressed. 


required  to  maximize  the  gain.  In  conventional  COILs  liquid- 
phase  chemistries  (reactions  between  basic  H202  and  Cl2) 
produce  the  02(,A)  with  high  yields  [4],  although  the  use 
of  liquid  peroxides  and  Cl2  gas  to  produce  the  02(‘  A)  creates 
challenges  for  storage  and  transport.  Recently,  efforts  have 
focused  on  generating  the  02(!A)  in  electrical  discharges 
(eCOILs)  due  to  the  increased  robustness  and  safety  of  the 
all  gas-phase  system  [1-10].  Laser  gain  and  oscillation  have 
been  demonstrated  in  eCOILs  by  Hicks  et  al  [  I  ]  and  Verdeyen 
et  al  [3]. 

A  challenge  in  eCOILs  is  to  produce  sufficiently  high 
yields  of  02(!A),  and  hence  laser  gain,  to  enable  the  laser 
transition  to  be  saturated  and  produce  high  power  [5].  In 
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this  regard,  recent  research  has  focused  on  tailoring  the 
discharge  parameters  [6,  7,  10]  and  using  additives  such 
as  CO,  H2,  D2  [8]  and  NO  [9]  to  improve  the  excitation 
efficiency  of  02(‘  A).  NO  and  NO2  are  also  used  as  additives 
to  control  Ihe  post-discharge  chemistry  [9].  In  fact,  all 
demonstrations  of  laser  gain  and  oscillation  to  date  have 
used  flowing  plasmas  in  He/02  mixtures  with  NO  as  an 
additive. 

eCOILs  differ  from  conventional  COILs  in  that  atomic 
oxygen  is  also  produced  in  the  electric  discharge  by  electron 
impact  dissociation  of  O2.  The  O  atoms  flow  downstream 
where  3-bo<iy  reactions  produce  O3  and  the  remaining  O 
atoms  may  react  with  the  injected  I2.  Beneficial  reactions  of 
O  with  I2  p  oduce  I  atoms,  thereby  eliminating  the  expense 
of  02(' A)  molecules  for  initiating  the  dissociating  reactions. 
Detrimental  reactions  involving  O  atoms  include  quenching 
of  I*  by  O  which  reduces  gain.  As  such,  management  of 
the  O  atom  density  is  important  to  optimize  these  opposing 
effects. 

Injection  of  NO  or  NO2  through  or  downstream  of  the 
discharge  in  eCOIL  systems  has  two  goals;  improving  the 
efficiency  of  direct  production  of  02(1A)  in  the  discharge 
by  electron  impact  and  management  of  the  O  atom  density 
downstream  of  the  plasma.  Including  NO  in  the  gas  stream 
flowing  through  the  discharge  has,  in  part,  the  goal  of 
improving  production  of  O^A).  NO,  having  a  lower 
ionization  potential  (9.26  eV)  than  either  02  or  He,  is  likely  to 
provide  mors  rapid  ionization  with  the  possibility  of  lowering 
the  operating  E/N  (electric  field/gas  number  density)  and 
electron  temperature,  Tc.  Lowering  Te  from  the  values 
typical  of  self-sustained  He/02  mixtures  is  advantageous  in 
more  efficiently  producing  02(*  A)  by  direct  electron  impact 
[6,8,  II].  Discharges  in  NO  have  also  been  known  to 
produce  02(‘A)  in  relatively  large  amounts  even  though 
02  may  not  be  present  as  a  feedstock  gas  [12].  However, 
when  flowing  NO  through  the  discharge,  some  of  the 
power  that  would  otherwise  be  available  to  excite  02  is 
dissipated  by  excitation  and  ionization  of  NO,  an  unwanted 
consequence. 

In  addition  to  possibly  improving  the  production  of 
02(‘ A),  NO  and  N02  are  potentially  effective  in  managing 
the  inventor/  of  O  atoms  through  direct  and  cyclic  reactions 
which  have  ihe  effect  of  converting  O  atoms  back  into  02.  In 
the  context  of  eCOIL  systems  where  I2  is  injected,  O  atoms 
are  beneficial  by  dissociating  I2  (and  product  species  10)  to 
form  I  atom:;  which  are  then  pumped  to  I*  by  collisions  with 
02(*  A).  Th<  O  atoms  are  detrimental  by  quenching  P.  Totally 
eliminating  ()  atoms  is  therefore  not  necessarily  beneficial.  NO 
and  N02  also  have  secondary  effects  in  that  they  react  with  I 
atoms  forming  intermediary  species  such  as  INO,  and  IN02 
which  further  react  with  I  to  reform  I2  [  1 3]. 

In  this  paper,  we  report  on  results  from  a  computational 
investigation  of  the  consequences  of  NO  and  N02  additives  on 
flowing  He/02  plasmas  and  their  afterglows  with  I2  injection 
in  the  context  of  eCOIL  systems.  These  investigations  were 
conducted  using  plug  flow  and  2-dimensional  (2D)  plasma- 
hydrodynamics  models.  Although  for  most  conditions  the 
addition  of  NO  to  the  inlet  gas  stream  reduced  the  density  of 
02  ( 1  A),  it  ultimately  increased  the  densities  of  I*  downstream 
of  the  discharge  through  management  of  the  O  atom  density. 


The  inlet  NO  mole  fraction  also  typically  increased  the  extent 
of  the  region  over  which  positive  laser  gain  could  be  achieved. 
This  can  be  particularly  useful  in  high  speed  flows  where 
mixing  lengths  are  longer.  Injection  of  N02  in  the  post- 
discharge  flow  can  help  in  rapidly  scavenging  ()  atoms  in  two 
body  reactions  (as  compared  with  scavenging  by  NO  which 
proceeds  by  a  3-body  mechanism).  In  general,  addition  of 
N02  in  the  post-discharge  region  improves  laser  gain  [9],  The 
downside  to  N02  injection  is  a  rise  in  gas  temperature  due  to 
the  exothermicity  of  the  reactions  between  N02  and  O. 

The  models  used  in  this  investigation  are  described  in 
section  2  followed  by  a  discussion  of  the  reaction  mechanism 
in  sections  3.  In  section  4  the  consequences  of  NO  flowing 
through  the  discharge  and  its  influence  on  downstream  kinetics 
when  I2  is  injected  are  described.  In  section  5,  the  post¬ 
discharge  kinetics  with  N02  addition  upstream  of  I2  injection 
are  discussed.  Our  concluding  remarks  are  in  section  6. 

2.  Description  of  models 

This  investigation  was  conducted  using  a  plug  flow  model, 
GlobalKIN ,  and  a  2D  plasma  hydrodynamics  model, 
nonPDPSIM.  GlobalKIN  has  been  previously  described  in 
[11,14]  and  so  will  only  be  briefly  discussed  here.  GlobalKIN 
consists  of  a  volume  averaged  plasma  chemistry  module 
and  an  electron  energy  transport  module.  The  plasma 
chemistry  module  provides  the  time  rate  of  change  of  species 
densities  based  on  gas-phase  chemistry  and  surface  reactions. 
Electron  temperature,  Te,  and  average  gas  temperature,  7g, 
are  also  solved  for  by  integrating  their  respective  conservation 
equations.  The  electron  energy  transport  module  consists 
of  a  solution  of  Boltzmann’s  equation  for  the  electron 
energy  distribution  (EED)  which  provides  electron  impact 
rate  coefficients  based  on  the  EEDs  and  fundamental  cross- 
sections. 

For  plasmas  flowing  through  cylinders  having  large  aspect 
ratios,  transport  to  the  radial  surfaces  is  taken  into  account  by 
using  a  diffusion  length.  By  simultaneously  calculating  the 
axial  speed  of  the  flow  based  on  constant  pressure,  change 
in  enthalpy,  species  densities,  conservation  of  mass  and  gas 
temperature,  the  integration  in  time  is  mapped  to  axial  position. 
The  resulting  rate  equations  are  integrated  in  time  using  a  stiff 
ordinary  differential  equation  solver.  To  address  the  use  of 
additives,  GlobalKIN  was  modified  to  enable  the  downstream 
injection  of  gases  into  the  flow.  Injection  nozzles  were  treated 
as  point  sources  of  mass,  axial  momentum  and  enthalpy.  Power 
deposition  as  a  function  of  axial  position  must  be  specified 
in  GlobalKIN.  This  distribution  was  estimated  by  averaging 
the  power  deposition  obtained  from  the  nonPDPSIM  over  the 
cross-section  of  the  tube. 

nonPDPSIM  has  been  discussed  in  detail  in  [  1 5, 1 6],  and  so 
will  only  be  briefly  described  here  along  with  pertinent  updates 
to  the  model.  Continuity  equations  for  charged  gas-phase 
species,  surface  charges  and  Poisson’s  equation  for  the  electric 
potential  are  simultaneously  implicitly  integrated  in  time. 
Updates  of  these  quantities  are  followed,  in  a  time  splicing 
manner,  with  updates  of  Tc  and  neutral  species  densities  using 
a  modified  form  of  the  compressible  Navier-Stokes  equations 
for  continuity,  momentum  and  energy  (gas  temperature)  which 
accounts  for  interactions  with  the  plasma.  A  circuit  model  was 
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used  to  interface  the  plasma  with  metal  surfaces  connected  to  a 
voltage  souice  and  circuit  elements  (e.g.  ballast  resistor).  The 
power  supply  voltage  was  adjusted  to  deposit  the  specified 
power  in  the  plasma. 

To  address  the  longer  time  scales  for  the  injection  and 
mixing  of  additives  downstream  the  following  method  was 
used.  It  was  determined  and  validated  that  the  injected  gases 
were  sufficiently  far  downstream  from  the  plasma  that  their 
presence  did  not  affect  the  plasma  properties.  As  such, 
plasma  quartities  (for  example,  electron  density  (nc),  Tc,  rates 
of  electron  impact  reactions,  and  ionization  sources)  were 
calculated  until  a  quasi-steady  state  was  reached  after  which 
these  quantities  were  ‘frozen’  in  the  upstream  plasma  zone. 
Following  tie  establishment  of  the  plasma  properties  only 
the  neutral  fluid  equations  were  integrated  (accounting  for 
injection  of  gases  downstream)  which  did  not  affect  quantities 
upstream  in  the  plasma.  This  enabled  a  longer  time-step 
downstream  where  the  injection  kinetics  occurred. 

3.  Reaction  mechanism 

The  species  included  in  the  reaction  mechanism  are  e.  He, 
He(2S),  He*,  Oj,  02(v),  02(‘  A),  02(‘E),  OJ,  OJ,  O,  0('D), 
O('S),  0~,  0\  03,  OJ,  NO,  NO*.  N02,  NO",  NOj,  N, 
N2,  I2,  IJ,  l(2P3/2),  I(2P|/2),  IO,  INO  and  IN02.  He(2P) 
was  included  for  purposes  of  electron  energy  loss  in  solving 
Boltzmann’s  equation  but  was  lumped  with  He(2S)  in  the 
chemical  kinetics.  C>2(v)  represents  the  sum  of  the  first 
four  vibrational  levels  of  O2.  The  vibrational  levels  of 
N2  were  not  included  as  individual  species  in  the  reaction 
mechanism  as  N2  densities  were  small  and  predominantly 
produced  downstream  of  the  discharge  where  electron  impact 
processes  are  negligible.  Excitation  of  the  first  five  vibrational 
levels  of  NO  (As  =  0.23,  0.46,  0.69,  0.91  and  1.13)  and 
the  electronic  states  NO(a2E)  (As  =  5.48  eV),  NO(c2FI) 
(As  =  6.5  eV)  and  NO (fc4Z)  (As  =  7.58 eV)  were  included 
for  purposes  of  energy  loss  collisions  with  ground  state  NO. 
The  densities  of  these  states  were  not  explicitly  tracked  in  the 
model.  The  negative  ions  NO"  and  NOj  were  included  only  in 
GlobalKIN  after  confirming  that  their  exclusion  does  not  make 
a  significant  change  to  the  reaction  kinetics  in  nonPDPSIM. 
Their  exclusion  from  nonPDPSIM  was  for  the  purpose  of 
speeding  the  calculation. 

The  reaction  mechanism  builds  upon  that  previously 
developed  for  He/02  discharges  and  which  is  discussed  in  [  1 4]. 
The  processes  that  were  added  to  that  reaction  mechanism  to 
account  for  the  injection  of  NO,  NO2  and  I2  are  listed  in  table  I . 

The  majority  of  02(‘ A)  is  produced  by  electron  impact 
excitation  o'  ground  state  O2  in  the  plasma  zone,  and 
secondarily  by  excitation  of  02('£)  followed  by  quenching 
to  O2 ( 1  A)  by  collisions  with  O  atoms. 


e  +  O2  — ►  02(^  A)  +  e, 

a) 

e  +  02  -*  02('  £)  +  e, 

(2) 

o-*-o2(,5:)  02(,A)  +  0, 

(3) 

e  +  O2  -►  O  +  O  + e, 

(4) 

e  +  02  ->  0(,D)  +  0  +  e. 

(5) 

Electron  impact  dissociation  of  O2  (equations  (4)  and  (5))  is  a 
significant  scurce  of  O  atoms  and,  in  part,  the  motivation  for 


injection  of  NO  and  NO2  to  manage  the  resulting  flow  of  O 
atoms.  For  self-sustaining  discharges,  roughly  one  O  atom  is 
generated  for  every  02(‘  A)  produced  in  the  plasma  zone.  By 
lowering  Tc  to  1-1 .5  eV  from  self-sustaining  values  of  2-3  eV 
by,  for  example,  using  the  spiker-sustainer  methods,  the  rate 
of  excitation  of  02(*  A)  can  be  increased  relative  to  electron 
impact  dissociation  [II]. 

The  negative  ion  chemistry  of  He/02  mixtures  is  initiated 
by  dissociative  attachment  and  3-body  reactions, 

e  +  02  — ►  CT  +  O,  (6) 

e  +  02  +  M  -♦  OJ  +  M,  (*  =  3.6  x  10-31  T~0 5 cm6  s_l), 

(7) 

where  M  is  a  third  body  and  k  is  the  rate  coefficient.  Oj 
subsequently  charge  exchanges  with  O  atoms  to  form  O-  or 
associatively  detaches  to  form  O3.  A  similar  process  occurs 
with  O”  to  reform  O2, 

0J  +  0-*02  +  0",  (*  =  1.5  X  lO-^cmV),  (8) 

OJ  +  O  ->•  03  +  e,  (*  =  1.5  X  10-'° cm’s-1),  (9) 

0_+0-*02+e,  (k  =  2  x  10-locm3  s_l).  (10) 

Through  dissociative  detachment  reactions,  both  O"  and  Oj 
react  with  O2  ( 1 A )  to  form  O3 ,  and  to  quench  O2  ( 1 A ) , 

O-  +  O2 ( 1  A)  -*03  +  e,  {k  =  3  x  10-,° cm3  s”1), 

(ID 

Oj  +  O2 ( 1 A )  — ►  02  +  O2  +  e,  (k  =  2  x  10“,0cm3  s_l). 

(12) 

An  important  motivation  for  flowing  NO  through  the  plasma 
zone  is  that  it  has  a  smaller  ionization  potential  than  O2  and 
He  and  so  should  produce  more  ionization  at  a  lower  7*c, 

e  +  NO  — ►  NO*  +  e  +  e,  (At  =  9.26eV).  (13) 

This  ionization  does,  however,  come  at  the  cost  of  channelling 
discharge  power  into  NO  (e.g.  vibrational  excitation  and 
dissociation),  power  that  might  otherwise  be  deposited  into 
O2.  This  power  therefore  does  not  directly  lead  to  generation 
of  02(*  A).  Since  NO*  has  the  smallest  ionization  potential 
among  all  the  atoms  and  molecules  in  this  mechanism, 
significant  amounts  of  NO*  are  formed  by  charge-exchange 
reactions  with  He*,  O*  and  OJ.  Dissociative  recombination  of 
NO*  with  electrons  can  produce  N,  O  and  O(’D),  as  well  as 
contribute  to  gas  heating  through  Frank-Condon  processes, 

e  +  NO*  -*•  N  +  O,  (k  =  3.4  x  l(T77;~05cm3  s_l, 

Al/  =  —2.77  eV),  (14) 

-+  N  +  0('D)  (k  =  6.0  x  10~8r-°5cmV, 

AH  =  — 0.8eV),  (15) 

where  A//  is  the  change  in  enthalpy  of  the  recombination 
event. 

NO  also  intercepts  the  negative  ion  chemistry  through 
charge-exchange  reactions  and  associative  attachment, 

CT  +  NO  -*•  N02  +  e,  (*  =  2.8  x  10-|0cm3s_l),  (16) 
CT  +  N02  -+  NOJ  +  O,  (*  =  1.2  x  lO-’cmV).  (17) 

These  reactions  are  potentially  beneficial  in  two  ways.  They 
increase  the  electron  density  which  can  lower  Tc  into  a  more 
favourable  range  for  exciting  O2OA)  and  they  reduce  the 
density  of  O-  which  would  otherwise  quench  02(‘  A)  to  form 
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Table  1.  Processes  added  to  He/02  reaction  mechanism  with  NO,  N02  and  I2. 


Species 


He 

NO 

e 

He* 

no2 

He* 

O 

N 

o2 

O('D) 

N2 

o2‘ 

O('S) 

h 

O* 

o2 

I 

O- 

02(v) 

I* 

O3' 

o2('A) 

10 

NO* 

c^CD 

INO 

NO- 

o3 

IN02 

no2 

Reaction 

Rate  coefficient* 

Reference 

e  +  NO  -►  N  +  O- 

b 

[IS] 

e  +  NO  ->  N  +  O  +  e 

b 

[19] 

e  +  NO  — ►  NO*  +  e  +  e 

b 

[20] 

e  +  NO  -*N  +  0*  +  e  +  e 

b 

[20] 

e  +  NO*  -►  N  +  O 

3.4  x  10"’  T~°* 

[21] 

e  +  NO*  -►  N  +  Of'D) 

6.0  x  10-’  T-0i 

[21] 

e  +  NO  +  He  ->  NO"  +  He 

1.0  x  10-”  re-°5 

[22] 

e  +  NO  +  NO  -►  NO"  +  NO 

6.5  x  10-JI  7;-05 

[22] 

He*  +  NO  -►  He  +  NO*  +  e 

2.5  x  10-'° 

c 

He*  ♦  NO  -►  He  +  NO* 

1.6  x  10-’ 

[23] 

He*  -1-  NO  -►  He  +  N  +  0* 

4.2  x  10*'° 

[23] 

O*  +  NO  O  +  NO* 

1.7  x  10- 12 

[23] 

O*  +  NO  ->  02  +  NO* 

4.5  x  10- 

[23] 

OJ  +  N  -►  NO*  +  O 

1.2  x  10“ 10 

[23] 

NO".  NOJ  +  M*  -*•  NO.  N02  +  Md 

2.0  X  I0-7  (r,/300)-' 

c 

NO*  +  M"  -►  NO  Me 

2.0  x  10-7  (T,/300)-' 

c 

O-  +  NO  -►  N02  +  e 

2.8  x  ro'°  (r,/300)°s 

[23] 

O"  +  N  -►  NO  +  e 

2.2  X  IO-|0(r,/300)0  5 

[23] 

O"  +  N02  -►  NO"  +  e 

1.2  x  10-9(T,/300)0> 

[23] 

02  +  N  — ►  N02  +  e 

4.0  X  10  ,o(rt/30O)°' 

[23] 

02  +  NOj  ->  N02"  +  02 

7.0  X  10" 10  (Tt/300)o  > 

[23] 

O'  +  NO  -►  NO"  +  02 

2.2  x  10- 12  (r,/300)°  5 

[23] 

NO"  +  NO  -►  NO  +  NO  +  e 

5.0  x  10- 12  (7,/300)°5 

[23] 

NO"  +  He  — ►  NO  +  He  +  e 

2.4  x  I0-'3  (r,/300)°5 

[23] 

NO"  +  N02  -+  NO  +  NO- 

7.4  x  IO"10  (7",/300)°5 

[23] 

NO"  +  02  -*  NO  +  02 

5.0  x  10-|0(r,/300)°5 

[23] 

N02  +  03  -►  O3  +  NOj 

9.4  x  10-'2  (r,/300)°5 

1.4  x  10-|2e-'3,o'r. 

[23] 

NO  +  03  -*>  N02  +  02 

[24] 

N0  +  02('A)  -►  NO+O2 

3.5  x  10- 17 

[23] 

NO  +  O('D)  — ►  N  +  02 

5.0  x  10-15 

[26] 

NO  +  O('D)  -*  NO  +  O 

4.0  x  10  " 

[27] 

NO  +  O  +  M  -►  NOj  +  Mf 

1.0  x  10-31  (T./300)~16 

4.2  x  10'12  e273/r» 

[24] 

N02  +  O  ->  NO  +  02 

[24] 

N02  +  0('D)-*  N02  +  0 

3.2  x  10- 10 

[28] 

NOi  +  Of'D)-*  N0  +  02 

3.0  x  10-‘° 

[28] 

N  +  Oi  -►  NO  +  O 

2.4  x  10- 11  e-S320'r. 

[29] 

N  +  N02  -►  NO  +  NO 

6.1  x  10- 12 

[30] 

N  +  N02  -►  N2  +  0  +  O 

2.4  x  10- 12 

[30] 

N  +  NOj  -►  N2  +  Oj 

1.8  x  10~12 

[30] 

N  +  O  +  M  — ►  NO  +  Mf 

5.5  x  IO-33  el55/T‘ 

[31] 

02(*  A)  +  02  +  O  — ►  02  +  02  +  0 

1.0  x  10- 32 

[10] 

OzCEJ  +  Iz  ->  O2  +  I  +  I 

2.8  x  10"" 

[91 

02('L)  +  l2  -►02(,A)  +  I2 

2.3  x  10'" 

[9] 

02(i£)  +  I2  -►  02  +  I2 

6.0  x  10-'2 

[9] 

02  ( 1 A )  +  I2  -►  02  +  15 

7.0  x  I0-'5 

[9] 

02(' A)  +  I2  ►  02  +  I2 

5.0  x  IO'16 

[9] 

O  +  I2  -►  IO  +  I 

1.4  x  10"10 

[9] 

He  +  I2  — ►  He  +  I2 

9.8  x  10-'2 

[9] 

02  +  I2  — ►  02  +  I2 

4.9  x  10-'2 

[9] 

02(*A)  +  I2  — ►  o2  +  1+  I 

3.0  x  10- 10 

[9] 

02(,A)  +  I-^02  +  I* 

7.7  x  10-"  (r,/300)-'° 

[9] 

o2(,a)  +  i-^o2  +  i 

1.0  x  I015 

[9] 

O3  + 1  -*>  o2  + 10 

2.0  x  10-"  e-,907r* 

[9] 

He  +  r  -*  He  +  I 

5.0  x  10-’* 

[9] 

o2('A)  +  r  -►  o2('5:)  +  i 

8.4  x  10“ 15  (7,/300)3  8  e,<M/r> 

[9] 
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Table  1.  Continued. 


Reaction 

Rate  coefficient* 

Reference 

02('A)  +  l*  -►  o2  +  1 

1.10  X  IO*13 

V) 

o2  +  p  — ►  o;  + 1 

1.0  x  10“ 10  (7g/300)_l  0  e~4a3'Tt 

19] 

O  +  P  ->0  +  1 

8.0  x  IO'12 

[9] 

NO  +  P  NO  +  I 

1.2  x  IO'13 

[9] 

no2  +  P  -►  no2  +  I 

8.5  x  IO"14 

[9] 

i2  + 1*  -►  i;+ 1 

1.40  x  |o~,J  el600/7» 

19] 

10  +  10  -►  02  +  I  +  I 

8.2  x  IO" " 

19] 

NO  +  io  -►  no2  +  I 

4.3  x  IO'12  e-m'T> 

[32] 

0  +  10  -+  Oj  +  I 

1.4  x  IO'10 

[9] 

0+I0-^02(,A)  +  I 

1.5  x  I0-" 

[9] 

I  +  INO  — ►  I2  +  NO 

1.6  x  IO-‘° 

[13] 

I  +  IN02  ->  I2  +  NOj 

8.3  x  IO'" 

[13] 

INO  +  INO  -*  I2  +  NO  +  NO 

8.4  x  IO-'1  e-2m/Tf 

[13] 

IN02  +  IN02  -►  I2  +  N02  +  no2 

2.9  x  io-"  e-2#00/ri 

[13] 

I  + 1  +  I2  -►  I2  +  I2 

3.6  x  IO"30 

[9] 

I  +  I  +  He  -►  I2  +  He 

3.8  x  IO'33 

[9] 

I  +  I  +  o2  — ►  i2  +  o2 

3.3  x  IO"32 

[9] 

I  +  p  + 12  ->  I  + 1  + 12 

3.6  x  IO' 30 

[9] 

i  +  i  +  o2  -►  I2  +  02('A) 

3.7  x  IO"33 

[9] 

I  +  NO  +  He  — ►  INO  -»■  He 

6.0  x  10- 33  (T,/300)-' 0 

[13] 

I  +  NO  +  02  — ►  INO  +  02 

1.6  x  IO" 32 

[13] 

I  +  N02  +  He  -►  IN02  +  He 

1.5  x  IO’31  (Tg/300)-1 0 

[13] 

I  +  no2  +  Oj  -►  ino2  +  02 

2.6  x  I0-” 

[13] 

P  -►  1 

10 

[33] 

a  Rate  coefficient  in  cm3  s_l  for  2-body  reactions,  and  cm6  s  'for  3-body  reactions 
and  s~l  for  radiation  reactions. 

h  Rate  coefficients  calculated  using  cross-section  data  from  the  indicated  reference. 
c  Estimated. 

d  Where  M  is  one  of  the  cations  O*,  Oj,  He*  or  NO*. 
c  Where  M  is  one  of  the  anions  O- ,  Oj  or  O3 . 

f  Where  M  is  one  of  the  major  neutral  species  He,  02,  02(v),  02(‘  A),  O  or  O3. 


O3.  The  negative  ions  NO-  and  NOj  were  not  included 
in  the  2D  model  (to  increase  computational  speed),  whereas 
they  were  included  in  the  plug  flow  model.  Computational 
experiments  were  conducted  using  the  plug  flow  model  to 
quantify  the  effects  of  the  NO”  and  NOj  and  are  discussed 
below. 

One  of  1  he  primary  motivations  of  injecting  NO  and  N02 
is  in  their  potential  for  managing  of  the  O  atom  inventory. 
Much  of  this  chemistry  is  cyclic.  NO  reacts  with  O  to  form 
NO2  which  then  further  reacts  with  O  to  regenerate  02, 

NO  +  O  +  M  -*•  N02  +  M,  (i=I.Ox  l(T3l(7i/300r16, 

AW  =  -3.18eV),  (18) 

N02  +  0->  NO  +  02.  (k  =  4.2  x  10~12 exp(273/rg), 

AW  =  -2.0eV).  (19) 

This  is  an  important  reaction  chain  that  can  reduce  the 
inventory  of  O  atoms  and  so  eliminate  a  quencher  of  the  upper 
laser  level,  I*.  At  low  pressures  where  the  3-body  density  is 
low,  the  rate  of  reaction  of  O  with  N02  occurs  at  a  higher 
rate  than  with  NO.  Hence  addition  of  N02  in  the  downstream 
region  may  be  preferred  over  that  of  NO  if  the  residence  time  in 
the  flow  tube  is  a  limiting  factor.  These  reactions  are,  however, 
exothermic  and  so  can  increase  the  gas  temperature  which  is 
generally  not  beneficial. 

Laser  gain  is  ultimately  achieved  by  injection  of  I2 
downstream  of  the  plasma  zone,  and  by  its  reacting  with 
O,  02(I  A),  and  02('E)  to  create  I  atoms  and  to  pump  the 
upper  laser  level.  In  a  conventional  COIL  there  is  a  negligible 


inventory  of  O  atoms  and  02(J  E),  and  so  the  dissociation  of 
I2  is  dominantly  by  02(*  A)  in  a  two-step  process, 

02('A)  +  I2  -»  02  +  i;,  (Jfc  =7  x  KTl5cm3s_l),  (20) 
02('A)  +  i;  -►  02  +  I  +  I,  (t  =  3x  10-locm3s_l).  (21) 

In  eCOIL,  the  presence  of  02(*  E)  and  O  helps  in  dissociating 
I2  and  producing  I  atoms.  A  reaction  intermediate  10  is  also 
helpful  in  this  regard, 

02(’ E)  +  12  ->  02  +  I  +  I,  (*  =  2.8  x  10_llcmV), 

(22) 

0  +  I2->I0  +  I,  (k  =  1.4  x  lO'^cm’s"1),  (23) 

0  +  I0-xI  +  02,  (k  =  1.4  x  10-|0cm3s-1).  (24) 

The  laser  pumping  reaction  by  02('A)  in  collisions  with  I 
is  favoured  at  lower  temperatures  over  its  endothermic  back 
reaction, 

o2('a)  +  i-x  o2  +  r, 

(. k  =  7.7  x  10-"(rg/300)_l  cm3  s"1),  (25) 

o2  +  r  -+  o2('a)  +  i, 

(k  =  1.0  x  10~lo(7'g/300rle~403/T«).  (26) 

To  suppress  the  back  reaction,  the  laser  cavity  is  typically 
placed  in  a  supersonically  expanded  flow  to  lower  the 
translational  temperature.  As  such,  the  threshold  yield,  T,h,  of 
02(‘ A)  required  for  positive  optical  gain  in  an  undissociated 
flow  of  02  is  [  1 7], 

KM1  A)]  1 

y  _  1  7 3  _  _  ('ll} 

[02]  “  1  +  1.5exp(401/7g)' 
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For  example,  in  an  undissociated  flow  of  O2,  the  threshold 
yield  at  room  temperature  is  15%  whereas  at  180K  the 
threshold  yield  is  6%. 

Quench  ng  of  I*  occurs  dominantly  by  reactions  with  O 
and  O2  (the  backward  reaction),  and  to  a  lesser  extent  with 
02(*  A)  and  NO, 

o  +  r  -» 0  +  1, 


0 k  =  8.0  x  10~12  cm3  s_l). 

(28) 

NO  +  I*  -*•  NO  +  I, 

(k  =  1.2  x  10‘l3cm3s~'). 

(29) 

NO  and  NO-/  also  help  in  removing  ground  state  I  atoms  and 
so  aid  in  maintaining  the  inversion, 

I  +  NO  +  O2  -+  IN0  +  02i  ( k  =  1.6  x  10-32  cm6  s_l), 

(30) 

I  +  NO  +  He  -+  INO  +  He,  (*  =  6.0  x  10"33  cm6  s_1), 

(31) 

I  +  N02  +  O.  -+  1N02  +  02,  (*  =  2.6  X  10"31  cm6  s_l), 

(32) 

I  +  NO2  +  HC-+  INOj  +  He,  (k  =  6.0  x  10"33 cm6 s_l). 

(33) 

The  INOj,  species  in  turn  abstract  I  atoms  to  release  I2  and  NO* 
back  into  the  flow, 

INO  +  I  — »  I2  +  NO,  (k  =  1.6  x  10"locmV),  (34) 
IN02  +  I  -+  I2  +  NOj,  (*  =  8.3  x  l<rllcmV).  (35) 

Typically,  the  injection  of  I2  occurs  during  or  following  a 
supersonic  expansion  of  the  gas  to  velocities  equivalent  to 
Mach  2  or  M  ich  3  to  lower  Tg  to  decrease  the  yields  of  A) 
required  to  achieve  positive  gain.  Accurately  simulating  a 
supersonic  expansion  is  difficult  in  our  modelling  platform. 
In  order  to  piovide  a  best  case  estimate  for  laser  gain  we  used 
for  the  temperature  of  the  reactions  in  equations  (25)  and  (26) 
a  value  appropriate  for  a  Mach  2  flow.  This  temperature  can 
be  approximated  from 


where  T0  is  the  stagnation  (or  tank)  temperature,  y  is  the  ratio 
of  specific  heats  and  M  is  the  Mach  number.  For  M  =  2 
and  Tq  =  30)  K,  this  ratio  is  ^2.18,  and  would  lead  to  a  gas 
temperature  of  137  K.  The  reaction  mechanism  was  validated 
by  comparing,  with  experimental  data  from  Carroll  et  al  [34], 
as  discussed  l?elow. 

4.  Consequences  of  NO  in  the  inlet  flow 

A  schematic  of  the  cylindrically  symmetric,  6  cm  diameter  flow 
tube  used  in  tlis  study  is  shown  in  figure  1 .  In  the  base  case,  a 
power  of  40  W  was  capacitively  coupled  using  ring  electrodes 
operated  at  25  MHz.  A  mixture  of  He/02/NO  at  3  Torr  entered 
through  the  i  ilet  at  a  flow  rate  of  6  slpm  which  corresponds 
to  an  axial  speed  of  ^890  cm  s-1.  The  flow  consisted  of  30% 
O2  with  the  balance  divided  between  He  and  NO.  The  NO 
mole  fraction  was  varied  from  0-10%.  Two  injection  nozzles 
are  located  downstream.  (In  the  2D  cylindrically  symmetric 
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Figure  1.  Geometry  of  the  cylindrically  symmetric  reactor.  The 
flow  enters  from  the  top  and  is  pumped  from  the  bottom.  Discharge 
power  at  25  MHz  is  capacitively  coupled  through  ring  electrodes. 
Two  nozzles  downstream  of  the  discharge  inject  mixtures  of  He, 
NO,  N02  and/or  I2. 


geometry  these  nozzles  appear  to  be  rings.)  The  first  nozzle 
at  51.5  cm  was  used  to  inject  a  mixture  of  He/N0/N02-  The 
second  nozzle  at  64.5  cm  injected  a  mixture  of  He/I2.  At  the 
outlet,  axial  gradients  were  assumed  to  be  zero  and  the  exit 
speeds  were  adjusted  to  maintain  constant  pressure  and  mass 
flux. 

Plasma  characteristics  (nc,  TCy  negative  total  negative  ion 
density  M total  positive  ion  density  M*  and  power  density) 
obtained  with  nonPDPSIM  are  shown  in  the  vicinity  of  the 
electrodes  in  figure  2  for  the  base  case  having  an  inlet  flow 
of  He/02/NO  =  67/30/3.  Tc  in  the  bulk  of  the  discharge 
was  2.3  eV  and  the  peak  electron  density  was  9.0  x  109cm-3. 
Since  NO  and  O2  are  attaching  species,  the  negative  ion  density 
(maximum  of  1.1  x  1010  cm-3)  is  commensurate  with  nc.  The 
electron  density  is  fairly  symmetric  between  the  electrodes 
and  relatively  uniform.  The  power  deposition  is  moderately 
higher  near  the  upstream  electrode  due  to  gas  heating  which 
reduces  the  neutral  densities  downstream.  This  is  not  a  general 
result  as  higher  power  deposition  and  more  rarefaction  can 
produce  regions  of  locally  intense  power  deposition  near  the 
downstream  electrode. 

The  densities  of  the  neutral  species  O,  C^1  A),  I  and  I*  and 
Tg  for  the  base  case  are  shown  in  figure  3.  A  flow  of  36  seem 
of  pure  NO  was  injected  through  the  first  nozzle  and  a  flow 
of  100  seem  of  He/l2  =  99/1  was  injected  through  the  second 
nozzle.  Having  few  quenchers  at  these  pressures,  the  C^1  A) 
accumulates  as  the  gas  flows  through  the  discharge  reaching 
a  maximum  value  of  1.35  x  10l5cm~3.  The  (^(‘A)  density 
remains  nearly  constant  thereafter  until  the  injection  point  for 
I2.  Reactions  of  C^1  A)  with  I2  and  I  (the  latter  being  the  laser 
pumping  reaction  producing  I*)  reduce  its  density  by  a  factor 
of  five  to  2.7  x  1014  cm-3. 
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Figure  2.  Ba»e  case  plasma  characteristics  from  the  2D  model  for 
3Torr,  40  W.  He/02/N0  =  67/30/3  and  a  flow  rate  of  6slpm.  (a) 
Electron  density,  ( b )  electron  temperature,  (c)  sum  of  negative  ion 
densities,  ( d )  sum  of  positive  ion  densities  and  (e)  power  density. 
The  scales  arc  linear  for  Tt  and  power  (0  to  maximum),  and  the 
densities  are  plotted  on  a  2-decade  log  scale. 
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Similar  to  CM1  A),  the  density  of  O  atoms  accumulate 
from  electron  impact  dissociation  of  O2  passing  through  the 
discharge  (with  a  small  amount  of  depletion  of  the  O  in  forming 
O3)  until  the  injection  point  of  NO  at  the  first  nozzle.  The  O 
densities  decrease  from  a  peak  value  of  1.4  x  1015  to  4.0x 
1014  cm"3  downstream  of  the  NO  injection  due  to  its  reaction 
with  NO.  The  remaining  O  atoms  flow  to  the  injection  point  of 
I2  where  the)  are  further  depleted  in  dissociating  reactions  with 
I2.  The  I2  is  nearly  completely  dissociated  by  reactions  with 
the  O  atoms  to  form  I.  Pumping  reactions  between  02(!A) 
and  1  produ;e  I*.  The  region  of  positive  gain  (where  the 
inversion  de  isity  G  =  [/"]  -  l/2[/]  >  0)  is  a  narrow  band 
downstream  of  the  I2  injection  point,  with  a  peak  value  of 
G  =  2.6  x  10n  cm"3. 

7g  has  two  local  maxima.  The  first  is  due  to  discharge 
Joule  heating,  Frank-Condon  heating  and  exothermic 
reactions  of  NO  with  O  reaching  340  K  adjacent  to  the 
downstream  electrode.  As  the  walls  are  held  at  300  K,  the  gas 
rapidly  cools  by  thermal  conduction.  A  second  local  maximum 
occurs  downstream  of  the  first  injection  nozzle  (up  to  348  K) 
due  to  additional  exothermic  reactions  of  O  and  the  newly 
injected  NO. 

To  investigate  the  consequences  of  NO  in  the  inlet  flow 
over  a  wider  parameter  space,  the  plug  flow  model  GlobalK/N 
was  used.  We  first  addressed  the  importance  of  NO-  and  NOJ 
in  the  reaction  mechanism  and  their  effects  on  the  densities  of 
electrons,  O?  (*  A)  and  O2 ( 1 5D ) .  The  densities  of  these  species 
arc  shown  along  the  axis  of  the  discharge  in  figure  4  for  a 
3Torr  mixture  of  He/02/NO  =  60/30/10  with  40  W  power 
deposition.  A  high  mole  fraction  of  NO  was  used  as  an  extreme 


Figure  3.  Neutral  species  densities  and  Tt  for  the  base  case.  36  seem 
of  pure  NO  was  injected  through  the  first  nozzle  and  100  seem  of 
He/I2  =  99/ 1  was  injected  through  the  second  nozzle.  Densities  of 
(a)  02(‘  A),  (b)  O,  (c)  I,  ( d )  P,  (e)  inversion  density  ((I*]-0.5[I])  and 
(/)  gas  temperature.  The  scales  are  linear  except  for  02(I  A)  and  O 
which  are  plotted  on  a  2-decade  log  scale.  Injection  of  NO  decreases 
the  flow  of  O  and  injection  of  I2  consumes  02 (*  A).  Exothermic 
reactions  at  both  injection  points  produce  local  maxima  in  T%. 

case.  When  including  NO-  and  NOJ  the  peak  electron  density 
decreased  approximately  10%  from  1.4  x  10,0cm-3.  The 
reduction  in  nc  results  in  large  part  from  the  attachment  of 
electrons  to  NO  through  3-body  reactions  similar  to  those  for 
OJ  formation  (equation  (7))  and  dissociative  attachment  to 
NO.  The  changes  in  densities  of  02(,A)  and  02(*£)  were 
relatively  small  (<5%)  and  could  be  attributed  to  the  fact 
that  nc  decreases  and  Te  increases  modestly  (0.05-0. 1  eV  from 
%2.0eV)  with  the  inclusion  of  NO“  and  NO^.  (Recall  that 
electron  impact  excitation  of  02( 1  A)  is  maximum  for  Tc  =  1- 
1 .5  eV.)  The  density  of  O  atoms  was  relatively  independent  on 
the  inclusion  of  NO-  and  NOJ .  Based  on  these  trends,  we  can 
expect  results  from  nonPDPS/M  (which  do  not  include  NO” 
and  NOj  in  the  reaction  mechanism)  to  over-predict  02(' A) 
densities  by  a  few  per  cent. 

The  consequences  of  NO  in  the  inlet  flow  on  the  maximum 
densities  of  charged  species,  Tgl  Tc,  and  power  deposition  into 
different  species  by  electron  impact  are  shown  in  figure  5  for 
3  Torr  and  40  W.  A  modest  increase  in  nc,  from  1.1  x  10locm-3 
to  1.25  x  10locm"3,  occurs  as  the  NO  mole  fraction  is 
increased  from  0  to  10%  due  to  the  higher  rates  of  ionization 
with  NO  in  spite  of  a  decrease  in  the  positive  and  negative 
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Figure  4.  Dej>endence  of  plasma  characteristics  on  the  inclusion  of 
NO"  and  NO;.  for  3Torr,  40  W,  He/02/N0  =  60/30/10  and  6slpm. 
Densities  of  (<i)  electrons,  ( b )  O^1  A)  and  (c)  O7C  E).  The  densities 
of  O  did  not  slow  any  noticeable  variation.  These  results  are  from 
the  plug  flow  model.  Inclusion  of  these  negative  ions  produces 
nominal  changes  in  excited  state  densities. 


ion  densities.  The  increase  in  nc  results  from  the  density  of 
negative  ions  decreasing  more  rapidly  than  the  positive  ions 
when  increasing  the  mole  fraction  of  NO.  The  positive  ion 
density  decreases  from  2.3  x  10locm“3  to  2.0  x  10,0cm”3 
while  the  decrease  in  negative  ion  density  is  1.2  x  10,0cm-3 
to  0.75  x  10,0cm-3.  Tc  decreases  from  2.2  to  1.9eV  as  the 
NO  mole  fraction  is  increased  from  0  to  10%.  The  reduction 
in  Te  is  due  to  the  additional  power  loss  resulting  from  electron 
impact  on  NO  when  substituting  NO  for  He  and  the  reduction  in 
the  average  i  Dnization  potential  thereby  enabling  the  discharge 
to  be  self-sustained  with  a  lower  value  of  Te.  This  trend 
is  partially  offset  by  the  additional  attachment  probability 
represented  by  NO. 

While  the  total  ion  densities  gradually  decrease  with 
increasing  NO  mole  fraction,  the  densities  of  individual  ions 


Figure  5.  Plasma  and  gas  characteristics  while  varying  the  NO  mole 
fraction  for  3Torr,  40  W  and  a  6slpm  flow  of 
He/02/NO  =  70  —  x/30/x.  (a)  Maximum  nc,  ( b )  maximum 
positive-  and  negative-ion  densities,  (c)  maximum  Tg  and  ( d ) 
fractional  power  deposition  into  02  and  NO,  and  in  producing 
02(‘ A)  and  02(*E).  These  results  are  from  the  plug  flow  model 
and  from  calculations  of  EEDs.  In  spite  of  NO  consuming  a  large 
fraction  of  the  power,  the  decrease  in  Tc  increases  the  rate  of 
excitation  of  02(*  A)  and  so  there  is  not  a  large  change  in  power 
dissipated  in  exciting  02(l  A). 

show  more  variation.  In  the  absence  of  NO,  OJ  forms  the 
majority  (>99%)  of  the  positive  ions.  However  with  as  little 
as  1-2%  of  NO,  the  majority  ion  is  NO+.  This  is  due  to 
charge  exchange  reactions  between  OJ  and  NO.  The  dominant 
negative  ion  in  the  absence  of  NO  is  O-,  and  as  the  NO  mole 
fraction  is  increased,  the  density  of  NOJ  ions  increases,  a  result 
of  charge  exchange  between  O-  and  NO2  (equation  ( 1  7)). 

Tg  increases  with  NO  mole  fraction  primarily  due  to 
the  exothermicity  of  reactions  between  O,  NO  and  NO2 
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(equations  ( I 8)  and  ( 1 9)).  In  the  absence  of  NO,  the  maximum 
value  of  Tg  is  328  K  (increasing  from  the  inlet  value  of  300  K). 
This  value  increases  to  362 K  with  10%  NO.  The  peak  gas 
temperatures  occur  just  downstream  of  the  discharge  region. 
With  large  NO  mole  fractions,  the  contribution  of  Frank- 
Condon  gas  heating  due  to  dissociation  of  NO  can  be  more 
than  a  few  per  cent. 

The  fractional  power  deposited  into  NO  increases  nearly 
linearly  witi  the  NO  mole  fraction,  increasing  to  25%  for 
an  NO  mole  fraction  of  10%.  (Note  that  these  values  were 
computed  using  the  actual  mole  fractions  of  species  in  the 
discharge  and  not  the  inlet  conditions.)  Correspondingly,  the 
power  depoj  ition  into  O2  decreases,  from  more  than  85%  in  the 
absence  of  NO  to  66%  for  an  NO  mole  fraction  of  1 0%.  (Little 
power  is  dissipated  by  He.)  In  spite  of  the  fraction  of  discharge 
power  dissipated  in  O2  decreasing,  the  power  expended  in 
exciting  02( 1  A),  cr,  was  relatively  constant,  reaching  a  shallow 
maximum  of  8.7%  for  a  NO  mole  fraction  of  5%.  As  the  NO 
mole  fractio  1  increases,  Tc  decreases  from  2.2  to  1 .9  eV  thereby 
increasing  the  efficiency  of  02(*A)  production  by  electron 
impact  [  1 1 J.  Fortuitously,  a  combination  of  decreasing  power 
dissipated  by  O2  and  an  increase  in  the  efficiency  of  exciting  of 
02(‘  A)  maintains  the  total  power  expended  in  exciting  02^  A) 
relatively  constant.  This  is  not  necessarily  a  general  result  but 
lowering  Tc  is  generally  beneficial.  Similar  trends  are  seen  for 
O2  ( 1 E)  as  the  fractional  power  expended  in  exciting  02(!  E) 
increases  slightly  from  3.2%  to  3.9%  as  the  NO  mole  fraction 
increases  to  10%. 

The  consequences  of  NO  mole  fraction  in  the  inlet  flow  on 
the  densities  of  O,  O2  ( 1  A)  and  O2  ( 1 E ),  and  the  yield  of  O2  ( 1  A) 
along  the  axis  of  the  discharge  are  shown  in  figure  6  for  the 
base  case  of  3Torr,  40  W  and  30%  O2.  (Note  that  this  yield 
is  based  on  he  mole  fraction  equivalent  of  02  for  all  oxygen 
species  (II])  A  flow  of  100  seem  of  Hc/l2  =  99/1  is  injected 
through  the  second  nozzle.  In  the  absence  of  NO,  the  density 
of  02(!  A)  increases  to  1 .2  x  1015  cm“3due  to  electron  impact 
excitation  as  the  inlet  flow  passes  through  the  plasma  zone.  The 
increase  in  the  density  of  02(!  A)  to  1 .75  x  1015  cm-3  after  the 
discharge  zene  results  largely  from  quenching  of  a  density  of 
0.4  x  1015  cm-3  of  O2 ( *  E)  by  O  atoms  and  excitation  transfer 
from  0('D).The  density  of  02(*  A)  upstream  of  the^  injection 
point  monotDnically  decreases  with  the  addition  of  NO,  from 
1.75  x  10,5cm”3  without  NO  to  1.1  x  10l5cm-3  with  10% 
NO.  A  portion  of  this  decrease  results  from  an  increase  in  the 
Tg  when  adding  NO  though  the  yield  ofC^A)  also  decreases 
from  5.9%  when  adding  NO. 

The  production  of  C^C1  A)  is  dominated  by  electron  impact 
on  O2  and  so  is  dependent  on  Tt  and  nc.  Since  the  power 
deposition  producing  (^(‘A)  is  nearly  constant,  secondary 
processes  m  jst  be  responsible  for  the  decrease  in  the  yield  of 
02(’  A)  with  increasing  NO.  One  such  secondary  process  is  the 
production  of  02(' A)  by  excitation  transfer  between  0(*D) 
and  02.  Thi  majority  of  0(‘D)  is  produced  by  dissociative 
excitation  of  O2  (equation  (5)).  Reactions  of  0(‘D)  with  O2 
produces  its  excited  states, 

0(,D)  +  02  -►  0  +  02('E), 

(*  =  2.6  x  10“M  exp(67/Tg) cm3  s_l),  (37) 

0(‘D)  +  02  0  +  02(‘A), 

( k  =  1.5  x  10" 12  exp(67/Tg)cm3  s-1).  (38) 
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Figure  6.  Densities  of  oxygen  species  and  yield  of  O2  ( 1  A)  along 
the  axis  of  the  tube  while  varying  NO  mole  fraction  for  3  Torr  and 
40  W  for  a  6slpm  flow  of  Hc/02/NO  =  70  -  x/30/x.  100  seem  of 
He/I2  =  99/1  is  injected  through  the  second  nozzle,  (a)  02(*  A),  ( b ) 
yield  of  02('  A),  (c)  O  and  ( d )  02 (*  E).  These  results  are  from  the 
plug  flow  model.  Increasing  flows  of  NO  decrease  the  flow  of  O 
atoms. 
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With  the  decrease  inTc  with  addition  of  NO,  the  production 
of  O('D)  by  electron  impact  dissociation  of  O2  decreases 
and  so  these  secondary  sources  of  02(*A)  and  02(,£)  also 
decrease.  Furthermore,  both  NO  and  NO2  are  quenchers 
ofO('D), 


NO  +  CX'D)  -V  NO  +  O, 

(*  =  4.0  X  10_n  cm3  s-1), 

(39) 

N02  +  0('C) 

no  +  o2. 

(k  =  3.0  x  10"locmV), 

(40) 

N02  +  0('C) 

-►  no2  +  0, 

(k  =  3.2  x  10“ 10  cm3  s_l). 

(41) 

So  the  addition  of  NO  reduces  the  production  of  O('D)  and 
increases  its  rate  of  quenching,  thereby  reducing  the  production 
of  02(‘ A)  oy  excitation  transfer.  The  importance  of  the 
quenching  reactions  of  O('D)  is  demonstrated  by  excluding 
the  reaction?  in  equations  (39)— <4 1 )  from  the  mechanism.  For 
10%  NO  in  the  flow,  the  yield  of  02(!A)  improved  from 
4.6%  to  5.3%  when  quenching  of  O('D)  is  eliminated.  The 
quenching  of  02(jA)  by  NO  has  a  small  rate  coefficient 
(3.5  x  10“n  cm3  s"1)  and  so  does  not  significantly  contribute 
to  the  loss  of  02(‘  A). 

The  density  of  O  atoms  increases  monotonically  through 
the  discharge  zone  to  a  maximum  value  of  2.5  x  10,5cm~3 
in  the  absence  of  NO.  (The  small  increase  in  the  density  of 
O  after  the  c  ischarge  is  mainly  due  to  gas  cooling.)  With  the 
addition  of  140,  the  density  of  O  atoms  decreases  throughout 
the  flow-tube  and.  in  particular,  downstream  of  the  discharge. 
This  decrease  is  due  to  both  the  reduction  in  the  rate  of 
electron  impact  dissociation  of  O2  by  the  decrease  in  Tc  and 
the  formation  of  NO2  in  reactions  with  NO.  In  the  post- 
discharge  region  in  the  absence  of  NO,  the  majority  of  02('  £) 
is  converted  to  02(‘  A)  through  collisions  with  O  and  O3.  Due 
to  the  reduciion  of  O  atoms  with  increasing  NO,  the  rates  of 
quenching  of  02(,E)  to  02('A)  are  also  smaller,  leading  to 
higher  densities  of  02(‘  E)  in  the  afterglow.  Higher  O2 ( 1 ) 
densities  are  not  necessarily  bad  since  they  help  in  dissociating 
I2  but  maintaining  those  densities  does  result  in  lower  densities 
of  02('  A)  that  directly  pump  P. 

The  densities  of  I  and  I*,  and  the  optical  gain  at  1.315 /x 
are  shown  in  figure  7  for  the  conditions  of  figure  6  ( 1  seem 
of  I2  injecte  d  in  a  100  seem,  He/l2  =  99/1  flow).  The 
gain  was  given  by  or([P]  —  0.5(1])  where  a  is  the  stimulated 
emission  cross-section.  At  pressures  of  less  than  tens  of  Torr, 
Doppler  broadening  dominates  over  pressure  broadening,  and 
so  the  stimulated  emission  cross-section  can  be  approximated 
by  [35], 

a  =  1.33  x  10-I67--|/Jcm2.  (42) 

In  the  absence  of  NO,  the  density  of  I  increases  from 
1 .75  x  1013  :m“3  at  the  I2  injection  point  to  3.2  x  1013  cm-3 
downstream  due  in  large  part  to  the  reaction  of  O  atoms  with 
I2.  The  density  of  I*  is  maximum  at  the  injection  point  at 
1 .45  x  1013  cm-3  and  decreases  to  negligible  values  by  the  end 
of  the  flow  tube  due  to  the  depletion  of  02(‘ A)  (the  species 
responsible  for  the  pumping  reaction)  and  quenching  by  O 
atoms.  Since  the  lifetime  of  P  (125/xs  for  quenching  by  O) 
is  short  com  Dared  with  flow  times,  the  density  of  I*  does  not 
appreciably  accumulate  in  the  discharge  and  its  density  is  a 


b)  Flow  Distance  (cm) 


c)  Flow  Distance  (cm) 


Figure  7.  Densities  of  atomic  iodine  species  as  a  function  of  inlet 
NO  mole  fraction  for  3  Torr,  40  W,  HeAI^/NO  =  70  -  jc/30/jc  and 
6slpm.  100  seem  of  He/I2  =  99/1  is  injected  through  the  second 
nozzle.  Densities  of  (a)  I,  (b)  I*  and  (c)  optical  gain.  These  results 
arc  from  the  plug  flow  model.  Low  flows  of  NO  produce  the 
maximum  peak  gain  whereas  high  flows  of  NO  produce  larger 
plumes  of  positive  gain. 

reflection  of  instantaneous  production  and  quenching  rates. 
The  end  result  is  a  peak  gain  of  4.2  x  10“5cm_l  within  a 
centimetre  of  the  I2  injection  point.  Note  that  the  rate  of 
quenching  of  I*  by  O  is  faster  than  that  due  to  spontaneous 
emission  (0.1  s)  [33]. 

With  injection  of  I2  the  density  of  02('  A)  decreases  by 
virtue  of  excitation  transfer  and  dissociative  excitation  of  I2, 
and  reactions  with  I  which  pumps  I*.  In  the  absence  of  NO, 
the  density  of  02('A)  is  fully  depleted  by  the  reactions.  In 
the  absence  of  NO  there  is  also  a  large  density  of  O  atoms 
at  the  injection  point  which  react  with  I2  producing  IO  and  I. 
Since  the  rate  of  dissociation  of  I2  by  02(!A)  and  02(,£)  is 
slower  than  by  O  atoms  there  is  a  larger  density  of  I  available 
for  O2  ( 1 A )  to  react  with,  and  so  the  O2  ( 1 A )  is  rapidly  depleted. 
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As  the  flow  of  NO  is  increased,  the  flow  of  O  atoms  decreases, 
resulting  in  <i  lower  rate  of  dissociation  of  I2  and  fewer  I  atoms. 
With  the  lov/er  density  of  I  the  reactivity  of  02(!  A)  is  lower; 
and  so  its  density  decreases  less  rapidly  after  injection  of  I2. 

The  production  of  P  results  from  reactions  of  02(IA) 
with  I  while  the  quenching  of  P  is  largely  due  to  collisions 
with  O  atons.  Increasing  the  NO  mole  fraction  decreases 
the  O  atom  density  downstream  thereby  lowering  the  rates  of 
quenching  of  I*  and  extending  its  plume  beyond  the  I2  injection 
point.  However,  having  too  large  an  NO  mole  fraction  results  in 
too  low  rate*  of  dissociation  of  I2  and  hence  poor  utilization  of 
O2 ( 1  A)-  Th'2  end  result  is  that  the  region  over  which  positive 
gain  can  be  sustained  is  maximum  for  an  intermediate  mole 
fraction  of  NO  of  3%. 

The  purpose  of  flowing  NO  (or  injecting  NO2)  is  largely  to 
manage  the  O  atom  density  and  so  control  the  quenching  of  I* 
by  O.  A  sens  itivity  study  was  conducted  of  the  rate  coefficient 
for  this  quenching  reaction  from  the  nominal  value  of  8  x 
10_,2cm3s'  ‘.  The  densities  of  02(*A)  and  P,  and  gain  are 
shown  in  figure  8  while  varying  this  rate  coefficient  from  0  to 
1 .0  x  1 0~ 1 1  <  m3  s“ 1 .  In  the  absence  of  quenching,  the  densities 
of  P  and  O2  (‘  A)  do  not  significantly  decrease  downstream  of 
the  I2  injection  point.  In  the  absence  of  quenching  by  O,  the 
predominam  quencher  of  P  is  O2  and  this  quenching  produces 
02(,A).  I’he  forward  and  backward  pumping  reactions 
(equations  (25)  and  (26))  reach  an  equilibrium  where  the  P  and 
02(l  A)  denudes  gradually  decrease  due  to  minor  quenchers 
of  P  (such  as  NO  and  02(‘  A))  and  radiative  relaxation.  When 
increasing  the  rate  coefficient  for  quenching  P,  the  density  of 
P  decreases  proportionately.  By  removing  I*  in  this  manner 
the  rate  of  the  backward  reaction  with  O2  decreases  and  so 
does  the  density  of  02(,A).  The  extent  of  positive  gain  is 
progressively  limited  to  the  few  cm  beyond  the  I2  injection 
point  as  the  quenching  of  P  increases,  though  the  peak  value 
of  gain  is  not  particularly  sensitive  to  the  rate  coefficient. 

The  densities  of  02(!  A)  and  O,  and  Tg  2  cm  upstream  of 
the  second  nazzle  are  shown  in  figure  9  as  a  function  of  power 
deposition  a  id  NO  mole  fraction  in  the  inlet  flow.  In  general, 
Tg  increases  with  power  deposition  and  with  NO  addition 
reaching  a  maximum  of  480  K  with  400  W  power  deposition 
and  10%  NO  in  the  inlet  flow.  Increasing  power  deposition 
produces  more  electron  impact  dissociation  of  O2  and  NO, 
producing  larger  densities  of  O  atoms.  However,  increasing 
NO  mole  fractions  decreases  the  density  of  O  atoms  by  virtue 
of  scavenging  by  NO  and  NO2  and  decreasing  7C. 

The  dersity  of  02(*A)  decreases  with  NO  addition,  as 
discussed  above,  and  increases  with  power  deposition.  The 
saturation  in  the  density  O2 ( 1  A)  at  6.5  x  1015  cm-3  at  higher 
powers  is  due  in  part  to  the  depletion  of  O2  by  electron  impact 
dissociation  and  in  part  to  gas  heating.  For  example,  with 
400  W  and  0%  NO,  the  fractional  dissociation  of  O2  is  68%. 
The  addition  of  NO  reduces  the  depletion  of  O2  by  both 
reducing  the  rate  of  O2  dissociation  and  by  recycling  O  atoms 
back  to  O2.  For  example,  for  10%  NO  addition  and  400  W,  the 
fractional  dissociation  decreases  to  36%.  The  yield  of  02(*  A) 
saturates  with  power  at  18%  due  to  the  depletion  of  O2. 

The  densities  of  I  and  P,  and  gain  as  function  of  power 
deposition  and  NO  mole  fraction  are  shown  in  figure  1 0  at  a 
location  of  2  cm  downstream  from  the  I2  injection  point.  This 
mixing  length  of  2  cm  was  chosen  based  on  previous  studies 


Figure  8.  Sensitivities  of  the  value  of  the  rate  coefficient  for  the 
quenching  reaction  between  O  and  P  on  the  post-discharge  kinetics 
for  3  Torn,  40  W,  He/02/NO  =  69/30/1  and  6slpm.  100  seem  of 
He/I2  =  99/1  is  injected  through  the  second  nozzle.  Densities  of  (a) 
02('  A),  ( b )  I*  and  ( c )  optical  gain.  These  results  are  from  the  plug 
flow  model. 


for  similar  flow  conditions  [36].  The  inlet  flow  has  30%  02 
and  the  flow  injected  through  the  second  nozzle  is  1  seem  of  I2 
in  a  1 00  seem  flow  of  He/l2  =  99/ 1 .  At  low  power  deposition, 
the  density  of  I  is  large,  (2.2-2. 3)  x  10I3cm“3,  because  the 
yields  of  CM1  A)  are  low  enough  that  the  pumping  of  I  to  P 
is  slow.  As  the  yield  of  C^A)  increases  at  higher  powers, 
the  pumping  reactions  reduce  the  density  of  I  and  increase  that 
of  P,  leading  to  an  increase  in  gain.  However,  at  large  power 
deposition,  the  density  of  O  increases  whereas  that  of  C^1  A) 
saturates.  This  leads  to  an  increased  rate  of  quenching  of  P  by 
O  which  reduces  the  density  of  P  and  increases  that  of  I.  The 
end  result  is  a  reduction  in  gain. 

Increasing  NO  in  the  inlet  flow  reduces  the  density  of 
O  atoms  at  the  injection  point,  thereby  reducing  the  rate  of 
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a)  Power  (W) 


d)  Power  (W) 


Figure  9.  Coisequences  of  power  deposition  and  NO  mole  fraction 
on  the  neutral  gas  properties  for  3Torr,  He/02/NO  =  70  —  jc / 30/ jc 
and  6slpm.  (<i)  7*g,  (fc)  density  of  O,  (c)  density  of  02('  A)  and  ( d ) 
yield  of  02('  A).  CM1  A)  densities  saturate  with  power  due  to  a  high 
degree  of  dissociation  of  O2.  These  results  are  from  the  plug  flow 
model  and  the  values  are  for  2  cm  upstream  of  the  second  nozzle. 

dissociation  of  I2  and  the  density  of  I  atoms.  The  reduction 
in  O  by  addition  of  NO  also  reduces  the  quenching  of  I*  by 
O,  making  the  gain  predominantly  dependent  on  the  yield  of 
O2  (*  A).  Asa  result,  for  3%  NO,  the  gain  decreases  from  6.0 


Figure  10.  Consequences  of  power  deposition  and  NO  mole 
fraction  on  iodine  species  and  gain  for  3Torr, 

He/02/NO  =  70  —  jc/30/jc  and  6slpm.  Densities  of  (a)  I  and  ( b )  I* 
and  (c)  optical  gain.  These  results  are  from  the  plug  flow  model  and 
the  values  are  for  2  cm  downstream  of  the  second  nozzle.  Gain  is 
maximum  at  high  power  only  with  high  flow  rates  of  NO. 

x  10-5  to  4.3  x  10_5cm_l  between  150  and  400  W,  whereas 
in  the  absence  of  NO,  the  gain  reduces  by  nearly  a  factor  of  3. 

The  eCOIL  system  may  operate  in  either  power  limited  or 
iodine  limited  modes.  In  the  iodine  limited  mode,  the  flow 
of  02(,A)  and  O  generated  by  the  discharge  fully  utilizes 
the  injected  flow  of  I2,  and  so  gain  saturates  with  increasing 
power.  To  some  degree  (quenching  of  I*  and  depletion  of  02 
aside),  this  is  the  mode  that  applies  to  the  results  shown  in 
figure  10.  In  the  power  limited  mode,  the  flow  of  C^C1  A)  and 
O  is  insufficient  to  fully  utilize  the  injected  flow  of  I2,  and  so 
gain  saturates  with  flow  of  I2. 

These  modes  of  operation  are  demonstrated  by  varying 
power  deposition  and  I2  flow  rate.  The  densities  of  I  and  I*, 
and  gain  are  shown  2  cm  downstream  of  the  I2  injection  point 
in  figure  II  for  a  6slpm  inlet  flow  of  He/02/NO  =  67/30/3 
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Figure  11.  Iodine  species  densities  and  gain  while  varying  the  I2 
flow  rate  through  the  second  nozzle.  Densities  of  (a)  I  and  ( b )  I\ 
and  (c)  optical  gain  for  3%  NO  in  the  inlet  flow  and  ( d )  gain  for 
80  W  power  deposition.  These  results  are  from  the  plug  flow  model 
and  the  value;  are  for  2  cm  downstream  of  the  second  nozzle. 
Varying  power  and  I2  flow  rate  switches  between  the  power  limited 
and  I2  limited  regimes. 


while  varying  power  and  I2  flow  rate.  The  densities  of  I  and  I* 
increase  nearly  linearly  with  the  flow  of  I2  and  increase  sub- 
linearly  with  power  deposition.  For  example,  with  a  flow  of 
3  seem  of  I2,  the  density  of  I  is  4.6  x  10,3cm-3  with  40  W 


and  only  5.6  x  1013  cm"3  with  120  W.  For  these  conditions  the 
production  of  1  is  limited  by  the  availability  of  I2  as  the  injected 
I2  is  nearly  completly  dissociated  by  reactions  with  O,  A) 
and  O2  ( *  E )  for  all  power  depositions.  The  production  of  P  is 
limited  by  the  yield  of  C^1  A)  which  does  not  scale  linearly 
with  power  at  higher  power  depositions.  As  a  consequence  the 
density  of  I*  saturates  with  power  deposition  with  large  flows 
of  I2.  This  saturation  is  exacerbated  at  high  power  depositions 
(or  low  NO  mole  fractions)  by  large  densities  of  O  which 
quench  P. 

With  a  power  deposition  of  40  W,  gain  increases  with  flow 
rate  of  I2  (shown  in  figure  1 1(c))  saturating  at  3  x  10-5  cm-1 
at  2  seem,  a  consequence  of  being  in  a  power  limited  regime. 
As  the  power  increases  the  gain  increases,  and  the  I2  flow  rate 
at  which  the  maximum  gain  is  obtained  also  increases.  For 
example,  for  80  W,  the  gain  saturates  at  9  x  10-5cm_l  for 
an  I2  flow  rate  of  3  seem.  The  higher  power  deposition  is 
able  to  better  utilize  the  increased  flow  of  I2.  As  the  power 
further  increases,  the  peak  gain  increases  only  moderately  and 
decreases  above  200  W  as  the  system  transitions  to  an  I2  limited 
regime  and  the  densities  of  A)  saturate. 

The  consequences  of  flow  rates  of  NO  and  I2  on  gain  are 
shown  for  80  W  in  figure  1 1  (d).  For  a  fixed  power,  increasing 
the  mole  fraction  of  NO  at  the  inlet  reduces  the  density  of  O 
atoms  which  reduces  the  quenching  of  P.  Increasing  the  NO 
mole  fraction  increases  the  range  of  I2  flow  rates  over  which 
the  positive  gain  can  be  achieved  and  increases  the  flow  rate 
of  I2  at  which  the  maximum  gain  is  obtained.  For  example, 
the  maximum  gain  in  the  absence  of  NO  is  4.2  x  10-5  cm"1 
for  an  I2  flow  rate  of  0.8  seem.  As  the  rate  of  quenching  of 
P  decreases  with  increasing  NO  flow,  gain  increases  to  more 
than  1 .2  x  10-4  cm-1  for  an  inlet  mole  fraction  of  5%  NO  and 
flow  rate  of  3  seem  of  I2. 

The  consequences  of  flow  dynamics  on  the  densities  of 
02(‘A)  and  I*  are  shown  in  figure  12  with  results  from 
nonPDPSIM.  The  densities  of  C^A)  and  P  are  shown  for 
80  W  and  a  6slpm  inlet  flow  of  He/02/NO  =  68/30/2.  A 
100  seem  flow  of  He/l2  was  injected  from  the  second  nozzle 
with  the  I2  flow  varied  from  0.5  to  3  seem.  As  the  flow 
passes  through  the  discharge  zone,  electron  impact  produces 
02(!  A)  at  large  radius  first  where  the  plasmadensity  is  highest. 
Diffusion  homogenizes  the  density  of  2.5  x  1015  cm-3  across 
the  radius  by  10- 12  cm  downstream  of  the  plasma  zone.  At  and 
after  the  I2  injection  point,  the  densities  of  O2  ( 1 A )  decrease  due 
to  excitation  transfer  to  I2  and  pumping  of  P,  first  at  the  outer 
radius  where  the  12  is  injected  and  on  axis  4-5  cm  downstream. 
With  low  flow  rates  of  I2  (^0.5  seem),  the  (^('A)  is  not 
significantly  depleted  whereas  with  flow  rates  ^  1  seem,  the 
flow  of02(,A)  is  largely  consumed  by  reactions  with  I2  and  I. 

The  densities  of  P  are  maximum  adjacent  to  the  second 
nozzle  where  densities  of  both  C^A)  and  I2  are  largest. 
The  P  diffuses  radially  to  the  centre  of  the  discharge  4-5  cm 
downstream  of  the  injection  point.  For  small  flow  rates  of  I2, 
the  densities  of  1  are  smaller,  and  hence  consumption  of  O2  ( 1 A ) 
in  pumping  the  P  is  gradual,  leading  to  a  longer  distance  over 
which  significant  amounts  of  I*  are  present.  At  high  flow  rates 
of  I2,  CM1  A)  is  largely  consumed  in  the  vicinity  of  the  injection 
point  and  so  I*  does  not  extend  appreciably  beyond  that  point. 

The  densities  of  I  and  P,  and  the  gain  obtained  with 
nonPDPSIM  at  a  radius  of  1.5  cm  are  shown  in  figure  13  for 
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Figure  12.  Etcited  state  densities  while  varying  the  I2  flow  rate 
through  the  second  nozzle,  (a)  02('  A)  and  ( b )  P.  These  results, 
obtained  with  nonPDPSIM ,  are  for  3Torr,  80  W  power  deposition, 
and  6slpm  flow  of  He/02/NO  =  68/30/2. 

the  conditions  of  figure  1 2.  Similar  to  the  results  obtained 
with  GlobaiKIN ,  the  densities  of  I  increase  with  increasing 
flow  rate  of  I2  indicating  that  the  system  is  in  an  I2  limited 
regime.  The  densities  of  F  also  increase  with  increasing  I2 
mole  fraction  reaching  a  maximum  3-4  cm  downstream  of 
the  injection  point.  P  is  rapidly  quenched  by  the  O  atoms 
whose  dens  ties  are  commensurate  with  CM*  A).  At  lower 
I2  mole  fractions,  the  densities  of  P  do  not  decrease  rapidly 
downstream  of  the  injection  point  due  to  the  availability  of 
O2 ( *  A )  to  continue  to  pump  1  to  P. 

The  maximum  gain  (at  a  radius  of  1 .5  cm)  is  1 0"4  cm" 1  for 
an  I2  flow  of  3  seem.  The  axial  extent  of  positive  gain  is  limited 
by  the  depletion  of  C^C1  A).  Lower  flow  rates  of  I2  produce 
lower  gain  b  Jt  the  axial  extent  of  gain  is  greater.  A  comparison 
of  maximum  gain  obtained  with  GlobaiKIN  and  nonPDPSIM 
as  a  function  of  flow  rate  of  I2  is  shown  in  figure  1 3(d).  The 
predicted  gains  are  commensurate  except  at  large  flow  rates  of 
I2  where  the  gain  with  GlobaiKIN  is  significantly  higher  than 
with  the  2d  model.  The  lack  of  axial  transport  in  the  plug  flow 
model  produces  artificially  high  rates  of  reaction  between  O 
and  02(*  A)  with  1  and  I2,  which  increases  predicted  gain. 
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Figure  13.  Iodine  species  densities  and  gain  at  a  radius  of  1.5  cm  as 
a  function  of  the  flow  rate  of  I2  through  the  second  nozzle.  Densities 
of  (a)  I  and  (b)  I*,  (c)  gain  and  (d)  comparison  of  the  plug  flow  and 
2D  models  for  gain  2  cm  downstream  of  the  second  nozzle.  These 
results  are  from  nonPDPSIM  for  3  Tou,  80  W  power  deposition  and 
6slpm  flow  of  He/02/NO  =  68/30/2.  1 00  seem  of  He/12  is  flowed 
through  the  second  nozzle. 


5  NO2  injection  downstream  of  the  plasma  zone  and  prior  to  the  addition  of 

1 2  [9,34].  Comparisons  between  predictions  from  GlobaiKIN 
In  experimental  demonstrations  of  laser  oscillation  and  and  experiments  by  Carroll  et  al  [34]  using  this  strategy  are 
significant  gain  in  eCOIL  systems,  NO2  was  usually  injected  shown  in  figure  14.  The  conditions  are  a  lOTorr,  26.9  slpm 
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Figure  14.  Gain  in  a  subsonic  flow  for  10  Torr,  26.9 slpm  of 
He/02  =  80/20,  25-800  W,  with  10.8  seem  of  I2  injected  through 
the  second  nozzle,  (a)  Osccm  N02,  (b)  672  seem  N02  and  (c) 

1344  seem  N()2  through  the  first  nozzle.  Experimental  values  are 
from  [34], 

(20 mmols-1)  inlet  flow  of  He/02  =  80/20  and  25-800 W 
followed  by  injection  of  0-1 344 seem  (0-1  mmols-1)  N02 
and  injection  of  10.7  seem  (0.008  mmol  s-1)  of  I2,  equivalent 
to  few  per  cent  of  the  02(,A)  flow  rate.  The  diameter  of 
the  reactor  s  4.9  cm  and  I2  injection  is  20  cm  downstream 
of  the  N02  injection  point.  The  power  deposition  spans 
nearly  25  cm  due  to  a  larger  separation  between  electrodes  in 
the  experiment.  The  experimental  measurements  were  made 
10  cm  downstream  of  the  I2  injection  point  in  a  subsonic  (high 
gas  temperature  flow)  and  so  gains  are  negative. 

With  the  exception  of  low  powers  and  low  flow  rates 
of  N02,  the  experimental  trends  are  captured  by  GlobalKIN. 
Addition  of  N02  prior  to  injection  of  I2  scavenges  some  of 
the  O  atoms  in  the  flow  and  so  reduces  the  quenching  of  1* 
by  O  atoms.  Higher  powers  produce  larger  flows  of  O  atoms 
as  well  as  more  02(l  A)  but  the  quenching  of  F  dominates. 
Increasing  the  N02  flow  rate  increases  the  scavenging  of  O 
atoms  and  extends  the  power  prior  to  transitioning  to  large 
negative  gain. 


Figure  15.  Consequences  of  injection  of  N02  through  the  first 
nozzle  for  3  Torr,  40  W,  He/02/NO  =  67/30/3  and  6  slpm.  36  seem 
of  He/N02  mixture  is  injected  through  the  first  nozzle.  Densities  of 
(a)  O  and  (b)  02(*  A)  and  (c)  yield  of  02(‘  A).  Addition  of  N02 
rapidly  consumes  the  O  atoms.  These  results  are  from  the 
GlobalKIN ,  the  plug  flow  model.  Values  from  the  2D  nonPDPSIM 
are  shown  without  N02  injection. 

Since  N02  is  more  effective  than  NO  in  scavenging  of 
O  atoms,  we  investigated  N02  injection  through  the  first 
nozzle.  The  conditions  are  a  3  Torr,  6  slpm  inlet  flow  of 
He/02/NO  =  67/30/3  and  discharge  power  of  40  W.  A 
He/N02  flow  of  36  seem  was  injected  through  the  first  nozzle 
with  the  fraction  of  N02  being  varied.  As  before  a  100  seem 
flow  of  He/12  =  99/1  was  injected  through  the  second  nozzle. 

The  consequences  of  N02  flow  rate  through  the  first 
nozzle  on  the  densities  of  O  and  02(*  A),  and  yield  of  02(‘  A) 
are  shown  in  figure  15.  These  results  are  from  GlobalKIN 
with  a  result  from  nonPDPSIM  without  N02  injection  for 
comparison.  The  injection  of  N02  produces  a  decrease  in 
the  O  atom  density  due  to  the  titration  of  O  by  N02  and  the 
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Figure  16.  Gain  when  including  NO  in  the  inlet  flow  and  injecting 
NO2  through  he  first  nozzle  for  3Torr,  40  W, 

He/02/N0  =  70  -  jc/30/jc  and  6slpm.  100 seem  of  He/I2  =  99/1 
is  injected  through  the  second  nozzle,  (a)  Gain  as  a  function  of  NO 
mole  fraction  and  (b)  gain  as  a  function  of  N02  flow  rate.  These 
results  are  fro ti  the  plug  flow  model.  Maximum  gain  is  obtained  at 
low  NO  flow  rates  and  high  N02  injection. 

conversion  of  O  to  O2,  respectively.  As  the  NO2  flow  rate 
increases  to  36  seem  the  O  atom  densities  decrease  by  a  factor 
of  3-4  just  downstream  of  the  injection  point  due  to  the  more 
rapid  rate  of  reaction  with  O  (compared  with  that  of  NO).  The 
heat  of  reaction  between  NO2  and  O  locally  increases  the  gas 
temperature  leading  to  a  reduction  in  the  density  of  (^('A) 
near  the  first  nozzle  due  to  rarefaction.  Note,  however,  that  the 
yield  of  CM1  A)  is  not  affected  by  NO2  injection  because  NO2 
does  not  appreciably  quench  (^(’A).  The  reduction  in  the 
densities  of  O  with  NO2  injection  implies  that  the  quenching 
of  I*  by  O  is  'educed.  This  produces  a  reduction  in  the  amount 
of  CM1  A)  used  in  pumping  the  I*.  Hence,  downstream  of  the 
second  nozzle,  the  yield  of  02(*  A)  is  higher  for  larger  NO2 
mole  fractions.  Results  from  GlobalKIN  and  nonPDPSIM  are 
in  general  agreement  except  downstream  of  the  12  injection 
point  due  to  the  artificially  higher  rates  of  reaction  upon 
injection  of  I2. 

Gain  is  shown  in  figure  1 6  for  2  cm  downstream  from  the 
second  nozz  e  as  a  function  of  NO  mole  fraction  in  the  inlet 
flow  and  flow-  rate  of  NO2  through  the  second  nozzle.  The  flow 
conditions  aie  3  Torr  and  40  W  power  deposition.  Having  NO 
in  the  inlet  flow  affects  the  production  of  CM1  A)  as  discussed 
above  as  well  as  managing  the  O  atoms  density.  Injection 
of  NO2  dow  istream  of  the  discharge  largely  only  affects  the 
density  of  O  atoms  (and  gas  temperature).  As  such,  at  low 
values  of  NO  flow,  injection  of  NO2  is  effective  in  managing 
the  O  atom  density  and  larger  flow  rates  tend  to  maximize  gain 
by  reducing  quenching  of  P  by  O  atoms.  At  large  flow  rates 
of  NO,  the  management  of  O  atoms  is  dominated  by  reactions 
with  NO,  and  so  the  injection  of  NO2  is  less  effective.  Since 


there  is  a  deleterious  effect  on  C^1  A)  production  by  having 
large  flows  of  NO  through  the  discharge,  managing  the  O  atom 
density  with  injection  of  N02  is  likely  the  optimum  strategy. 

6.  Concluding  remarks 

The  consequences  of  NO  in  the  inlet  flow  of  a  He/02  plasma 
and  its  flowing  afterglow,  and  NO2  and  I2  injection  on  the 
post-discharge  kinetics  of  the  eCOIL  were  investigated  using 
plug  flow  and  2D  models.  The  addition  of  NO  to  the  inlet 
flow  through  the  discharge  produces  a  reduction  in  Tc  and  a 
modest  increase  in  nc  resulting  in  the  densities  and  yields  of 
02(l  A)  being  generally  lower  with  NO.  Including  NO  in  the 
flow  reduces  the  density  of  O  atoms  both  by  a  reduction  in  the 
electron  impact  dissociation  of  O2  and  by  exothermic  reactions 
of  O  with  NO.  This  proves  beneficial  to  improving  optical  gain 
by  reducing  the  quenching  of  I*  by  O  atoms.  At  higher  power 
deposition,  the  dissociation  of  O2  saturates  the  yield  of  A). 
By  virtue  of  adding  NO  to  the  inlet  flow,  the  reduction  in  Tc 
reduces  the  rate  of  dissociation  of  02.  Even  though  the  yields 
of  02(‘  A)  were  generally  lower,  the  optical  gain  was  generally 
higher  when  the  NO  mole  inlet  mole  fraction  was  between  1% 
and  3%. 

The  eCOIL  system  can  operate  in  power  limited  and 
I2  limited  regimes.  At  low  flow  rates,  I2  is  nearly  totally 
dissociated  and  so  the  densities  of  I  depend  largely  on  the 
I2  flow  rate.  Upon  increasing  the  flow  rate  of  I2,  the  system 
transitions  to  a  power  limited  regime  and  higher  powers  are 
required  to  optimize  gain.  Small  flows  of  N02  in  the  post¬ 
discharge  region  can  be  used  to  fine  tune  the  gain.  The  addition 
of  NO2  rapidly  consumes  O  atoms  without  significantly 
changing  other  parameters  (other  than  Tg)  and  so  increases 
the  optical  gain.  The  injection  of  NO2  was  most  effective  at 
low  flow  rates  of  NO.  In  general,  management  of  the  O  atom 
density  is  critical  to  optimizing  gain  due  to  its  rapid  rate  of 
quenching  of  I*. 
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Chemical  oxygen-iodine  lasers  (COILs)  oscillate  on  the  ^  Ml  *  ^3/2  transition  of  atomic  iodine  at 
1.315  /zm  by  a  series  of  excitation  transfers  from  02(,A).  In  electrically  excited  COILs  (eCOILs), 
the  02(,A)  is  produced  in  a  flowing  plasma,  typically  He/02,  at  a  few  to  tens  of  Ton*.  Many  system 
issues  motivate  operating  eCOILs  at  higher  pressures  to  obtain  larger  absolute  densities  of  02(’A) 
for  a  given  yield  and  to  provide  higher  back  pressure  for  expansion.  In  this  paper,  we  discuss  results 
from  a  computational  investigation  of  02('A)  production  in  flowing  plasmas  sustained  at  moderate 
pressures  (^50  Torr).  Power  deposition  and  flow  rates  were  scaled  such  that  in  the  absence  of 
second  order  effects,  yield  should  be  constant  and  absolute  02(!A)  production  should  scale  linearly 
with  pressure.  We  found  in  many  cases  that  absolute  02( 1  A)  production  scaled  sublincarly  with 
pressure.  Ozone  is  found  to  be  one  of  the  major  quenchers  of  02('A)  and  its  production  increases 
with  pressure.  Gas  heating  also  increases  with  increasing  pressure  due  to  exothermic  three-body 
reactions.  The  gas  heating  reduces  03  production,  increases  03  destruction  and,  for  certain 
condition!;,  restores  yields.  With  increasing  pressure  and  increasing  absolute  densities  of  atomic 
oxygen  and  pooling  reactions  of  02(IA),  quenching  by  these  species  also  becomes  important, 
though  th2  influence  of  O-atom  quenching  can  be  controlled  by  managing  the  density  of  O  atoms 
with  additives.  The  yield  of  02(‘A)  is  also  determined  by  discharge  stability  which  becomes 
problemalic  at  higher  pressure.  ©  2007  American  Institute  of  Physics.  [DOI:  10.1063/1.2743878] 


I.  INTRODUCTION! 

Chemical  oxygen-iodine  lasers  arc  being  investigated 
because  of  their  ability  to  be  delivered  through  optical  fibers 
(1.315  /an),  their  highly  scalable  continuous  wave  power, 
and  favorable  material  interaction  properties.1-^  Operation  of 
the  chemical  oxygen-iodine  laser  (COIL)  is  based  on  an  elec¬ 
tronic  transition  be  tween  the  spin-orbit  levels  of  the  ground 
state  configuration  of  the  iodine  atom,  \(2P]/2)  — • 

The  upper  level  is  sopulated  by  near  resonant  energy  transfer 
(quantum  defect  '-219  cm"1)  from  Oo(‘A)  to  ground  state 
K  2Pvi)-  The  COIL,  generates  ground  state  I  atoms  by  disso¬ 
ciative  excitation  transfer  from  02(!A)  to  I2.  Typically, 
02(!A)  is  produced  in  an  external  chemical  reactor  by  a  gas- 
liquid  reaction  between  gaseous  chlorine  and  a  basic  hydro¬ 
gen  peroxide  solution  producing  yields  approaching  10()% 
of  the  oxygen  emerging  in  the  02(!A)  state.  The  long  life¬ 
time  of  02(!A)  (spontaneous  lifetime  of  about  60  min)  and 
robustness  against  quenching  enables  transport  over  long  dis¬ 
tances  to  the  laser  cavity. 

In  electrically  excited  COILs  (eCOILs),  the  02(lA)  is 
produced  in  a  flowing  plasma,  typically  He/02,  at  a  few  to 
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tens  of  Torr.  Helium  is  often  used  as  a  buffer  gas  to  lower  the 
discharge  gas  temperature  and  to  provide  a  more  optimum 
electron  temperature  for  excitation  of  02(!A).  Unlike  the 
conventional  COIL,  the  electric  discharge  variant  also  pro¬ 
duces  atomic  oxygen  by  electron-impact  dissociation  of  mo¬ 
lecular  oxygen.  Although  the  atomic  oxygen  is  a  quencher  of 
the  upper  laser  level,  it  can  also  aid  in  the  dissociation  of  I2. 
Recent  and  ongoing  investigations  have  shown  that  substan¬ 
tial  yields  of  02(!A)  can  be  generated  by  an  appropriately 

“  |  -i 

tailored  electric  discharge'  in  mixtures  of  02  with  rare  gas 
diluents.  Scaling  studies  of  non-self-sustained  discharges 
where  ionization  is  supported  by  an  external  source  suggest 
that  CM1  A)  yields  approaching  25%-3()%  might  be  possible. 
This  is  accomplished  by  tailoring  the  electron  temperature  to 
more  closely  match  the  energy  range  where  the  excitation 
cross  sections  for  02(*A)  are  largest.1""14 

The  three  main  components  of  COIL  devices  are  the 
02(!A)  generator,  a  supersonic  nozzle,  and  the  laser  cavity. 
The  primary  purpose  of  the  supersonic  nozzle  is  to  lower  the 
temperature  in  the  laser  cavity,  thereby  reducing  the  likeli¬ 
hood  of  the  energy  transfer  process  from  02(!A)  to  I(2P|/2) 
proceeding  in  the  reverse  direction.  Recently,  positive  gain 
and  laser  oscillation  in  atomic  iodine  in  a  supersonic  flow 
optical  cavity  have  been  reported  resulting  from  electric  dis¬ 
charge  produced  02(,A).IS"21  An  overview  of  02('A)  gen¬ 
eration  in  plasmas  is  given  in  Ref.  22. 

©  2007  American  Institute  of  Physics 
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Many  system  issues  motivate  operating  eCOILs  at 
higher  pressures  to  obtain  larger  densities  of  02(1A)  for  a 
given  yield  and  to  provide  higher  back  pressure  for  expan¬ 
sion.  If  higher  order  effects  are  not  important,  and  param¬ 
eters  such  as  energy  deposition  per  molecule  are  maintained 
constant,  it  is  expected  that  absolute  02(,A)  production  [that 
is,  the  total  number  of  02(*A)  molecules  produced]  should 
scale  linearly  win  pressure,  thereby  providing  additional 
motivation  for  higher  pressure  operation. 

In  previous  work,  the  scaling  of  production  of  02(!A) 
was  computationally  investigated  using  global-kinetics,  one¬ 
dimensional,  and  wo-dimensional  (2D)  models.13 14  It  was 
found  that  the  yield  of  02(,A)  generally  scaled  linearly  with 
energy  deposition  at  low  pressures  (<10Torr)  until  the 
ground  state  02  is  depleted.  In  this  work,  we  report  on  results 
from  a  computational  investigation  of  radio  frequency  (rf) 
excited  flowing  He/02  plasmas  using  plug-flow  and  two- 
dimensional  (2D)  models.  The  emphasis  is  on  developing 
scaling  laws  for  CM1  A)  production  in  eCOIL  systems  when 
operating  at  higher  pressures.  We  found  that  although  yields 
may  decrease  with  increasing  pressure,  absolute  densities  of 
CM1  A)  typically  do  increase.  Ground  state  and  vibrationally 
excited  ozone  arc  major  quenchers  of  02(IA),  and  produc¬ 
tion  of  these  species  increase  with  pressure.  With  increasing 
absolute  densities  of  atomic  oxygen  and  pooling  reactions  of 
02(’A).  quenching  by  these  species  also  become  important, 
though  the  influence  of  O-atom  quenching  can  be  minimized 
by  managing  the  O-atom  density  with  additives.  The  yield  of 
02(1A)  is  also  determined  by  discharge  stability  which  be¬ 
comes  problematic  at  higher  pressures. 

The  models  and  reaction  mechanism  are  briefly  de¬ 
scribed  in  Sec.  II.  The  results  from  our  investigations  using 
the  plug-flow  model  are  discussed  in  Sec.  Ill,  followed  by  a 
discussion  of  results  obtained  from  the  2D  model  in  Sec.  IV. 
Concluding  remarks  are  in  Sec.  V. 


II.  DESCRIPTION  OF  THE  MODELS 

This  investigalion  was  performed  with  plug-flow  and  2D 
models.  The  plug-3ow  model  GLOBAL_KIN  addresses  gas- 
phase  chemistry  ard  transport,  solution  of  Boltzmann's  equa¬ 
tion  for  the  electron  energy  distribution,  and  equations  for 
axial  fluxes  of  mass,  momentum,  gas  energy,  and  electron 
energy.  The  mode  is  described  in  detail  in  Ref.  13  and  so 
will  be  only  briefly  discussed  here. 

GLOBAL_KIN  is  a  volume  averaged,  global-kinetics 
model  for  plasma  chemistry.  When  operated  in  a  plug-flow 
mode,  axial  transport  is  also  approximated.  GLOBAL_KIN 
contains  a  plasma  chemistry  module,  a  surface  kinetics  mod¬ 
ule.  and  an  electrcn  energy  transport  module.  In  the  plasma 
chemistry  module,  the  time  rate  of  change  of  species  and 
temperatures  (eleciron  and  gas)  are  obtained  integrating  their 
respective  conservation  equations.  Assuming  a  linear  axial 
flow  in  a  cylindrical  tube,  transport  to  radial  surfaces  is  in¬ 
cluded  by  using  a  diffusion  length.  The  reaction  of  fluxes  on 
surfaces  is  addressed  by  the  surface  kinetics  module  which 
utilizes  a  surface  site-balance  model.  Electron-impact  rate 
coefficients  arc  piovided  by  the  electron  energy  transport 
module  where  the  electron  energy  distributions  are  obtained 


by  solving  Boltzmann’s  equation.  By  simultaneously  calcu¬ 
lating  the  axial  speed  of  the  flow  based  on  constant  pressure, 
change  in  enthalpy,  species  densities,  and  gas  temperature, 
the  integration  in  time  is  mapped  to  axial  position.  Although 
computationally  fast,  the  weakness  of  this  method  is  neglect¬ 
ing  axial  transport  by  diffusion. 

The  2D  model  used  in  this  study,  nonPDPSIM.  is  a  mul¬ 
tifluid  hydrodynamics  simulation,  described  in  detail  in  Ref. 
14.  in  which  transport  equations  for  all  charged  and  neutral 
species  and  Poisson’s  equation  are  integrated  as  a  function  of 
time.  Poisson’s  equation  [Eq.  (I)],  transport  equations  for 
conservation  of  the  charged  species  [Eq.  (2)],  and  the  surface 
charge  balance  equation  [Eq.  (3)]  are  simultaneously  inte¬ 
grated  using  a  Newton  iteration  technique, 

-  V  •  (e,)t'r  V  <t>)  =  2)  Nflj  +  Ps •  ( I ) 


ON, 

Pi 


=  - V-  Tj  +  S,. 


(2) 


^  =  2  <?,(-  Vl',  +  Sj )  -  V  •  <o[-  V4>)].  (3) 

m  J 

where  £0.  En  ^  Pv  A^,  l'y,  cr,  Sr  and  are  the  permittivity 
of  free  space,  dielectric  constant,  electric  potential,  surface 
charge  density,  species  density  and  flux,  conductivity  of  solid 
materials,  sources,  and  charge,  respectively.  The  subscripts  j 
denote  gas-phase  species.  Updates  of  the  charged  particle 
densities  and  electric  potential  are  followed  by  an  implicit 
update  of  the  electron  temperature  by  solving  the  electron 
energy  equation  for  average  energy,  e. 

= <?r,  •  e  -  #.,2  a -  v  ■  (|«r, -  x,  v  r,).  (4) 

where  jkTr=E.  The  terms  in  Eq.  (4)  are  for  the  contribution 
from  Joule  heating,  elastic  and  inelastic  impact  processes 
with  heavy  neutrals  and  ions  with  energy  loss  /c;.  and  elec¬ 
tron  heat  flux  consisting  of  terms  for  electron  energy  flux 
(l'r)  and  a  conduction  (X,  is  the  electron  thermal  conductiv¬ 
ity).  The  electron  transport  coefficients  and  rate  coefficients 
for  bulk  electrons  as  a  function  of  Te  are  obtained  by  solving 
the  zero-dimensional  Boltzmann’s  equation  for  the  electron 
energy  distribution  to  capture  the  non-Maxwellian  nature  of 
the  electron  swarm. 

The  fluid  averaged  advective  velocity  v  is  obtained  by 
solving  a  modified  form  of  the  compressible  Navier-Stokes 
equations  in  which  momentum  transfer  from  ion  and  eleciron 
collisions;  and  acceleration  by  the  electric  field  are  included 
in  the  momentum  equation,  and  Joule  heating  is  included  in 
the  energy  equations. 

—  =  -V  (pv)  +  /\  (5) 

rft 

—  ^  =  -  Vp -  V  •  (pvv)  -  V  f  +  2 

<>'  , 

(6) 


Downloaded  01  Jul  2007  to  129.186.96.224.  Redistribution  subject  to  AIP  license  or  copyright,  see  http://jap.aip.org/jap/copyright.jsp 


1 23306-3  Babaeva,  Arakoni,  and  Kushner 


J  Appl.  Phys.  101,  123306  (2007) 


,lip(i'T)  =  _  Vi- *•  V  T+  pvcpT)  +  2  j;  •  E  -  2  RAH, 

(ft  j  j 

+  /?  V  V,  (7) 

where  P  represents  the  inlet  and  outlet  flows  (the  inlet  flow  is 
specified  while  the  output  flow  is  adjusted  to  maintain  a  con¬ 
stant  mass  flux),  p  is  the  total  mass  density,  p  is  the  thermo¬ 
dynamic  pressure,  r  is  the  viscosity  tensor,  cp  is  the  heal 
capacity,  k  is  the  species  averaged  thermal  conductivity,  /x  is 
the  mobility,  and  M  the  molecular  weight.  AH,  is  the  change 
enthalpy  due  to  reaction  /'  having  total  rate  Rr  The  reactions 
include  Frank-Cordon  heating  from  electron-impact  disso¬ 
ciation  of  molecules  as  well  as  conventional  chemical  reac¬ 
tions.  The  sums  (other  than  for  reactions)  are  over  all 
charged  and  neutral  species.  The  contributions  to  momentum 
from  charged  particles  include  those  of  electrons.  The  con¬ 
tributions  to  the  energy  equation  from  Joule  heating  include 
contributions  from  ions.  The  heat  transfer  from  electrons  is 
included  as  a  collisional  change  in  enthalpy.  The  relationship 
between  pressure,  density,  and  temperature  is  given  by  the 
ideal-gas  law.  The  numerical  grid  in  nonPDPSIM  uses  an 
unstructured,  cylir drically  symmetric  mesh  with  triangular 
elements. 

The  reaction  mechanism  for  He/02  plasmas  used  here  is 
essentially  the  same  as  that  described  in  Refs.  1 3  and  14,  and 
involves  reactions  in  the  gas-phase  discharge  and  afterglow 
as  well  as  recombination  and  quenching  reactions  on  the 
discharge  tube  wails.  The  species  in  the  model  are  ground- 
state  neutrals  02,  O,  03,  and  He:  02(v)  (the  first  four  vibra¬ 
tional  levels  of  02),  electronic  states  CL^A),  02(iS),  O (’d), 
CK’S).  and  He(2S)  and  the  ions  02+,  0\  O,",  O',  03“,  and 
He*. 

The  precursors  to  the  eCOIL  laser,  02('A)  (0.98  eV)  and 
CM1!)  (1.63  eV),  are  dominantly  produced  in  the  discharge 
by  direct  electron  mpact  with  the  ground  state. 

e  +  02  — »  02( 1  A)  +  e,  (8) 

c  +  02-02(,I)  +  <’.  (9) 

and  in  reactions  of  electronic  excitation  transfer  from  meta¬ 
stable  0(lD)  atoms, 

0('d)  +  02-02('2)  +  0.  (10) 

Production  of  02('l)  generally  also  results  in  generation  of 
0?(1A)  through  rapid  collisional  quenching  reactions  with 
atomic  oxygen  and  ozone.  The  atomic  oxygen  is  dominantly 
produced  by  electron-impact  dissociation  of  02, 


e  +  02  — *  O  +  O  +  e. 

(ID 

e  +  02  -*0(lD)  +  0  +  <’. 

(12) 

The  02(!A)  persists  far  into  the  afterglow  due  to  its  long 
radiative  lifetime,  where  the  most  significant  quenching  re¬ 
actions  are  collisions  with  O  atoms  and  03.  and  energy  pool¬ 
ing  with  02(,A).  llie  rate  of  quenching  of  02(!A)  by  colli¬ 
sions  with  the  walls  is  uncertain  due  to  the  variability  of  the 
quenching  probability  with  temperature  and  conditions  of  the 
wall.  Based  on  estimates  for  similar  conditions,  we  have  as¬ 
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FIG.  1.  Schematic  of  the  discharge  tube,  (a)  Computational  domain.  The 
flow  enters  uniformly  from  the  left  where  pressure,  speed,  and  temperature 
are  specified.  The  walls  arc  held  at  constant  temperature  and  mass  flux  is 
conserved  at  the  outlet,  (b)  Close-up  showing  ring  electrodes. 


signed  a  wall  quenching  coefficient  of  10"  \  which  results  in 
insignificant  wall  quenching  for  typical  eCOIL  conditions. 

The  effective  yield  of  02(JA)  is  defined  as  the  ratio  of 
the  combined  02(,A)  and  02(‘2)  densities  to  the  sum  of  the 
densities  of  all  oxygen-containing  species  on  a  molecular  02 
basis, 

_ [O/AIMO^V)] _ 

([02]+[02(v)]+f02('A)]+[02(  ’2)]  +  0.5(0]  +  1.5(0,]) ' 

(13) 

This  choice  of  yield  was  made  with  the  prior  knowledge  that 
the  majority  of  02('S)  is  quenched  directly  to  02( 1  A) .  Typi¬ 
cally  in  the  afterglow  region,  the  density  of  02( 1  i)  is  negli¬ 
gibly  small.  As  such.  Eq.  (13)  is  most  relevant  for  best  case 
energy  scaling. 

In  our  model,  the  first  four  vibrational  states  of  02  arc 
lumped  into  a  single  effective  vibrational  state  02(v)  having 
excitation  energy  of  0.19  eV.  In  addition  to  supcrclastic  elec¬ 
tron  collisions  with  02(v),  all  electron  impact  processes  in¬ 
cluded  for  ground  state  O?  were  also  included  for  02(v)  with 
the  threshold  energy  shifted  by  0.19  eV.  In  a  similar  manner, 
in  addition  to  collisional  quenching  of  02(v)  which  contrib¬ 
utes  to  gas  heating,  all  heavy  particle  reactions  included  for 
ground  state  02  were  also  included  for  02(v).  Although  we 
acknowledge  that  rate  coefficients  for  these  processes  may 
depend  on  the  vibrational  state,  there  is  limited  kinetic  avail¬ 
able  to  account  for  those  effects  and  so  they  have  not  been 
included  here.  Since  the  density  of  02(v)  is  typically  only 
10"2-10"3  that  of  ground  state  02  we  expect  that  the  error 
resulting  from  those  omissions  is  small. 

A  schematic  of  the  idealized  eCOIL  device  addressed  in 
this  study  is  shown  in  Fig.  I.  A  Hc/02  =  70/30  mixture  is 
flowed  through  a  cylindrically  symmetric  quartz  tube  60  cm 
in  length  and  6  cm  in  diameter.  A  rf  electric  discharge  is 
operated  between  two  ring  electrodes  2  cm  wide  with  centers 
separated  by  13  cm.  The  electrodes  are  powered  up  to  a  few 
kW  at  25  MHz.  The  flow  rate  for  each  pressure  corresponds 
to  an  average  axial  inlet  speed  of  985  cm/s.  Thus  the  species 
residence  time  is  approximately  the  same  for  each  pressure 
provided  we  neglect  gas  heating.  Our  investigations  were 
limited  to  the  region  of  the  reactor  prior  to  supersonic  expan- 
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sion  and  injection  of  I2.  The  wall  temperature  was  held  fixed 
at  300  K  assuming  active  water-jacket  cooling. 

III.  SCALING  OI:  02C±)  PRODUCTION  WITH 
PRESSURE:  PLUG-FLOW  MODELING 

The  basic  scaling  laws  of  02(,A)  generation  as  a  func¬ 
tion  of  pressure  ir  flowing  He/02  plasmas  were  first  inves¬ 
tigated  using  the  plug-flow  model  GLOBAL_KIN  for  the  ge¬ 
ometry  shown  in  Fig.  1.  To  eliminate  as  many  systematic 
trends  as  possible  while  changing  pressure,  two  parameters 
were  kept  constant:  energy  deposition  per  inlet  oxygen  mol¬ 
ecule  a  (cV/mole:ule)  and  flow  residence  time  in  the  dis¬ 
charge  tube,  r.  To  keep  r  constant,  the  flow  rate  is  increased 
in  proportion  to  pressure.  To  keep  a  constant,  power  is  in¬ 
creased  in  proportion  to  pressure  provided  r  is  constant. 
Computationally,  this  is  accomplished  by  computing  the 
power  deposition  as  a  cycle  averaged  product  of  voltage 
times  current  and  comparing  that  power  to  the  value  required 
to  keep  a  constart.  Every  few  rf  cycles  the  voltage  is  ad¬ 
justed  (increased  or  decreased  by  a  maximum  of  a  few  per¬ 
cent)  in  the  direction  to  deliver  the  desired  power.  This  pro¬ 
cess  is  continued  until  a  quasi-steady-state  is  reached. 

Keeping  these  parameters  constant,  if  there  are  no  sec¬ 
ond  order  effects,  yield  should  remain  constant  and  absolute 
CM1  A)  production  should  scale  linearly  with  pressure.  We 
also  expect  that  the  electron  temperature  T,  and  gas  tempera¬ 
ture  should  weakly  depend  upon  pressure.  The  former 
insensitivity  results  from  electron  losses  being  dominated  by 
volumetric  processes  (i.e.,  attachment  and  dissociative  re¬ 
combination)  and  so  the  rate  of  loss  by  diffusion  does  not 
significantly  affect  the  electron  accounting.  The  latter  insen¬ 
sitivity  results  from  thermal  conduction,  which  is  the  domi¬ 
nant  temperature  legulating  loss  process,  not  being  a  sensi¬ 
tive  function  of  pressure  in  the  range  we  investigated. 

At  least  two  second  order  effects  will  be  seen  to  be  im¬ 
portant:  three-body  reactions  and  discharge  stability.  The 
three-body  reactions  which  produce  quenchers  of  (^(’A), 
such  as  Ov  are  detrimental  due  to  the  resulting  reduction  in 
02('A)  density.  Three-body  reactions  of  the  sort  A  +  #  +  A/ 
—+AB  +  M  are  typically  exothermic  and  can  be  a  significant 
heating  source  as  these  reactions  begin  to  dominate  at  higher 
pressures.  The  resulting  increase  in  gas  temperature  usually 
has  a  negative  effect.  The  onset  of  discharge  instabilities  at 
higher  pressures  is  potentially  harmful  due  to  the  nonuniform 
power  deposition  that  results. 

Plasma  parameters  obtained  with  the  plug-flow  model 
for  pressures  from  3  to  50  Torr  for  Hc/02  =  70/30  and  a 
=0.3  cV/molecule  (40-670  W)  are  shown  in  Figs.  2-5.  The 
power  profile  is  sj>ecified  based  on  the  results  from  the  2D 
model.  When  keeping  a  constant,  the  electron  density 
weakly  depends  upon  pressure,  increasing  in  peak  value 
from  1  X  1010  to  1  9X  I0ll)  cm"3  from  3  to  50  Torr.  There  is 
a  small  decrease  in  T(.  (from  2.7  to  2.4  cV)  when  increasing 
pressure  due  to  a  reduction  in  diffusion  losses  of  electrons. 
The  decrease  in  7r  then  reduces  the  power  dissipation  per 
electron  which  necessitates  an  increase  in  electron  density. 

Electron-impact  dissociative  excitation  of  02  results  in 
approximately  10%  dissociation,  producing  peak  O-atom 


FFG.  2.  Plasma  and  flow  parameters  for  3-50  Ton  (0=0.3  eV/molecule. 
He/02=70/30,  r=3  cm),  (a)  Power,  (b)  electron  density,  and  (c)  atomic 
oxygen  density.  In  the  discharge  region  up  to  10%  of  the  molecular  oxygen 
is  dissociated. 

densities  of  2.2  X  1015  to  1.6X  1016  cm"3  from  3  to  50  Torr. 
This  increase  in  O-atom  density  is  sublincar  with  respect  to 
pressure  due  to  the  increased  efficiency  of  consumption  of  O 
atoms  by  the  three-body  collisions  that  form  ozone. 

0  +  02  +  A7  —  03  +  M.  (14) 

This  reaction  has  a  rate  coefficient  of  6 X  I0_34(7',,/ 
300)“2  8  cm6  s-1  which  decreases  the  formation  of  ozone 
with  increasing  gas  temperature.  The  main  reaction  for  ozone 
destruction, 

O  +  O3  — *  02  +  02 ,  (15) 

has  a  rate  coefficient,  8X  10-12 exp(-2060/7’;lf)  cm3s_l.  The 
increase  in  exothermic  three-body  reactions,  such  as  Eq. 
(14),  with  increasing  pressure  results  in  significant  gas  heat¬ 
ing.  In  spite  of  a  being  a  constant,  the  result  is  that  the  peak 
gas  temperature  increases  from  330  to  520  K  when  increas¬ 
ing  pressure  from  3  to  50  Torr.  Note  that  at  3  Torr,  the  peak 
in  7^  occurs  in  the  discharge  region  where  plasma  heating 
(mostly  ion  Joule  heating)  is  the  dominant  heating  mecha¬ 
nism.  At  50  Torr,  T 'x  peaks  downstream,  a  consequence  of 
heating  from  exothermic  three-body  reactions  which  occur 
dominantly  outside  the  plasma  zone. 

The  combination  of  increasing  production  of  O3  due  to 
three-body  processes  and  increasing  destruction  of  O3  due  to 
an  increase  in  gas  temperature  results  in  the  axial  depen¬ 
dence  of  O3,  as  shown  in  Fig.  3(a).  At  low  pressure,  the  peak 
in  03  density  occurs  downstream.  The  low  rate  of  three-body 
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Discharge 


FIG.  3.  Density  of  oxygen  species  for  3-50  Torr  (nr=0.3  eV/molecule. 
He/0>= 70/30.  r=3  cm),  (a)  Ov  (b)  O/AJ+O^'S)  and  (c)  combined 
02('A)  +  02(  *2)  yield.  The  peak  of  ozone  production  is  in  the  upstream  side 
of  the  discharge  region  where  the  gas  temperature  is  lower. 

reactions  and  low  gas  temperature  (low  rates  of  03  destruc¬ 
tion)  result  in  mo>t  of  the  03  being  produced  only  after  a 
finite  flow  time  <nd  this  density  accumulates  in  the  dis¬ 
charge.  (The  average  residence  time  in  the  flow  tube  is 
60  ms.)  At  higher  pressures,  the  rate  of  three-body  reactions 


FIG.  4  Density  of  oxsgen  species  (or=0.3  eV/molecule,  He/0>  =  70/30.  r 
=  3  cm)  as  a  function  of  pressure  at  the  exit  of  the  tube,  (a)  Oil1  A) 
+  0_.('X)  and  (b)  combined  02('A)+CM '2)  yield. 


FIG.  5.  Gas  temperature  and  yield  for  3-50  Ton  (<*=0.3  eV/molecule. 
He/O2  =  7()/30.  r=3  cm),  (a)  Enthalpy  from  three-body  recombination  heal¬ 
ing  is  included,  (b)  recombination  heating  is  excluded,  and  Ic)  combined 
02('A)  +  0?( ’X)  density  for  30  and  50  Ton*  with  and  without  recombination 
heating.  For  these  conditions  there  is  a  decrease  in  yield  when  the  gas 
becomes  colder  due  to  the  increased  production  of  Ov 

is  proportionally  higher  resulting  in  large  rates  of  formation 
of  03  in  the  discharge  prior  to  there  being  significant  gas 
heating.  For  example,  the  peak  03  density  is  2.2 
X  1015  cm-3  at  50  Torr.  The  accompanying  increase  in  gas 
temperature  downstream  results  in  increased  rates  of  destruc¬ 
tion  of  03  by  collisions  with  O  atoms  which  ultimately  re¬ 
duces  its  density.  The  end  result  is  that  the  downstream  den¬ 
sity  of  03  is  maximum  at  (I  -2)X  1014  cm"3  at  pressures  of 
10-20  Torr. 

The  03  density  is  an  important  consideration  due  to  its 
being  an  efficient  quencher  of  CM1  A), 

02( 1  A)  +  03  —  02  +  02  +  O .  (16) 

With  a  rate  coefficient  of  5.2 X  I0_l  1  exp(-2840/Tv) 
cm3  s"1,  the  rate  of  quenching  increases  with  the  increase  in 
gas  temperature  occurring  at  higher  pressures.  [Note  that  the 
rate  coefficient  for  quenching  of  02('A)  by  vibrationally  ex¬ 
cited  03  is  more  than  an  order  of  magnitude  larger  than  for 
the  ground  state."3]  The  increase  in  03  and  the  increase  in 
gas  temperature  both  increase  quenching  rates.  This  results 
in  the  density  of  0/  [denoting  the  sum  of  the  densities  of 
02(!A)  and  02(,S)]  increasing  sublincarly  with  increasing 
pressure  and  in  the  decrease  of  total  yield  as  a  function  of 
pressure,  as  shown  in  Fig.  4.  At  low  pressure 
(^10  Torr)  with  constant  a,  the  yield  is  nearly  constant  and 
so  the  absolute  production  of  02*  scales  with  pressure.  At 
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higher  pressures,  tie  production  and  yield  of  02*  peak  in  the 
plasma  zone.  The  density  of  02  decreases  thereafter,  a  con¬ 
sequence  of  downstream  quenching  at  the  higher  gas  tem¬ 
peratures  in  the  presence  of  03.  The  result  is  sublinear  scal¬ 
ing  with  pressure. 

At  higher  pressures  contributions  to  gas  heating  are 
dominated  by  three-body  recombination  healing  and  Frank- 
Condon  heating  due  to  dissociative  electronic  excitation.  To 
distinguish  between  these  two  contributions,  we  artificially 
excluded  contribu  ions  to  gas  heating  from  changes  in  en¬ 
thalpy  A/y  from  three-body  recombination  reactions.  The  re¬ 
sulting  gas  temperatures,  shown  in  Fig.  5(b),  are  lower  by 
50-70  K  at  higher  pressures  (30  and  50  Torr).  For  lower 
pressures  (3-10  Torr)  the  differences  in  gas  temperature  are 
small  and  heating  is  dominated  by  Frank-Condon  and  Joule 
heating  in  the  plasma  zone.  Decreasing  the  temperature  in 
this  manner  increases  the  03  density  which  translates  directly 
into  higher  quenching  and  lower  exit  densities  of  02  .  This 
results  in  lower  O,  yields,  as  shown  in  Fig.  5(c). 

The  hierarchy  of  the  main  quenchers  of  Oil1  A)  for  10, 
30,  and  50  Torr  is  shown  in  Fig.  6  for  the  base  case  condi¬ 
tions.  Quenching  limes  are  determined  by  l/(fc[Af]),  where 
[M]  is  density  of  the  quencher  and  k  is  the  rate  constant. 
More  rapid  and  important  quenching  reactions  are  denoted 
by  shorter  quenching  times.  By  the  exit  of  the  discharge  tube, 
O3  is  the  primary  quencher  for  pressures  >10  Torr.  For  pres¬ 
sures  higher  than  30  Torr,  quenching  times  by  03  become 
comparable  to  the  gas  residence  time  which  is  about  60  ms. 
Atomic  oxygen  is  the  second  most  important  quencher. 
Quenching  times  for  02  and  02(  ]A)  are  of  the  order  of  a  few 
seconds  and  not  important.  Note  that  the  O3  density,  espe¬ 
cially  at  high  pressures,  is  two  orders  of  magnitude  smaller 
than  atomic  oxygen  but  it  is  the  dominant  quencher. 

Another  second  order  effect  that  complicates  pressure 
scaling  when  changing  the  tube  size  is  heat  transfer  to  the 
walls.  For  example,  plasma  and  flow  parameters  for  tubes 
having  radii  of  3  and  1  cm  while  keeping  other  conditions 
equal  (a=0.9  eV/ molecule  and  60  ms  residence  time)  are 
shown  in  Figs.  7  aid  8.  In  the  larger  tube  and  with  the  larger 
diffusion  length  for  heat  transfer  to  the  walls,  gas  heating  at 
50  Torr  produces  temperatures  of  nearly  800  K.  As  a  conse¬ 
quence,  03  produced  in  large  quantities  upstream  in  the  dis¬ 
charge  is  rapidly  depleted  in  the  afterglow.  The  prior  domi¬ 
nance  of  O3  as  a  quencher  is  reduced  to  be  comparable  to 
that  of  atomic  oxygen.  Yield  for  02*  at  the  end  of  the  tube  for 
50  Torr  is  10%,  down  from  14%  at  3  Torr. 

Reducing  the  tube  radius  results  in  more  efficient  gas 
cooling  due  to  the*  shorter  radial  diffusion  length.  For  ex¬ 
ample.  the  gas  temperature  at  the  tube  exit  at  50  Torr  drops 
from  750  K  for  a  radius  of  3  cm  to  430  K  for  a  radius  of 
1  cm.  As  a  result,  there  is  increasing  production  of  03  down¬ 
stream  as  the  temperature  decreases,  ozone  again  becomes 
the  main  quencher  of  02('A),  and  yield  decreases  to  7%. 

Ozone  and  gas  temperature  management  are  clearly  im¬ 
portant  considerations  in  pressure  scaling.  For  example,  pres¬ 
sure  scaling  for  a  tube  with  a  radius  of  1  cm  was  repeated 
keeping  a=0.9  eV /molecule  while  increasing  the  flow  rate 
to  reduce  the  residence  time  by  a  factor  of  9  to  approxi¬ 
mately  7  ms.  The  results  are  shown  in  Fig.  9.  Since  a  is 


Discharge 


FIG.  6.  Quenching  times  for  02('A)  for  (a)  10,  (b)  30,  and  (c)  50  Torr  (rr 
=  0.3  eV/molecule,  Hc/0*=70/30.  r=3  cm).  Ozone  is  the  main  quencher 
in  the  postdischargc  downstream  region.  For  pressures  higher  than  30  Ton 
the  ozone  quenching  time  becomes  comparable  with  the  gas  residence  time. 

constant  the  maximum  gas  temperature  is  nearly  the  same  as 
the  prior  case.  However,  the  shorter  residence  lime  down¬ 
stream  reduces  the  amount  of  conductive  cooling  and  so  the 
exit  gas  temperature  increases  to  610  K  for  50  Torr.  The 
higher  gas  temperature  increases  the  rate  of  03  destruction 
and  decreases  the  rate  of  its  formation,  thereby  reducing  03 
densities.  The  lower  03  densities  and  shorter  residence  times 
reduce  the  likelihood  that  quenching  of  02(’A)  will  occur. 
As  a  result,  yield  is  restored  to  13.3%  at  50  Torr. 

Another  factor  that  is  potentially  important  to  pressure 
scaling  is  the  possibility  of  three-body  quenching  of  02('A) 
by  O  atoms.  Previous  studies  of  excited  state  kinetics  in  dis¬ 
charges  sustained  in  oxygen11'15* 16  indicated  lower  than  ex¬ 
pected  yields  of  O^'A)  at  higher  pressures  where  the 
O-atom  density  is  high.  The  authors"  proposed  that  oxygen 
atoms  in  a  three-body  collision  were  responsible  for  the 
quenching, 

02( 1  A)  +  O  +  02 — ►  02  +  02  +  O.  (17) 

Good  agreement  was  obtained  between  modeling  and  experi¬ 
ments  with  a  rate  coefficient  of  1  X  10"32  cm6  s”1.  This  pro¬ 
cess  is  in  addition  to  the  two-body  quenching  of  02(*A)  by  O 
atoms, 

02(,A)  +  0-02  +  0.  (18) 

having  a  rate  coefficient  of  1.1  X  lO"10  cm3  s“‘. 

The  densities  of  02*  and  yields  for  10  and  50  Torr  with 
and  without  three-body  O-atom  quenching  of  02(*A)  arc 
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c) 


PIG.  7.  Plasma  and  flow  parameters  for  3-50  Torr  (rr=0.9  eV/molecule, 
He/0?  =  70/30)  for  a  ndius  of  3  cm  and  a  residence  time  of  60  ms.  (a)  Gas 
temperature,  (b)  cknsity,  (c)  quenching  times  at  50  Torr,  and  (d)  com¬ 
bined  02( 1  A) O »( *51)  yield.  A  long  residence  time  results  in  higher  gas 
temperature  and  lower  ozone  density. 


FIG.  8.  Plasma  and  flow  parameters  for  3-50  Ton  (n-=0.9  eV/molecule. 
He/0:  =  70/30)  for  a  radius  of  1  cm  and  a  residence  time  of  60  ms.  (a)  Gas 
temperature,  (b)  Oy  density,  (c)  quenching  limes  at  50  Torr.  and  (d)  com¬ 
bined  02( 1  A) 02( ' — )  yield.  The  smaller  tube  radius  provides  for  higher 
rates  of  gas  cooling  which  restores  the  ozone  density  in  the  afterglow.  As 
such  the  yield  for  50  Torr  decreases  from  10.8%  to  7%  compared  to  a  3  cm 
tube. 


shown  in  Fig.  10.  Yields  and  02*  densities  decrease  by  ap¬ 
proximately  a  fac.or  of  2  at  50  Torr  when  including  three- 
body  quenching,  ^or  pressures  greater  than  20  Torr,  three- 
body  quenching  by  O  atoms  is  commensurate  with  that  by 
O},  as  shown  in  Fig.  10(d).  Resolution  of  the  value  and  tem¬ 
perature  dependence  of  three-body  quenching  of  02(!A)  by 
O  atoms  is  clearly  important  to  the  pressure  scaling  of 
cCOIL  devices.  Alternatively,  gas  mixtures  can  be  selected  in 
which  O  atoms  created  by  the  discharge  are  scavenged  and 
so  minimize  the  influence  of  O  atoms  on  quenching  of 
02(1A).  Oxygen  atoms  are  also  efficient  quenchers  of  the 
l(2PU2)  upper  laser  level,  and  so  minimizing  the  O-atom 
density  also  directly  aids  in  maintaining  the  inversion. 

Many  methods  have  been  proposed  to  scavenge  O  atoms 
and  so  improve  laser  performance,  such  as  coating  of  the 
tube  walls  with  z  HgO.15  A  practical  and  recently  imple¬ 
mented  method  is  to  use  small  amounts  of  NO  or  N02  either 
in  the  initial  gas  flow  or  as  additives  downstream. “  NO  and 


N02  both  rapidly  react  with  O  atoms  in  a  chain  reaction  that 
has  the  net  effect  of  converting  O  atoms  to  02.  These  reac¬ 
tions  arc 

NO  +  O  +  M  — >  N02  +  M .  (19) 

N02  +  0-*N0  +  02,  (20) 

with  rate  constants  of  I  X  1 0"3I( 7^/300)"'  6  cm6  s_l  and 
4.21  X  10~12  exp(-273/7^)  cm?s_I,  respectively.  By  control¬ 
ling  the  O-atom  density  in  this  manner,  quenching  by  O*  is 
also  regulated  as  O  atoms  are  the  precursor  to  O^  formation. 
NO  addition  has  the  added  benefit  of  removing  03  by 

NO  +  O,  — 02  +  N02.  (21) 

having  a  rate  coefficient  of  1 .4  X  I0"12  cxp(-l  3 10/Tv) 
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FIG.  9.  Plasma  and  flow  parameters  for  3-50  Torr  (rr=0.9  cV/molecule, 
He/02=70/30)  for  a  r  id  ins  of  I  cm  and  a  residence  lime  of  7  ms.  (a)  Gas 
temperature,  (b)  03  density,  (c)  quenching  times  at  50  Torr.  and  (d)  com¬ 
bined  CM1  Aj  +  C^f'l)  yield.  A  smaller  residence  time  results  in  lower  gas 
temperature  in  the  discharge  region,  but  prevents  rapid  cooling  in  the  after¬ 
glow.  Yield  is  restored  to  13.3%. 

The  effect  of  NO  addition  in  high  pressure  scaling  when 
including  three-body  quenching  by  O  atoms  was  investigated 
using  the  plug-flow  model.  Densities  of  oxygen  species  (O. 
O3,  and  O/  and  yield  are  shown  in  Fig.  1 1  as  a  function  of 
pressure  and  inlet  NO  mole  fraction.  These  quantities  are  at 
the  exit  of  the  flow  tube.  In  the  absence  of  NO,  O  and  O3 
densities  increase  with  pressure  and  saturate  around 
10-20  Torr.  With  these  large  densities  of  O  and  O3  as 
quenchers,  the  density  of  02*  reaches  a  maximum  between 
10  and  20  Torr  and  slowly  decreases  at  higher  pressures. 
This  quenching  results  in  a  monotonic  decrease  in  yield. 

With  addition  of  up  to  1%  of  NO,  the  exit  densities  of  O 
and  O3  at  pressures  above  20-30  Torr  are  essentially  elimi¬ 
nated.  The  reduction  in  these  densities  directly  translates  into 
less  quenching  and  higher  exit  densities  of  O^’A).  With  \% 
NO,  the  absolute  densities  of  O^'A)  increase  with  pressure 
and  yields  are  largely  restored.  Although  three-body  quench¬ 
ing  of  02('A)  is  potentially  an  important  process  at  high 


FIG.  10.  Density  of  oxygen  species  (o=0.3  cV/moIcculc.  Hc/02  =  70/30. 
r- 3  cni)  with  and  without  three-body  quenching  with  O  atoms  (a) 
CU'Al  +  O^'S)  for  10  and  50  Torr,  (b)  combined  Ojt'AHO/I)  yield  for 
10  and  50  Torr.  (c)  combined  02(  'AI  +  O^'X)  yield  as  a  function  of  pres¬ 
sure  at  the  exit,  and  (d)  quenching  limes.  Accounting  for  three-body  quench¬ 
ing  results  in  substantial  decrease  in  yield  (dashed  lines).  Symbols  0(2h)  and 

indicate  atomic  oxygen  in  two-body  and  three-body  reactions. 

pressures,  the  consequences  of  this  process  can  be  minimized 
by  managing  the  density  of  O  atoms  with  small  amounts  of 
additives,  such  as  NO  or  N02. 

IV.  CONSEQUENCES  OF  HYDRODYNAMICS  ON 
PRESSURE  SCALING 

The  results  obtained  using  the  plug-flow  model  provide 
insights  into  the  kinetic  processes  that  affect  O^A)  yields 
when  increasing  pressure.  In  this  section  additional  consid¬ 
erations  will  be  discussed  with  results  from  the  2D  model 
that  more  realistically  represent  the  electrical  circuitry,  elec¬ 
trode  losses,  uniformity,  and  flow  considerations. 

Another  second  order  effect  that  may  influence  02('A) 
production  is  the  discharge  constriction  commonly  observed 
at  higher  pressures  due  to  the  smaller  mean-free  path  of  clcc- 
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FIG  II.  Plasma  properties  at  the  exit  of  the  flow  tube  (« 
=  0.3  eV/molecule,  He/O:  =  7()/30.  r=  3  cm)  as  a  function  of  pressure  and 
NO  addition.  Densities  of  (a)  O,  (b)  03.  and  (c)  02(  'Aj  +  O^'S),  and  (d) 
combined  02( ,A)  +  02(,1)  yield.  Three-body  quenching  by  O  atoms  is  in¬ 
cluded.  Addition  of  small  mole  fractions  of  NO  regulate  the  O-alom  density 
and  restores  yield 


trons  and  gas  heating.  The  shorter  mean-free  paths  result  in 
electron-impact  excitation  occurring  dominantly  in  regions 
of  high  electric  field  (that  is,  electron  transport  is  local  as 
opposed  to  nonlocal).  Gas  heating  and  the  resulting  rarefac¬ 
tion  produce  nonuniformities  in  gas  density  which  even  with 
uniform  electric  fields  may  produce  nonuniform  values  of 
£7/V  (electric  held  divided  by  gas  number  density). 

For  example,  plasma  properties  and  the  spatial  distribu¬ 
tion  of  species  densities  are  compared  in  Fig.  12  for  a 
=  0.3  eV/molecule  for  moderate  (10  Torr.  133  W)  and  high 
(50  Torr,  670  W)  pressures.  These  results  are  averages  over 
the  25  MHz  cycle  in  the  quasi-steady-state.  (When  compar¬ 
ing  these  results  io  those  from  the  plug-flow  model,  recall 
that  the  values  ob.ained  from  the  plug-flow  model  are  vol¬ 
ume  averages.  The  peak  values  as  a  function  of  radius  in  the 
2D  model,  in  man)  cases  occurring  on  the  axis,  will  be  larger 
than  those  from  the  plug-flow  model.)  At  10  Torr  the  elec¬ 
tron  density  is  smoothly  distributed  between  the  electrodes 
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FIG.  12.  (Color)  Plasma  parameters  [electron  density,  electron  temperature, 
gas  temperature,  and  densities  of  O.  O*.  and  Od'A),  and  total  gas  mass 
density]  for  (a)  10  and  (b)  50  Torr  (He/02  =  70/30.  25  MHz  rt  excitation, 
«=0.3  eV/ molecule).  All  values  are  averaged  over  one  rf  period.  The  flow 
is  from  the  left.  The  scales  are  linear  with  zero  minimum  values  with  the 
exception  of  TK  unless  the  number  of  decades  is  indicated  for  log  plots.  The 
maximum  value  is  indicated  in  each  figure.  Discharge  constriction  at  the 
higher  pressure  occurs  near  the  downstream  electrode 


with  a  peak  value  of  1.6X  10,0citT\  The  cycle  averaged 
electron  temperature  is  1.9  eV  in  the  bulk  plasma,  peaking  at 
2.4  eV  near  the  electrodes.  The  dissociation  fraction  of  the 
inlet  02  is  4.5%,  producing  a  peak  O-atom  density  of  8.3 
X  1015  cm"3.  The  O  atoms  arc  progressively  converted  to  03 
as  they  flow  towards  the  exit,  producing  an  exit  density  of  03 
of  2.4  X  1014  cm'3.  The  increase  in  gas  temperature  is  about 
130  K  which  decreases  downstream  where  the  03  density  is 
highest.  Due  to  the  low  gas  temperature.  CM'A)  quenching 
(mostly  by  ozone)  is  negligible. 
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Al  50  Torr.  the  discharge  becomes  inhomogeneous  with 
the  highest  electron  density  of  5  X  I0I()  cm"3  near  the  down¬ 
stream  electrode.  The  peak  electron  temperature  drops  to 
1.85  eV.  Due  to  recombination  and  Frank-Condon  heating, 
the  peak  gas  temperature  increases  to  1700  K.  The  rarefac¬ 
tion  of  the  gas  near  the  downstream  electrode  allows  for 
higher  proportional  power  deposition  there.  Oxygen  mol¬ 
ecules  are  more  highly  dissociated  (10%  fractional  dissocia¬ 
tion)  producing  a  peak  O-atom  density  of  9.4  X  10,6cm_\ 
Due  to  the  high  temperature  near  the  downstream  electrode, 
03  produced  in  the  upstream  portion  of  the  discharge  is 
largely  destroyed  prior  to  leaving  the  discharge  zone.  The  03 
density  is  gradually  restored  downstream  near  the  tube  walls 
where  the  gas  is  coolest.  The  02(,A)  density  maximizes  in 
the  discharge  region  where  the  power  deposition  is  concen¬ 
trated  near  the  downstream  electrode. 

The  reduced  electric  field  E/N  (electric  field/gas  number 
density)  is  not  uniformly  distributed  over  the  discharge.  The 
maximum  values  of  E/N  occur  in  the  sheath  regions  near  the 
electrode.  However,  at  the  pressures  of  interest,  the  majority 
of  power  deposition  occurs  in  the  bulk  plasma  where  the 
electron  density  is  maximum.  In  these  volumes,  typical  val¬ 
ues  of  E/N  (averaged  over  a  rf  cycle)  are  10-25  Td  (1  Td 
=  1 0” 1 7  V  cm2)  wiih  there  being  a  tendency  for  the  E/N  to 
decrease  with  increasing  pressure  as  diffusion  losses  de¬ 
crease.  For  example,  for  €r=0.3  eV/molccule.  E/N  in  the 
bulk  plasma  decreased  from  25  Td  at  3  Torr  to  l3Td  at 
50  Torr.  (These  values  do  not  significantly  change  for  higher 
values  of  cr.)  Sirce  the  fractional  power  deposition  into 
02('A)  is  maximum  for  E/N^  10- 15  Td.  operating  at 
higher  pressures  is  kinetically  more  efficient  since  a  larger 
fraction  of  power  deposition  is  transferred  to  02(,A). 

Just  as  managing  the  O-atom  and  03  densities,  and  the 
gas  temperature  are  important  to  scaling  02(,A)  to  higher 
pressures,  so  is  managing  the  stability  of  the  discharge.  The 
transition  from  a  diffusive  and  homogeneous  discharge  to 
one  prone  to  instabilities  or  constricted  states  is  not  unusual 
when  operating  at  higher  pressures  in  electronegative 
gases. 12  ,5  |q  At  pressures  lower  than  those  associated  with 
streamer  propagation  (as  in  corona  discharges"4)  the  onset  of 
instabilities  or  constrictions  results  from  a  nonuniformity  that 
produces  local  regions  of  higher  power  deposition  or  elec¬ 
tron  heating  that  diffusion  is  not  able  to  disperse  across  a 
larger  volume.  That  is,  electron  transport  is  local  as  opposed 
to  nonlocal. 

In  the  geometry  investigated  here,  pressure  scaling  is 
ultimately  limited  by  an  instability  initiated  by  discharge 
constriction  which  occurs  between  40  and  50  Torr.  This 
trend  is  shown  in  Fig.  13  where  the  electron  density  and 
power  deposition  arc  plotted  for  pressures  from  3  to  50  Torr 
for  a=0.3  eV/ molecule.  The  corresponding  gas  densities  are 
shown  in  Fig.  14.  Power  deposition,  diffusively  distributed  at 
3  Torr,  gradually  constricts  about  the  downstream  electrode, 
producing  a  corresponding  constriction  in  the  electron  den¬ 
sity.  This  constriction  results,  in  part,  from  an  increasing  de¬ 
gree  of  rarefaction  in  the  total  gas  density  in  going  from 
upstream  to  downstream  as  the  gas  heats.  The  lower  gas 
density  near  the  downstream  electrode  enables  a  higher  con¬ 
ductivity  and  higher  local  power  deposition. 
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FIG.  13.  (Color)  Plasma  parameters  for  pressures  from  3  to  50  Ton 
(He/0:  =  70/30.  25  MHz  rf  excitation.  or= 0.3  eV/ molecule),  (a)  Power 
deposition  and  (b)  electron  density.  All  values  are  averaged  over  one  rf 
period.  The  flow  is  from  the  left.  The  scales  are  linear  with  zero  minimum 
values  unless  the  number  of  decades  is  indicated  for  log  plots.  The  maxi¬ 
mum  value  is  indicated  in  each  figure.  Discharge  constriction  occurs  near 
the  downstream  electrode  due  to  rarefaction  of  the  gas. 


A  measure  of  the  onset  of  a  constriction  instability  is  the 
maximum  value  of  the  pressure  normalized  power  deposi¬ 
tion,  P  (W/cm3  Torr).  When  keeping  a  constant  and  increas¬ 
ing  pressure,  and  in  the  absence  of  constrictions.  P  should 
also  remain  a  constant.  The  maximum  value  of  P  in  the 
discharge,  shown  in  Fig.  15.  has  a  gradual  increase  from 
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FIG.  14.  (Color)  Total  gas  density  for  pressures  from  3  to  50  Torr 
(He/0:  =  7()/3().  25  MHz  rf  excitation,  nr=0.3  eV/molecule).  Significant 
rarefaction  of  the  gas  occurs  near  the  downstream  electrode. 
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PIG  15.  Maximum  value  of  the  pressure  normalized  power  deposition 
(W/em3Torr)  as  a  function  of  pressure  (He/02  =  70/30.  25  MHz  rf  excita¬ 
tion,  rr=0.3  cV/ molecule).  A  transition  from  a  uniform  and  homogeneous 
discharge  to  constricted  discharge  occurs  between  40  and  50  Toit. 

3  to  40  Torr,  folic  wed  by  a  rapid  increase  above  40  Torr. 
This  rapid  increase  indicates  a  constriction  in  the  discharge. 

The  axial  distributions  of  gas  temperature  and  neutral 
oxygen  species  are  shown  in  Figs.  16  and  17  for  pressures 
from  10  to  50  Ton  for  c*=0.3  and  0.9  eV/molecule  (133  and 


FIG.  16.  Axial  distribution  of  (a)  gas  temperature,  (b)  density,  (c) 
02( 'A)  +  02(  '~)  density,  and  (c)  combined  02( 'AJ+O^'S)  yields  (He/02 
=70/30.  25  MHz  rf  excitation)  for  o=0.3  eV/moleeulc. 


400  W).  The  higher  gas  temperatures  obtained  at  high  pres¬ 
sures  result  in  depletion  of  ozone  in  the  afterglow.  The  maxi¬ 
mum  yields  (15%-20%  for  a=0.3  and  25%-30%  for  a 
=0.9)  are  obtained  at  50  Torr  at  the  edge  of  the  discharge 
zone.  The  yields  then  drop  below  5%  for  cr=0.3  and  below 
10%  for  a=0.9  due  to  quenching. 

The  hierarchy  of  the  main  quenchers  of  O^'A)  for  10, 
30.  and  50  Torr  obtained  with  the  2D  model  is  shown  in 
Fig.  18  (c*=0.3  cV/molcculc)  and  Fig.  19  (£*=0.9  eV/ 
molecule).  For  the  lower  energy  input  (0.3  cV/ molecule).  Oj 
is  the  main  quencher  in  the  afterglow  up  to  30-40  Torr.  This 
result  is  basically  the  same  as  predicted  by  the  plug-flow 
modeling.  With  increasing  pressure  and  increasing  absolute 
densities  of  atomic  oxygen  and  02(IA),  quenching  by  these 
species  begins  to  dominate.  The  pooling  reaction. 

02(,A)  +  02(,A)-»02(,5:)  +  02,  (22) 

generates  02( 1 S)  as  a  product  that  is  quenched  to  form 
02(!A).  For  the  higher  energy  input  (0.9  eV/molecule). 
atomic  oxygen  and  CM1  A)  become  the  main  quenchers  in 


FIG.  17.  Axial  distribution  of  (a)  gas  temperature,  (b)  density,  (c) 
02(lA)+0:(lii)  density,  and  (c)  combined  02(  'A)  +  02( '!)  yields  (He/Or 
=  70/30.  25  MHz  if  excitation)  for  0=0.9  eV/molecule. 
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FIG.  18.  Quenching  limes  for  Oi^A)  for  (a)  10.  (b)  30.  and  (c)  50  Torr 
(He/02  =  70/30,  25  MHz  rf  excitation)  for  /r=0.3  eV/molecule.  Up  to 
30-40  Torr  ozone  is  the  main  quencher  in  the  afterglow.  With  increasing 
pressure  and  increasing  absolute  densities  of  atomic  oxygen  and  pooling 
reactions  of  OT'A).  quenching  by  these  species  begins  to  dominate. 

the  afterglow  even  for  lower  pressures,  as  shown  in  Fig.  19. 

A  summary  of  pressure  scaling  for  O^A)  production  is 
shown  in  Fig.  20.  02(JA)  densities  and  yields  obtained  from 
the  2D  model  as  a  function  of  power  (for  a=0.3,  0.6,  and 
0.9  eV/ molecule)  are  plotted  for  different  pressures.  These 
values  arc  taken  on  the  tube  axis  near  the  exit  plane  at 
58  cm.  Absolute  CL^A)  production  generally  increases  with 
power  for  any  given  pressure;  however,  it  is  not  necessary 
maximum  at  the  highest  pressure.  The  yield  of  02C&)  also 
increases  with  power  for  any  given  pressure;  however,  the 
yield  at  a  given  power  generally  decreases  with  pressure. 
These  trends  are  largely  a  result  of  the  increase  in  gas  heat¬ 
ing  and  influence  of  03,  O,  and  CLf'A)  quenching  at  higher 
pressures,  and  discharge  constriction  resulting  from  nonuni¬ 
form  rarefaction  along  the  flow  tube. 

V.  CONCLUDING  REMARKS 

Pressure  scaling  of  the  absolute  densities  and  yields  of 
02(‘A)  in  flowing  He/02  plasmas  was  investigated  using 
plug-flow  and  2D  plasma  hydrodynamics  models.  We  found 
that  the  densities  and  yields  of  02(IA)  can  have  significant 
sublinear  scaling  with  pressure.  Although  yields  may  de¬ 
crease  with  incrcas  ng  pressure,  absolute  densities  of  02('A) 
typically  do  increase.  Although  these  results  depend  on  the 
layout  of  the  electrodes  and  the  aspect  ratio  of  the  flow  tube 
more  general  conclusions  can  be  made.  Obtaining  high 
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FIG.  19.  Quenching  times  for  O.f'A)  for  (a)  10,  (b)  30.  and  (c)  50  Torr 
(He/02  =  70/30,  25  MHz  rf  excitation)  for  «=0.9  eV/molecule  Quenching 
by  atomic  oxygen  and  pooling  reactions  of  02('A)  dominate 

yields  of  02(,A)  will  require  careful  management  of  the  03 
density.  Left  unchecked,  quenching  of  02( 1  A)  by  the  03  pro¬ 
duced  at  high  pressure  is  a  rate  limiting  step.  At  higher  en¬ 
ergy  densities  and  pressures,  quenching  by  O  atoms  and  cn- 


FIG.  20.  Summary  of  pressure  scaling  of  (a)  0:( 1  A)  densities  and  (b)  yields 
as  a  function  of  power  (He/0:  =  70/30.  25  MHz  rf  excitation).  Values  are  on 
the  axis  at  58  cm.  CKl'A)  production  increases  with  increasing  eV/molecule 
while  yields  at  higher  pressures  may  decrease. 
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ergy  pooling  with  02(,A)  may  also  become  important.  This 
is  particularly  the  case  for  possible  three-body  quenching 
involving  O  atoms,  though  the  importance  of  this  reaction 
can  be  minimized  by  managing  the  density  of  O  atoms  with 
additives  such  as  NO.  Pressure  scaling  also  requires  manage¬ 
ment  of  the  gas  temperature,  as  the  rate  of  exothermic  re¬ 
combination  react  ons  rapidly  increase.  Although  the  de¬ 
struction  of  O3  at  higher  temperatures  is  beneficial,  in 
general  an  intermediate  gas  temperature  is  likely  preferred. 
Discharge  stability  must  also  be  managed  and  this  will  be 
highly  temperature  dependent.  Nonuniform  excursions  of  gas 
temperature  in  the  plasma  zone  will  produce  nonuniform 
power  deposition  and  constrictions  near  electrodes.  As  such, 
low  aspect  ratio,  transverse  electrodes  allowing  for  more  ag¬ 
gressive  gas  cooling  and  more  distributed  electric  fields  will 
likely  be  required  or  high  pressure  operation. 
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